Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,815,859	Δ82 bp		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3815859	3815940	82	54 [0]	[0] 54	[rph]–[rph]	[rph],[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	0 (0.000)	51 (0.990)	41/260	0.2	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
	?	NC_000913	3815941 =	0 (0.000)					pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815720‑3816000                                                                                                                                           |                                                                                                                                               ccTGCATGCTCAGCGCAAATGCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTCCTTTTCTACAGACCAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                                                            <  2:121425‑M1/143‑5 (MQ=255)   tgcttggtAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                          <  1:754770‑M1/143‑7 (MQ=255)   tgcttgccAAGCGCAAATGCAAGAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                          <  2:286083‑M1/143‑7 (MQ=255)   tccttgctaaCCCCAAATTCAATAAACTGCCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                          <  1:538554‑M1/141‑7 (MQ=255)          tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACGAACAGGGGCCGCAACAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                                                                                                   >  1:468463‑M1/1‑130 (MQ=255)           aaGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGTTCCTCCTCTTACCTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAATACCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                                                                                  <  2:496817‑M1/143‑15 (MQ=255)                  acttccATCAAGTGGCGCGGATATGGTTTCGTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                           <  2:977156‑M1/141‑22 (MQ=255)                  aaTCCAATAGCCTCGCGCTGACATGGTTCCATGCCTTCCCTCCTCATCTTACTCTTCTACAGACAAACAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                           <  1:286294‑M1/143‑22 (MQ=255)                   aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                          >  1:251805‑M1/1‑121 (MQ=255)                       aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                      >  1:399128‑M1/1‑117 (MQ=255)                        aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCTTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCACGGCGTTCTGCcctcacctgcacttcactgctgctccc                                                                                                                     >  2:746108‑M1/1‑116 (MQ=255)                           acctgGCGCTGATATCTTCTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                                                                                                  <  1:133817‑M1/141‑31 (MQ=255)                           aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                                                                                                  >  2:320471‑M1/1‑113 (MQ=255)                            aCTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTCATACAGTCAAAAAAAAGGCGACTCATGAGTCGCCTTAAAAATCATTTTGCCAGAGGGGCGTTGTTTcgtcccctgcacttcaatgatgcgcccgtct                                                                                                                 <  2:613835‑M1/143‑32 (MQ=255)                              tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCACCGCCTTCTCCcgtcccctgcacttcgctgatcccgccgccttc                                                                                                               >  1:723889‑M1/1‑110 (MQ=255)                               ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  2:701477‑M1/143‑35 (MQ=255)                               ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:819170‑M1/143‑35 (MQ=255)                               ggCGCTGATATGGTTTCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAGAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:468491‑M1/143‑35 (MQ=255)                                 cgcTGATATGGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                            <  1:946509‑M1/143‑37 (MQ=255)                                    tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                         <  1:588786‑M1/143‑40 (MQ=255)                                            tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgaggcccccgccttcgctcatcaccaagt                                                                                                 >  2:343793‑M1/1‑96 (MQ=255)                                              tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                               >  2:806109‑M1/1‑94 (MQ=255)                                                    cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccttcacttcaatgatgcgcccgtcttcggtcatcactacgttcatatcc                                                                                         >  2:138721‑M1/1‑88 (MQ=255)                                                            ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACACATCCGCTGATGAGTGTCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                                 <  1:232572‑M1/143‑64 (MQ=255)                                                                     cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                                                        <  2:533454‑M1/143‑73 (MQ=255)                                                                      ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTCAAAATCAGATTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcaccactacgttcatgtccatctctcccgcacagacc                                                                       >  2:385865‑M1/1‑70 (MQ=255)                                                                      ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                                                       <  2:251805‑M1/143‑74 (MQ=255)                                                                           tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAGAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagtcttcaac                                                                  >  1:653602‑M1/1‑65 (MQ=255)                                                                           tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                                  >  1:531255‑M1/1‑65 (MQ=255)                                                                              aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcctcaacgta                                                               >  2:579945‑M1/1‑62 (MQ=255)                                                                                  aaaaaaaaGGCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccaaaccgctcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                           <  2:492704‑M1/143‑86 (MQ=255)                                                                                   aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcca                                                          <  1:507054‑M1/143‑87 (MQ=255)                                                                                     aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtgttccaga                                                        >  1:305827‑M1/1‑55 (MQ=255)                                                                                     aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga                                                        <  2:95432‑M1/143‑89 (MQ=255)                                                                                      ttaaggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                       >  1:295653‑M1/3‑54 (MQ=255)                                                                                      aaaaGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                       <  1:686377‑M1/143‑90 (MQ=255)                                                                                             gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac                                                <  1:806109‑M1/143‑97 (MQ=255)                                                                                              aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagccttcaacgtattccagatcgcaaacc                                               >  2:580866‑M1/1‑46 (MQ=255)                                                                                                    caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttctactggttcctgatcgcacacggcgtcg                                         >  1:536519‑M1/1‑40 (MQ=255)                                                                                                           cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                                  >  2:570269‑M1/1‑33 (MQ=255)                                                                                                              aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtgttccagatcgcaaaccgcttcgccgttcacaa                               <  2:531255‑M1/143‑114 (MQ=255)                                                                                                              aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgcccttcacca                               >  2:703624‑M1/1‑30 (MQ=255)                                                                                                                 aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaaatcgcaaaccgcttcgccgttcacaattc                            >  1:410874‑M1/1‑27 (MQ=255)                                                                                                                  ttcggtttgCCCGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                           <  2:987951‑M1/139‑118 (MQ=255)                                                                                                                  aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                           <  1:574432‑M1/143‑118 (MQ=255)                                                                                                                      gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaca                       >  1:176670‑M1/1‑22 (MQ=255)                                                                                                                      gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagccttcaacctattccaaaccgcaaatcgcttcgcccttcacccttcccaca                       >  2:803402‑M1/1‑22 (MQ=255)                                                                                                                          gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaactccgccagcac                   >  2:1361‑M1/1‑18 (MQ=255)                                                                                                                            cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaact                 >  1:250871‑M1/1‑16 (MQ=255)                                                                                                                            cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcaccgttcacaattccgacagaaact                 >  1:751792‑M1/1‑16 (MQ=255)                                                                                                                                gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg             <  2:250871‑M1/143‑132 (MQ=255)                                                                                                                                       tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccat      >  1:893309‑M1/1‑5 (MQ=255)                                                                                                                                           ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccct  >  1:26200‑M1/1‑1 (MQ=255)                                                                                                                                           ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagaggctccaacgtattccagatcgcaaaccgcttcgccgttcacgatgccgccagggaccgcgcccaccatcctt  >  1:213472‑M1/1‑1 (MQ=255)                                                                                                                                           |                                                                                                                                               CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815720‑3816000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑