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	33 [0]	[0] 33	[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)	33 (0.880)	32/264	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

CTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGC  >  NC_000913/3815814‑3816072                                                                                                                                |                                                                                                                                    gcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagccgcctcaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCa                                                                                                                      >  1:724207‑M2/128‑143 (MQ=255)       tttcaataacccggcgctgatatggtttcatgccttcgctccgcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                <  1:643036‑M2/22‑1 (MQ=255)        ttcaataaactggcgctgatatggtttcatgccttcgctcatgatcttacttttgtagagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgc                                                                                                               >  1:124354‑M2/121‑143 (MQ=255)         tcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagccaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATCATgcg                                                                                                              >  1:493297‑M2/120‑143 (MQ=255)            ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                           <  2:264459‑M2/27‑1 (MQ=255)              aaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGGACTTCAGTGGTGCCCCCGt                                                                                                         >  1:391160‑M2/115‑143 (MQ=255)                                  tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                     <  1:590644‑M2/49‑1 (MQ=255)                                      gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGAAGTTCTACTTTTGGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                 >  1:495704‑M2/91‑143 (MQ=255)                                         ttcgctcctcatcttacttttctacagacaaaaaaaatgctactcatcagtcgacttaaaaatcagtttgccagcgccgcctactacCGTCCCCGGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCATGTCgg                                                                              >  2:484870‑M2/88‑143 (MQ=255)                                                 tcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcaatttgccaccgccgccttctgcCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGCCATCCCTACGTTCACGCCGGTGCCTGCg                                                                      >  2:234165‑M2/80‑143 (MQ=255)                                                  catcttactgttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                     <  2:395361‑M2/65‑1 (MQ=255)                                                        acttttctacagacaaaaaaaaggcgacgcatcagtcgccttaaaactcattttgcctgcggcgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                               <  2:352320‑M2/71‑1 (MQ=255)                                                         cttatgaacagagaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCTATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                              <  1:232792‑M2/72‑1 (MQ=255)                                                          ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                             >  1:396083‑M2/71‑143 (MQ=255)                                                            ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGGCGGTCTCTGCGGCAGAGTCTTc                                                           >  1:220643‑M2/69‑143 (MQ=255)                                                                acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                       >  1:417489‑M2/65‑143 (MQ=255)                                                                  agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagttggccagcaccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                     >  1:543122‑M2/63‑143 (MQ=255)                                                                   gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTAt                                                    >  1:226382‑M2/62‑143 (MQ=255)                                                                      cagaaaaaggcgactcatcagtcgccttaaaagtcagtttgccagcgccgccttctgcCTTGCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                 <  2:634293‑M2/85‑1 (MQ=255)                                                                         aaaaatgcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCCGa                                              >  2:560173‑M2/56‑143 (MQ=255)                                                                          aaaaggctactcatcagtagccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                             <  1:492388‑M2/89‑1 (MQ=255)                                                                           agaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTTCACTTCTATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                            <  2:634682‑M2/90‑1 (MQ=255)                                                                           aaaggcgcctcctcagtcgccttaaagatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGGTTCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                            <  2:101056‑M2/90‑1 (MQ=255)                                                                                    tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTATGTGCGGCTGAGTCTTCAACGTATTCCAGATCGCCAACcgt                                   >  1:495683‑M2/45‑142 (MQ=255)                                                                                      atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCAATCAACAGACATCGCTACGCCCACCTCGTTCTTTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                 >  1:200573‑M2/43‑143 (MQ=255)                                                                                       tcagtcgccttaaaaatcagtttgccagcgccgccttcagcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                <  1:54499‑M2/102‑1 (MQ=255)                                                                                           tcgccttaaacatccgtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                            <  2:561969‑M2/106‑1 (MQ=255)                                                                                                   aagaccagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAt                    <  1:437739‑M2/114‑1 (MQ=255)                                                                                                     aatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCCGATCGCAAACCGCTCCGCCCTTCCCAATTc                  >  2:166095‑M2/28‑143 (MQ=255)                                                                                                       ccagttcgccagcgccgccttctgcCGTCCCCTGCACTTCAACGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCGGTCTCTGCGGCACAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                <  2:462052‑M2/118‑1 (MQ=255)                                                                                                         agtttgcctgcgccgccttctgcCGTCCCCTGCACGTCAATGGTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc              <  2:727005‑M2/120‑1 (MQ=255)                                                                                                                 agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGAAACTg      >  1:506522‑M2/16‑143 (MQ=255)                                                                                                                     ccgccttctgcCCTCCCCTGCACTTCATTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGc  <  2:194192‑M2/132‑1 (MQ=255)                                                                                                                                |                                                                                                                                    CTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGC  >  NC_000913/3815814‑3816072

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