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	39 [0]	[0] 40	[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)	37 (0.840)	31/264	0.4	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

GCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815723‑3816000                                                                                                                                        |                                                                                                                                               gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                                                         >  1:801842‑M1/1‑136 (MQ=255)  cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCcgtcccct                                                                                                                                        >  1:479183‑M1/1‑135 (MQ=255)     gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCCACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                                                     >  2:179432‑M1/1‑132 (MQ=255)             caactgtaataaactgGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                             >  2:148460‑M1/8‑124 (MQ=255)             cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                             >  1:212036‑M1/1‑124 (MQ=255)                aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                          >  2:392147‑M1/1‑121 (MQ=255)                           tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgatgcgcccgtcttc                                                                                                               >  1:785981‑M1/1‑110 (MQ=255)                            ggCGCTGATCTGGTTTCATGCCTTCGCCCCTCCGCTTACTCTTCTACAGACGAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  2:538983‑M1/143‑35 (MQ=255)                            ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:245648‑M1/143‑35 (MQ=255)                              cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                            <  1:561828‑M1/143‑37 (MQ=255)                                          ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                                <  2:316564‑M1/143‑49 (MQ=255)                                                       tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctct                                                                                   >  1:57512‑M1/1‑82 (MQ=255)                                                        cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                                                                  >  1:174695‑M1/1‑81 (MQ=255)                                                             tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttctttcatcactacgttcatgtcggtctctgcggca                                                                             <  2:703710‑M1/143‑68 (MQ=255)                                                             tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                             >  1:129265‑M1/1‑76 (MQ=255)                                                              cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttagtatcggtctctgcggcag                                                                            >  2:133783‑M1/1‑75 (MQ=255)                                                               ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                           >  1:707563‑M1/1‑74 (MQ=255)                                                                       cTACAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                                                                   <  2:819348‑M1/143‑78 (MQ=255)                                                                         aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                                                 >  1:688422‑M1/1‑64 (MQ=255)                                                                         aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtccctacggcagagtcttcaacg                                                                 >  2:528983‑M1/1‑64 (MQ=255)                                                                              cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                            >  1:495392‑M1/1‑59 (MQ=255)                                                                              cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcacctgcacttcaatgatgcgcccgtcttcggtcctccctacgttcatgtcggtcactgcggcggagccttcgacgtagtc                                                            >  1:221069‑M1/1‑59 (MQ=255)                                                                                   ttaaggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggccgagtcttcaacgtattccagat                                                       >  1:455674‑M1/3‑54 (MQ=255)                                                                                          gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccactgcacttcaatgatgcgcccgtcttcagtgatcactaagatcatgtcggtctctgcggcagagtcttcaaagtattccagatcccaaga                                                >  2:100175‑M1/1‑47 (MQ=255)                                                                                             tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTGTGCcgtcgcctggacttcaatgatgcgccagactacggtcatccctacgttcatggcggactctgcgggagcgtctgcaacgtattccagggcgcaaaccgc                                             >  1:177726‑M1/1‑44 (MQ=255)                                                                                               atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                           >  1:680811‑M1/1‑42 (MQ=255)                                                                                               atcaGCCTCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                           <  1:833332‑M1/143‑102 (MQ=255)                                                                                                  aGTCGCCTTAAAACTCAGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgc                                        <  1:30477‑M1/143‑105 (MQ=255)                                                                                                      gCCTTAAAAATCATTTTGCCAGCGCCGCTTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgtt                                    <  2:38552‑M1/143‑109 (MQ=255)                                                                                                       ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                   >  1:727125‑M1/1‑34 (MQ=255)                                                                                                        cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacggtactgtcggtctctccggcagcgacgtcaacgcattccagatcgccaaccgcctccccgctca                                  >  2:697893‑M1/1‑33 (MQ=255)                                                                                                              attcggtttgcCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                            <  2:243536‑M1/138‑117 (MQ=255)                                                                                                              aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctacactgcaatgatgcgcccgtcttcggtcatcactacgctcatgccggcctctgcagcagagtcttcaacgtattccagatcccaaacagcttcgccgttcccagttc                            >  2:344814‑M1/1‑27 (MQ=255)                                                                                                                  actttgcAGGCGCCGCCGTCTGCggccccctgcccttaacttgtgcgcccgtcttcggtcatcactacgttcaggtcggcctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                        <  2:810230‑M1/141‑121 (MQ=255)                                                                                                                  aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgacgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                        <  1:231081‑M1/143‑121 (MQ=255)                                                                                                                          aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcacgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgagagaaactg                >  1:440475‑M1/1‑15 (MQ=255)                                                                                                                             gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg             >  1:128979‑M1/1‑12 (MQ=255)                                                                                                                                      tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcgcaaaccgcttagcccttcacaatgcctaccgaaccgcccgctaccatcc    >  2:490960‑M1/1‑3 (MQ=255)                                                                                                                                        ccgtcccctgcacttcaatgatgcgcccgtgttcggtcatcactccgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccct  <  2:311302‑M1/143‑143 (MQ=255)                                                                                                                                        |                                                                                                                                               GCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815723‑3816000

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