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	42 [0]	[0] 42	[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)	41 (1.190)	35/264	0.1	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

ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAT  >  NC_000913/3815799‑3816078                                                                                                                                               |                                                                                                                                          ccacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcacgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                           <  1:530746‑M2/1‑1 (MQ=255) acacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                           <  1:141101‑M2/1‑1 (MQ=255)    cctgcttgctaagcgcaaattcaagaaactggcgctgatatggtttcacgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                        <  1:544325‑M2/4‑1 (MQ=255)      tgcttgctaagcgcaaattcaatacactggcgctgatatggtttcatgcattcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcccccttctgcCGTccc                                                                                                                                      >  1:590363‑M2/138‑143 (MQ=255)      tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcactttgccagcgccgccttctgcCGTccc                                                                                                                                      >  1:448629‑M2/138‑143 (MQ=255)                  gacaattcaacagactggcgctgatatggtgtcatgccttcgctgctcatcttactttgctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                          <  2:170078‑M2/18‑1 (MQ=255)                    aaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgcattccgcCGTCACCTGCACTCCAatga                                                                                                                        >  2:501882‑M2/124‑143 (MQ=255)                        tcaataaagtggcgctgatatggtttcatgccttcgctcctcatcttgcttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcg                                                                                                                    <  2:183305‑M2/24‑1 (MQ=255)                             caactcgccctgatatggtttcatgccttcgctcctcatcttacttttctacacacaaataaaaggcgactcaccagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                                                               <  2:340727‑M2/29‑1 (MQ=255)                             aaactggcgctgatatggtttcatgccttcgctcctcaccttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                                                               <  2:373699‑M2/29‑1 (MQ=255)                                 tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgcctcatcagccgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTCCACTTCAATGATGCGCCCGTCTTc                                                                                                           >  1:162584‑M2/111‑143 (MQ=255)                                        gatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcccccttctgcCGTCCCCCGCCCCTCCATGCTGCACCCGtcgtcggcactg                                                                                                    >  1:512126‑M2/104‑138 (MQ=255)                                        gatatgattccacgtcttcgctcctcaccttacttttctacagacacaaaaaaggcgactcatcagtcgccttaaaaatcactttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                    <  1:245195‑M2/40‑1 (MQ=255)                                           atggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGACGCGCCCGCCCCCGCTCAGCACt                                                                                                 >  2:540741‑M2/101‑143 (MQ=255)                                            tggtttcatgccttcgctcctcatcttacgcgtctacagacaaaaaagacgcgactcatcagtcgccttaaaaatcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                <  2:264628‑M2/44‑1 (MQ=255)                                               tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                             >  1:630228‑M2/97‑143 (MQ=255)                                                        ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATATCgg                                                                                    >  2:827‑M2/88‑143 (MQ=255)                                                                   tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAGTGATGCCCCCGCCTCCGGCCCTCACTACGTTCATGTGGCGCCCTGCGGCa                                                                         >  2:580961‑M2/77‑143 (MQ=255)                                                                     ttactgttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                       <  2:571818‑M2/69‑1 (MQ=255)                                                                             ctacagacaaaaaaaaggcgactcagcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTTCCCTGCACTTCAATGATGCGCCCGTCGTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                               <  1:488570‑M2/77‑1 (MQ=255)                                                                              ttaagacagaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                              <  2:67181‑M2/78‑1 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                              <  2:179268‑M2/78‑1 (MQ=255)                                                                               acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTCCGTCCGTGTCGGTCTCTGCGGCAGACTCTTCAACg                                                             >  1:294045‑M2/65‑143 (MQ=255)                                                                                 agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCCTGATGCGCCCGTCTTCGGTCTTCCCTACGTTCATGTCGGTCACTACAGCAGAATATTCACCGTa                                                           >  1:520714‑M2/63‑143 (MQ=255)                                                                                 agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                           <  1:368458‑M2/81‑1 (MQ=255)                                                                                   acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGAACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGCCTACCCGCCTGAGGCTTcaacgactt                                                         >  2:261947‑M2/61‑139 (MQ=255)                                                                                     aaaaaaaaggcgactcatcagtcgccttaaaaatcagattgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCGTCACTACGTTCATGTCGCTCTCTGCGGCAGAGACATCAACGTATTcc                                                       >  2:321331‑M2/59‑143 (MQ=255)                                                                                        aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCATCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACCATCATGACGTTGTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                    <  1:641466‑M2/88‑1 (MQ=255)                                                                                           aaggcgactcatcagtcgccttaaaaatcagtttgccagcgcccccttccgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg                                                 >  1:28522‑M2/53‑143 (MQ=255)                                                                                                  ctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAGTCCAGATCGCAAACCg                                          >  2:124630‑M2/46‑143 (MQ=255)                                                                                                     atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                       >  1:361203‑M2/43‑143 (MQ=255)                                                                                                     atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATAGCAAACCGctg                                       >  1:361209‑M2/43‑142 (MQ=255)                                                                                                     atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCATGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                       >  2:348953‑M2/43‑143 (MQ=255)                                                                                                      tcagtcgtgttaaaaatcagtgtgccagcggcgccttctgcCGTCCCCTGCACTTCAGTGATGCGCCCGTCTTCGGTCATCATTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                      <  2:176655‑M2/102‑1 (MQ=255)                                                                                                          tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCCTCAACGTCTTCCAGATCGCAGACCGCTTCgccc                                  >  1:44258‑M2/38‑142 (MQ=255)                                                                                                              cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTCTTCCAGCTCGCACACCGCTTCGCCTTTca                              >  1:542207‑M2/34‑143 (MQ=255)                                                                                                               ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                             >  1:477569‑M2/33‑143 (MQ=255)                                                                                                                     atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                       >  1:481144‑M2/27‑143 (MQ=255)                                                                                                                          tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCAcacttcctccag                  >  1:578193‑M2/22‑138 (MQ=255)                                                                                                                           ttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCCGATCGCGAACCGCTTCGCCGTTCCCACTTCCGccaca                 >  1:165433‑M2/21‑141 (MQ=255)                                                                                                                             gccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACagaga               >  1:269685‑M2/19‑141 (MQ=255)                                                                                                                                          tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGACAGAAACTGCGGCTACCAt  >  1:400691‑M2/6‑143 (MQ=255)                                                                                                                                               |                                                                                                                                          ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAT  >  NC_000913/3815799‑3816078

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