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	44 [0]	[0] 43	[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)	39 (0.830)	36/262	0.3	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

AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815801‑3816080                                                                                                                                             |                                                                                                                                            acctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttaattttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgacagcgccgccttctgcCGt                                                                                                                                           >  2:280442‑M2/141‑143 (MQ=255)  cctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                          <  2:684938‑M2/4‑1 (MQ=255)                         ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                   <  1:829654‑M2/27‑1 (MQ=255)                         ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                   >  1:741915‑M2/117‑143 (MQ=255)                          taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggagacagatcagtagccttaaaaataagtttgccagggcagccttctgcAGTCCCCTGCACTTCAATGATGCGACCg                                                                                                                  >  1:88677‑M2/116‑143 (MQ=255)                                ggtgctgatatgctttcgtgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                            <  2:215817‑M2/34‑1 (MQ=255)                                 gcgctgatatggtttcatgccttagctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGCCAGCCTTCgg                                                                                                           >  2:924230‑M2/109‑143 (MQ=255)                                   gctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatccgtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACATACATGCTGCGCCCGCCCTCGgtg                                                                                                         >  2:267292‑M2/107‑142 (MQ=255)                                              atactgcctcctctcctcatcttccttttcgacacacaaaaaaaacgcgactcatcattcgcctgaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                              <  2:799388‑M2/48‑1 (MQ=255)                                               tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                             <  1:99171‑M2/49‑1 (MQ=255)                                                  tgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCCGTCCTCCCTACGTtcacg                                                                                          >  2:923525‑M2/92‑141 (MQ=255)                                                    ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTCCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGTc                                                                                        >  1:684869‑M2/90‑143 (MQ=255)                                                     cttctctcctcatctgacttcgctacagacgaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                       <  2:872412‑M2/55‑1 (MQ=255)                                                        cgctcctcatcttacttttctacagacaaaaaaaaggcgactcataagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                                                                    >  1:840227‑M2/86‑143 (MQ=255)                                                            cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgtGGTGTGCCAGACGCCAATGATGCGCCCGTCTTCGGTAATCACTACGTTAATGTCGGGCTATg                                                                                >  1:428436‑M2/82‑143 (MQ=255)                                                             ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATCTCGGTCTCTGc                                                                               >  2:203927‑M2/81‑143 (MQ=255)                                                                 tcttacttttaaacggaaagaaaaaaggcgactcatcagtcgccttaaaaatcggtttgccagcgccgccttctgcCGTCCCCCGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                           <  1:326653‑M2/67‑1 (MQ=255)                                                                  ctcccttttctacagacaaaaaaaaggcgactcaccagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                          <  1:293500‑M2/68‑1 (MQ=255)                                                                   tactgttctacagacaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                         <  1:571508‑M2/69‑1 (MQ=255)                                                                     acttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGGACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGt                                                                       >  2:108322‑M2/73‑143 (MQ=255)                                                                          tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                                                                  <  2:840227‑M2/76‑1 (MQ=255)                                                                               agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCTGTCTCTGCGGCAGAGTCTTCAACGTa                                                             >  1:20251‑M2/63‑143 (MQ=255)                                                                               agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccctctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCGTGGCGGCCACGGCCGCAGGGTCTTCAccgga                                                             >  2:144810‑M2/63‑141 (MQ=255)                                                                                        aaaggcgtctcatcagtcgccttaaaaatcagtttgccagcgccgctttctgtCGATCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                    <  1:834130‑M2/90‑1 (MQ=255)                                                                                                    tcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                        >  1:69304‑M2/42‑143 (MQ=255)                                                                                                       gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGcc                                     >  1:75232‑M2/39‑143 (MQ=255)                                                                                                             ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                               >  1:917196‑M2/33‑143 (MQ=255)                                                                                                             tgaaaaaccaggttgccgtcgccgccttctgcCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                               <  2:774733‑M2/111‑1 (MQ=255)                                                                                                             cgaaaaatccgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                               <  2:592785‑M2/111‑1 (MQ=255)                                                                                                              taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTGAGCCGTTccca                              >  2:140278‑M2/32‑143 (MQ=255)                                                                                                                 aaatccgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCAGAGTTTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtt                           <  1:535296‑M2/115‑1 (MQ=255)                                                                                                                  aatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                          <  1:828453‑M2/116‑1 (MQ=255)                                                                                                                     cagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                       <  1:564674‑M2/119‑1 (MQ=255)                                                                                                                      agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATAGCAAACCGCTTCGCCGTTCTCAATTACGAc                      >  1:47039‑M2/24‑143 (MQ=255)                                                                                                                      agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGTCATCACTACGTTCAGGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                      <  1:374499‑M2/120‑1 (MQ=255)                                                                                                                       gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATCCCAGATCTCAAACCGCTTCGCCGTTCACAATTCCGTCa                     >  1:368057‑M2/23‑143 (MQ=255)                                                                                                                        tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAg                    >  1:736888‑M2/22‑143 (MQ=255)                                                                                                                        tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGCATTCCAGATCGCAAACCGCCTCGCCGTTCCCAACTCCgaccg                    >  2:234604‑M2/22‑141 (MQ=255)                                                                                                                               gcgccgacttctgcCGTCCCATGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCTAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc             <  2:238404‑M2/129‑1 (MQ=255)                                                                                                                                   cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGACAGAACCTGCGGCt         >  2:587790‑M2/11‑143 (MQ=255)                                                                                                                                     ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGCATTCCAGATCGCAAACCGCTTCGCCGTTGCCAATTCCGCCAGAACATGCGgctgc       >  1:302752‑M2/9‑141 (MQ=255)                                                                                                                                       ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa     >  1:453660‑M2/7‑143 (MQ=255)                                                                                                                                          tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc  <  2:741915‑M2/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                            AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815801‑3816080

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