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	54 [0]	[0] 54	[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)	51 (0.990)	41/260	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

AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815802‑3816082                                                                                                                                            |                                                                                                                                              cctgcatgctcagcgcaaatgcaataaactggcgctgatatggtttcatgccttcgctcctcctcttccttttctacagaccaaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                            <  2:121425‑M2/4‑1 (MQ=255)   tgcttggtaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  1:754770‑M2/6‑1 (MQ=255)   tgcttgccaagcgcaaatgcaagaaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  2:286083‑M2/6‑1 (MQ=255)   tccttgctaaccccaaattcaataaactgccgctgatatggtttcatgccttcgctcctcatcttacttttctacagaccaaaaaaaggcgactcatcagtcgcctgaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  1:538554‑M2/6‑1 (MQ=255)          taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaacgaacaggggccgcaacagtcgccttaaacatcagtttgccagcgccgccttctgcCGTCCCCTGCACt                                                                                                                                   >  1:468463‑M2/131‑143 (MQ=255)           aagcgcaacttcaataaactggcgctgatatggtttcatgccttcgttcctcctcttacctttctccagacaaaaaaaaggcgactcatcagtcgccttaaataccagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                  <  2:496817‑M2/14‑1 (MQ=255)                  acttccatcaagtggcgcggatatggtttcgtgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                           <  2:977156‑M2/21‑1 (MQ=255)                  aatccaatagcctcgcgctgacatggttccatgccttccctcctcatcttactcttctacagacaaacaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                           <  1:286294‑M2/21‑1 (MQ=255)                   attcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                          >  1:251805‑M2/122‑143 (MQ=255)                       aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                      >  1:399128‑M2/118‑143 (MQ=255)                        ataaactggcgctgatatggtttcatgccttcgctcttcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcacggcgttctgcCCTCACCTGCACTTCACTGCTGCTccc                                                                                                                     >  2:746108‑M2/117‑143 (MQ=255)                           acctggcgctgatatcttctcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                  <  1:133817‑M2/30‑1 (MQ=255)                           aactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                  >  2:320471‑M2/114‑143 (MQ=255)                            acttgcgctgatatggtttcatgccttcgctcctcatcttactttcatacagtcaaaaaaaaggcgactcatgagtcgccttaaaaatcattttgccagaggggcgttgtttCGTCCCCTGCACTTCAATGATGCGCCCGTCt                                                                                                                 <  2:613835‑M2/31‑1 (MQ=255)                              tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaacaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcaccgccttctccCGTCCCCTGCACTTCGCTGATCCCGCCGCCTTc                                                                                                               >  1:723889‑M2/111‑143 (MQ=255)                               ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  2:701477‑M2/34‑1 (MQ=255)                               ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:819170‑M2/34‑1 (MQ=255)                               ggcgctgatatggtttcatgcctgcgctcctcatcttacttttctacagacaaagaaaaggcgactcatcagtcgccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:468491‑M2/34‑1 (MQ=255)                                 cgctgatatggtgtcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                            <  1:946509‑M2/36‑1 (MQ=255)                                    tgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                         <  1:588786‑M2/39‑1 (MQ=255)                                            tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGAGGCCCCCGCCTTCGCTcatcaccaagt                                                                                                 >  2:343793‑M2/97‑138 (MQ=255)                                              tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                               >  2:806109‑M2/95‑143 (MQ=255)                                                    cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTTCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATAtcc                                                                                         >  2:138721‑M2/89‑142 (MQ=255)                                                            ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttacacatccgctgatgagtgtcgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                 <  1:232572‑M2/63‑1 (MQ=255)                                                                     cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                        <  2:533454‑M2/72‑1 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttcaaaatcagattgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCACCACTACGTTCATGTCCATCTCTCCcgcacagacc                                                                       >  2:385865‑M2/71‑140 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                       <  2:251805‑M2/73‑1 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttagaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAAc                                                                  >  1:653602‑M2/66‑143 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                  >  1:531255‑M2/66‑143 (MQ=255)                                                                              agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCCTCAACGTa                                                               >  2:579945‑M2/63‑143 (MQ=255)                                                                                  aaaaaaaaggcgactcatcaggcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCAAACCGCTCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                           <  2:492704‑M2/85‑1 (MQ=255)                                                                                   aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                          <  1:507054‑M2/86‑1 (MQ=255)                                                                                     aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTGTTCCAGa                                                        >  1:305827‑M2/56‑143 (MQ=255)                                                                                     aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                        <  2:95432‑M2/88‑1 (MQ=255)                                                                                      ttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                       >  1:295653‑M2/55‑143 (MQ=255)                                                                                      aaaaggcgactcatcagtcgccttaacaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                       <  1:686377‑M2/89‑1 (MQ=255)                                                                                             gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc                                                <  1:806109‑M2/96‑1 (MQ=255)                                                                                              actcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAACGTATTCCAGATCGCAAAcc                                               >  2:580866‑M2/47‑143 (MQ=255)                                                                                                    cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCTACTGGTTCCTGATCGCACACGGCGTCg                                         >  1:536519‑M2/41‑143 (MQ=255)                                                                                                           cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                                  >  2:570269‑M2/34‑143 (MQ=255)                                                                                                              aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTGTTCCAGATCGCAAACCGCTTCGCCGTTCACaa                               <  2:531255‑M2/113‑1 (MQ=255)                                                                                                              aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCCttcacca                               >  2:703624‑M2/31‑141 (MQ=255)                                                                                                                 aatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAAATCGCAAACCGCTTCGCCGTTCACAATTc                            >  1:410874‑M2/28‑143 (MQ=255)                                                                                                                  ttcggtttgcccgcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                           <  2:987951‑M2/117‑1 (MQ=255)                                                                                                                  atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                           <  1:574432‑M2/117‑1 (MQ=255)                                                                                                                      gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACa                       >  1:176670‑M2/23‑143 (MQ=255)                                                                                                                      gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAACCTATTCCAAACCGCAAATCGCTTCGCCCTTCACCCTTCCcaca                       >  2:803402‑M2/23‑143 (MQ=255)                                                                                                                          gccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAACTCCGccagcac                   >  2:1361‑M2/19‑140 (MQ=255)                                                                                                                            cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                 >  1:250871‑M2/17‑143 (MQ=255)                                                                                                                            cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCACCGTTCACAATTCCGACAGAAACt                 >  1:751792‑M2/17‑143 (MQ=255)                                                                                                                                gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg             <  2:250871‑M2/131‑1 (MQ=255)                                                                                                                                       tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt      >  1:893309‑M2/6‑143 (MQ=255)                                                                                                                                           cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  >  1:26200‑M2/2‑143 (MQ=255)                                                                                                                                           cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCTCCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACGATGCCGCCAGGGACCGCGCCCACCatcctt  >  1:213472‑M2/2‑141 (MQ=255)                                                                                                                                            |                                                                                                                                              AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815802‑3816082

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