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	57 [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)	50 (0.710)	45/264	0.5	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

AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815801‑3816082                                                                                                                                             |                                                                                                                                              acctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcctattccttttatacagacaaaaaaaaggcgactcatgagtcgccttagaaatcagtttggccgcgccgccttctgcCGt                                                                                                                                             >  1:398467‑M2/141‑143 (MQ=255)  cctgcttgctaagcgcaaattcaataaaatggcgatggtatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                            <  1:690180‑M2/4‑1 (MQ=255)   ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGccc                                                                                                                                           >  1:641927‑M2/139‑143 (MQ=255)   ctgcttgctaagcgcaaattcaataaactggcgctgacatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                           >  2:853696‑M2/139‑143 (MQ=255)    tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgacgccttctgcCGTccc                                                                                                                                          >  1:658855‑M2/138‑143 (MQ=255)      cttgctacgcgcaaattccataaactggcgctgatatggttccatgccttcgctcctcatcttacttttctaccgacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                                        <  1:1094683‑M2/8‑1 (MQ=255)      cttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgcccctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                                        <  1:392863‑M2/8‑1 (MQ=255)         ggtgagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                     <  1:1186138‑M2/11‑1 (MQ=255)          ctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                                                                    >  1:150439‑M2/132‑143 (MQ=255)           taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgaatcatctgactcggtaaaaagcagtttgccagcgccgccttcggcCGTCACCTGCACt                                                                                                                                   >  1:624918‑M2/131‑143 (MQ=255)             agcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCGCTTc                                                                                                                                 >  1:1052031‑M2/129‑143 (MQ=255)               gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCTCTTCaa                                                                                                                               >  1:1007066‑M2/127‑143 (MQ=255)                gtaaattccatacactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                              <  2:412393‑M2/18‑1 (MQ=255)                        aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttcaaaatcactttgccagcgccgccttctgcCGTCCCCTGCACTTCAATTAtgcgct                                                                                                                      >  1:441002‑M2/118‑142 (MQ=255)                        aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGcc                                                                                                                      >  2:1059447‑M2/118‑143 (MQ=255)                                tgcgcttatatggtttcatgcctgcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  2:1210521‑M2/34‑1 (MQ=255)                                   gctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCATTGATCCGCCCGTCTTCGGtc                                                                                                           >  1:697742‑M2/107‑143 (MQ=255)                                           ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtccccttacaaatcagtttgccagcgccgccctctgcCGTCCCCTGCACTTGAATGATGCGCCCGTCTTCGCTCCTCCCTAc                                                                                                   >  1:222026‑M2/99‑143 (MQ=255)                                            gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                  >  2:97724‑M2/98‑143 (MQ=255)                                             tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCACTAcgc                                                                                                 >  2:910848‑M2/97‑142 (MQ=255)                                                 atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                             <  1:710584‑M2/51‑1 (MQ=255)                                                 atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGAAGCGGCCGTCATCGGTCATCACTACGctcct                                                                                             >  2:333427‑M2/93‑141 (MQ=255)                                                   tgccttcgctcctcatcttacttttttacagacaaaaaaaaggcgactcatcagtcgccttaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                           <  2:197030‑M2/53‑1 (MQ=255)                                                   gccttcgctcttcatcttacttttctacagataacaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                           <  2:1095720‑M2/53‑1 (MQ=255)                                                    ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTGCTCACTACGTTCATGTc                                                                                          >  1:1271238‑M2/90‑143 (MQ=255)                                                      ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                                                        >  1:909362‑M2/88‑143 (MQ=255)                                                               catctgacttttctacagacaaaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                               <  1:67139‑M2/65‑1 (MQ=255)                                                                   ttacttttgtagaggcaaaaagaaggcgactcatcagtcggcttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                           <  2:1160450‑M2/69‑1 (MQ=255)                                                                       ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGGCGGTCGCTCCGGCCGCGTCt                                                                       >  1:545718‑M2/71‑143 (MQ=255)                                                                         ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                     <  1:904453‑M2/75‑1 (MQ=255)                                                                            tacagacaaaaaacgggcgatgcatcagtcgccttaaaaatcagtttgccagcgccgccctctgcCGTTCCCTCCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                  <  2:1104637‑M2/78‑1 (MQ=255)                                                                                gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGCGCTCCGCGGCAGAGTCCTCAAggtct                                                              >  2:1199980‑M2/62‑141 (MQ=255)                                                                                 acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCAGAGTCCTCAACGTCtt                                                             >  2:449777‑M2/61‑143 (MQ=255)                                                                                    aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                          <  1:879909‑M2/86‑1 (MQ=255)                                                                                       caaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                       <  2:1078520‑M2/89‑1 (MQ=255)                                                                                           ggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCa                                                   >  2:380466‑M2/51‑143 (MQ=255)                                                                                            gcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCATCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCaa                                                  >  2:610144‑M2/50‑143 (MQ=255)                                                                                             cgactcatcctccgccttataaatcagtttgccagcgccgccttctgcCGTCTCCTGCACTTCAATGATGCGCCCCTCTTCGGTCATCACTACGTTCATGTCGGTCTCTTCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                                 <  2:922733‑M2/95‑1 (MQ=255)                                                                                                           cattaaaattcactttgccagcgccgccttctgcCGTCCCCTGCACCGCAACGATGCGCCCGTCTTCGTTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTGCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                   <  1:8372‑M2/109‑1 (MQ=255)                                                                                                               aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACaa                               <  1:384421‑M2/113‑1 (MQ=255)                                                                                                                 aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAATCTTAAACGTATTCCAAATCACAAAACGGTTCGCATTTGACAAtt                             >  1:129529‑M2/29‑143 (MQ=255)                                                                                                                 aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtt                             <  2:1273562‑M2/115‑1 (MQ=255)                                                                                                                     ctgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                         <  1:146738‑M2/119‑1 (MQ=255)                                                                                                                       gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTAATGCCGGTCTCTACGGCAGAGTCTTCAACGCCTTCCAGCGCGCCAACCGCTTCGCCCTTCACCATTCCGacc                       >  2:22589‑M2/23‑142 (MQ=255)                                                                                                                        tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGACCGCACACCGCTCCGCCCCTCACACGTCCCACAg                      >  2:838966‑M2/22‑143 (MQ=255)                                                                                                                              cgcgccgccttctgcCGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAAGGATTCCAGATCGCAGACCGCTTCGCCGTTTACAATTCCGACAGGAACTg                >  1:292019‑M2/16‑143 (MQ=255)                                                                                                                               gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc               <  1:1309531‑M2/129‑1 (MQ=255)                                                                                                                               gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGACATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTGACAATTCAGACAGAAACTGc               >  1:1072874‑M2/15‑143 (MQ=255)                                                                                                                                cgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg              <  1:427205‑M2/130‑1 (MQ=255)                                                                                                                                  ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGc            >  1:3017‑M2/12‑143 (MQ=255)                                                                                                                                       ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa       >  1:92552‑M2/7‑143 (MQ=255)                                                                                                                                           gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGACCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTCCCATccc   >  1:576869‑M2/3‑143 (MQ=255)                                                                                                                                            cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTCCCATCCCt  >  2:174854‑M2/2‑143 (MQ=255)                                                                                                                                            cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  >  2:998015‑M2/2‑143 (MQ=255)                                                                                                                                             |                                                                                                                                              AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815801‑3816082

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