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	49 [0]	[0] 49	[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)	46 (0.940)	38/264	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

AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815803‑3816080                                                                                                                                           |                                                                                                                                            ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcctacttttcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttccgccagcc                                                                                                                                         >  2:431617‑M2/139‑139 (MQ=255)   gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctccacctcttacacatcaacagacaaaaaaaaagcgactcatcagtcgccttaaaaatcagtttgccagcgcccccttctgcCGTcccc                                                                                                                                       >  1:632265‑M2/137‑143 (MQ=255)         taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcggccttcggcCGTCCCCTGCCCt                                                                                                                                 >  1:594000‑M2/131‑143 (MQ=255)         taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtagccttaaaaatcagtttgccagagcagccttctgcCGTCCCATGCACt                                                                                                                                 >  1:414129‑M2/131‑143 (MQ=255)          aagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                <  1:763498‑M2/14‑1 (MQ=255)          aagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                >  2:390152‑M2/130‑143 (MQ=255)               caaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                           <  1:74907‑M2/19‑1 (MQ=255)                aaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttacaaatcagtttgccagcgccgccttctgcCGCTCCACGCACATCAatga                                                                                                                          >  2:95443‑M2/124‑143 (MQ=255)                 aattcaataaactggcgctgatatggtttcatgccttcgctcctcaccttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGGCCCCTGCACTTCAatgat                                                                                                                         <  2:118864‑M2/21‑1 (MQ=255)                    tcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccctaaaaatcagtttgccagcgccgccttctgcCGGCCCCTGCACTTCAGTGATgcg                                                                                                                      >  1:309175‑M2/120‑143 (MQ=255)                    tcaataaactgacgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgacgtacaaatcaagttgccagagcacccttccgcCCTCCCCTACCCTTCAATGAtgca                                                                                                                      >  2:9213‑M2/120‑142 (MQ=255)                        taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagattgccagcgccgccttctgcCGTCCCCCGCACTTCAACGACGCgcccc                                                                                                                  >  2:146053‑M2/116‑142 (MQ=255)                          aactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                                >  2:208902‑M2/114‑143 (MQ=255)                             tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                             >  2:895059‑M2/111‑143 (MQ=255)                              ggtcctcatatcgtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                            <  1:793022‑M2/34‑1 (MQ=255)                              ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCACCCGTCTTCg                                                                                                            <  1:33669‑M2/34‑1 (MQ=255)                                cgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaacaaggcgactcatcagtcgccttaaaaatcagcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                          <  2:122189‑M2/36‑1 (MQ=255)                                       atggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaagacgactcatcagtcgccttaaaaatcagtttgcaagcgccgccttctgcCCTCCCCTGCACTTCAATGATGCTCCCGTCTTCGGTCATAACt                                                                                                   >  2:108869‑M2/101‑143 (MQ=255)                                         ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgcattatgtGGTGCTCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTAc                                                                                                 >  1:215140‑M2/99‑143 (MQ=255)                                           gttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaaccagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                               <  2:145317‑M2/47‑1 (MQ=255)                                             tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                             >  1:106794‑M2/95‑143 (MQ=255)                                                           ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtagcattaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTGTGc                                                                               >  1:264629‑M2/81‑143 (MQ=255)                                                            tcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                              <  1:876607‑M2/64‑1 (MQ=255)                                                                 ttgcttttcttcagacaaaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                         <  2:755420‑M2/69‑1 (MQ=255)                                                                 ttactttgctacagacaaaaaaaaggcgactgatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                         <  2:208722‑M2/69‑1 (MQ=255)                                                                  tacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcAGTCCACTGCACTGCAATGCTGCGCCCGTCTTCGATCACCACTCCGTCCATGTCGGtgtccgccgcccgg                                                                        >  2:172398‑M2/74‑139 (MQ=255)                                                                     ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                     >  1:749025‑M2/71‑143 (MQ=255)                                                                      ttttgacagacaaaaaaaaggcgcctcaccagtcgccttaaaagccagcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                    <  2:761044‑M2/74‑1 (MQ=255)                                                                          tacagacaaaaaaaaggcgactcatcggtcgccttaaaaaccagtgtgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                <  1:87818‑M2/78‑1 (MQ=255)                                                                             agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCGTCACTCCGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTa                                                             >  1:83557‑M2/63‑143 (MQ=255)                                                                               acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGCTCATCACTCCGTTCATGTCGGTCTCTCCGGCAGCGCCTTCAACGTCtt                                                           >  2:804120‑M2/61‑143 (MQ=255)                                                                                 aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGCCAGAGCCTTCAACGTATTcc                                                         >  1:796417‑M2/59‑143 (MQ=255)                                                                                  aaaaaaaggcgactcatcggtcgccttaaaaatcagtttgccagcgccgcctgctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                        <  2:392549‑M2/86‑1 (MQ=255)                                                                                   aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCcac                                                       >  2:541096‑M2/57‑142 (MQ=255)                                                                                   aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGCTCACGTCGCTCTCTGCGGCAGAGTCTTCAACGTCTTCcac                                                       >  2:173423‑M2/57‑142 (MQ=255)                                                                                      aaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                    <  1:501775‑M2/90‑1 (MQ=255)                                                                                              ctcatcagtcgccttaaaaatcagtttgccggcgccgccttctgcCGTCCCCTGCGCTTCAATGATGCGCCCGTCTGCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCCAACCg                                            >  1:524078‑M2/46‑143 (MQ=255)                                                                                                     gtcgccttaaaaatgggtttgccagcgccgccttctgcCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGcc                                     <  2:220206‑M2/105‑1 (MQ=255)                                                                                                      tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGTAAgacgcttgcccg                                    >  1:845086‑M2/38‑138 (MQ=255)                                                                                                      tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                    >  1:620096‑M2/38‑143 (MQ=255)                                                                                                         ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                 >  1:160617‑M2/35‑143 (MQ=255)                                                                                                           ttaaaaatcagtttgccagcgccgccttctgcCGTCTCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                               >  2:314701‑M2/33‑143 (MQ=255)                                                                                                                   cgttttggcagcgccgccttctgcCGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGTCGTTCACAATTTcgt                       <  2:936672‑M2/119‑2 (MQ=255)                                                                                                                    agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGCTCACTGCGACAGAGTCTTCAACGTATTGCAGATCGCAAACCGCGTCGCCGTTCTCAATCCAGAc                      >  2:827632‑M2/24‑143 (MQ=255)                                                                                                                              cgccgccttctgcCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg            <  1:390152‑M2/130‑1 (MQ=255)                                                                                                                               gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCCTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg           <  2:571981‑M2/131‑1 (MQ=255)                                                                                                                                   cgttctggCGTCTCCTGCACTTCAATGATGCGCCCGTCCGTGTGTAGCAATACGATCATGTCGGTCTTTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc       <  2:197967‑M2/135‑1 (MQ=255)                                                                                                                                     ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGGATTCACAAACGCAGACCGCTTCGCCGCTCACAATTCCGACAGAAACTGCGGCTaccc     >  1:439123‑M2/7‑142 (MQ=255)                                                                                                                                        tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc  >  1:35580‑M2/4‑143 (MQ=255)                                                                                                                                           |                                                                                                                                            AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815803‑3816080

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