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

AAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815805‑3816083                                                                                                                                         |                                                                                                                                               gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                                                                          >  1:801842‑M2/137‑143 (MQ=255)  cttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcccccttctgcCGTCCCCt                                                                                                                                         >  1:479183‑M2/136‑143 (MQ=255)     gctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggccactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                      >  2:179432‑M2/133‑143 (MQ=255)             caactgtaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                              >  2:148460‑M2/125‑143 (MQ=255)             caaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                              >  1:212036‑M2/125‑143 (MQ=255)                attcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                           >  2:392147‑M2/122‑143 (MQ=255)                           tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGCCCGTCTTc                                                                                                                >  1:785981‑M2/111‑143 (MQ=255)                            ggcgctgatctggtttcatgccttcgcccctccgcttactcttctacagacgaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  2:538983‑M2/34‑1 (MQ=255)                            ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  1:245648‑M2/34‑1 (MQ=255)                              cgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                             <  1:561828‑M2/36‑1 (MQ=255)                                          ttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                 <  2:316564‑M2/48‑1 (MQ=255)                                                       tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                    >  1:57512‑M2/83‑143 (MQ=255)                                                        cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                   >  1:174695‑M2/82‑143 (MQ=255)                                                             tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTTTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                              <  2:703710‑M2/67‑1 (MQ=255)                                                             tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                              >  1:129265‑M2/77‑143 (MQ=255)                                                              cttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTAGTATCGGTCTCTGCGGCag                                                                             >  2:133783‑M2/76‑143 (MQ=255)                                                               ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                            >  1:707563‑M2/75‑143 (MQ=255)                                                                       ctacagacaaaacaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                    <  2:819348‑M2/77‑1 (MQ=255)                                                                         acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                                  >  1:688422‑M2/65‑143 (MQ=255)                                                                         acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCCCTACGGCAGAGTCTTCAACg                                                                  >  2:528983‑M2/65‑143 (MQ=255)                                                                              caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                             >  1:495392‑M2/60‑143 (MQ=255)                                                                              caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCCCTACGTTCATGTCGGTCACTGCGGCGGAGCCTTCGACGTAGTc                                                             >  1:221069‑M2/60‑143 (MQ=255)                                                                                   ttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGAt                                                        >  1:455674‑M2/55‑143 (MQ=255)                                                                                          gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCACTGCACTTCAATGATGCGCCCGTCTTCAGTGATCACTAAGATCATGTCGGTCTCTGCGGCAGAGTCTTCAAAGTATTCCAGATCccaaga                                                 >  2:100175‑M2/48‑141 (MQ=255)                                                                                             tcatcagtcgccttaaaaatcagtttgccagcgccgccttgtgcCGTCGCCTGGACTTCAATGATGCGCCAGACTACGGTCATCCCTACGTTCATGGCGGACTCTGCGGGAGCGTCTGCAACGTATTCCAGGGCGCAAACCGc                                              >  1:177726‑M2/45‑143 (MQ=255)                                                                                               atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                            >  1:680811‑M2/43‑143 (MQ=255)                                                                                               atcagcctccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                            <  1:833332‑M2/101‑1 (MQ=255)                                                                                                  agtcgccttaaaactcagtctgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGc                                         <  1:30477‑M2/104‑1 (MQ=255)                                                                                                      gccttaaaaatcattttgccagcgccgctttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                     <  2:38552‑M2/108‑1 (MQ=255)                                                                                                       ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                    >  1:727125‑M2/35‑143 (MQ=255)                                                                                                        cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTACTGTCGGTCTCTCCGGCAGCGACGTCAACGCATTCCAGATCGCCAACCGCCTCCCCGCTca                                   >  2:697893‑M2/34‑143 (MQ=255)                                                                                                              attcggtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                             <  2:243536‑M2/116‑1 (MQ=255)                                                                                                              aatcagtttgccagcgccgccttctgcCGTCCCCTACACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGCCGGCCTCTGCAGCAGAGTCTTCAACGTATTCCAGATCCCAAACAGCTTCGCCGTTCCCAGTTc                             >  2:344814‑M2/28‑143 (MQ=255)                                                                                                                  agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGACGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                         <  1:231081‑M2/120‑1 (MQ=255)                                                                                                                  actttgcaggcgccgccgtctgcGGCCCCCTGCCCTTAACTTGTGCGCCCGTCTTCGGTCATCACTACGTTCAGGTCGGCCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                         <  2:810230‑M2/120‑1 (MQ=255)                                                                                                                          agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCACGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAGAGAAACTg                 >  1:440475‑M2/16‑143 (MQ=255)                                                                                                                             gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg              >  1:128979‑M2/13‑143 (MQ=255)                                                                                                                                      tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCAGATCGCAAACCGCTTAGCCCTTCACAATGCCTACCGAACCGCCCGCTACCATcc     >  2:490960‑M2/4‑143 (MQ=255)                                                                                                                                        cCGTCCCCTGCACTTCAATGATGCGCCCGTGTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt   <  2:311302‑M2/142‑1 (MQ=255)                                                                                                                                         cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt  <  1:696939/143‑1 (MQ=255)                                                                                                                                         |                                                                                                                                               AAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815805‑3816083

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: ‑