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	53 [0]	[0] 52	[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.770)	42/264	0.4	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                                                                                                                                             |                                                                                                                                              acctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                                                             >  2:1081333‑M2/141‑143 (MQ=255)    tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          >  1:576746‑M2/138‑143 (MQ=255)    tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          >  1:356168‑M2/138‑143 (MQ=255)         gctaagcgcaaatccaatacactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                     <  1:79727‑M2/11‑1 (MQ=255)          ctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                                                                    >  1:14049‑M2/132‑143 (MQ=255)           taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgacgcatcagtcgccttagaaatcagttagccagcgccgccttctcaCGTCCCCTGCACt                                                                                                                                   >  2:72787‑M2/131‑143 (MQ=255)            aagcgcaaattcaataaacgggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                  <  1:932898‑M2/14‑1 (MQ=255)               cgcaaattcaataaactggcgctgatatggttccatgcatccgctcctactttgacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                               <  1:789003‑M2/17‑1 (MQ=255)                  aaattcaataaactggcgctgatatggtttcatgccctcgctcatcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                                                            <  1:507472‑M2/20‑1 (MQ=255)                         ataagctgtcgctgttgtggtttcatgccttcgctctttatcttacttttctacagacaaaaaacaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                     <  2:560366‑M2/27‑1 (MQ=255)                         ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgacccatcagtcaccttaaaaatcagtttgcaagggccgccttctgcCATCCCCAACACTTCAATGACGCGccc                                                                                                                     >  2:480127‑M2/117‑143 (MQ=255)                                 gcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTCCCCCCGTCTTCgg                                                                                                             >  2:220206‑M2/109‑143 (MQ=255)                                   gctgatattgtttcatgccttcgctcctcatcttactcttctacagacaaaaaaaaggcgactcttcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                           <  1:107058‑M2/37‑1 (MQ=255)                                   gctgatatggttccctggctgcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                           <  2:452864‑M2/37‑1 (MQ=255)                                         atggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                     >  1:1005522‑M2/101‑143 (MQ=255)                                          gggtttcatgccctcgctccccctctcaaattctttcggtctaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                    <  1:896554‑M2/44‑1 (MQ=255)                                            gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                  <  1:452316‑M2/46‑1 (MQ=255)                                              ttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                >  1:515165‑M2/96‑143 (MQ=255)                                               tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                               <  1:882851‑M2/49‑1 (MQ=255)                                                catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACCACGTTCa                                                                                              >  1:271093‑M2/94‑143 (MQ=255)                                                 atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                             <  1:89840‑M2/51‑1 (MQ=255)                                                  tgccttcgctcctcctcttactcttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgcctcctgcCGTCCCCTGCACTTCAATGATGCCCCCGTCTTCGGTCATCACTACGTTCATg                                                                                            <  1:442986‑M2/52‑1 (MQ=255)                                                     cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                         >  2:450799‑M2/89‑143 (MQ=255)                                                           tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgagagcgccgccttctgcCGTCCCCTGCACTTCAATACTGCCCCCGTCGTCGGTCATCACTACGTTCACGTCGCTccct                                                                                   >  2:1124800‑M2/83‑143 (MQ=255)                                                            cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                  <  1:272833‑M2/62‑1 (MQ=255)                                                               catcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCgg                                                                               >  1:469653‑M2/79‑143 (MQ=255)                                                               catcgtacttttctacagtctataaaaatgcgactcatcagtcgccttaaaaatcattttgtcagcgccgccttctgcCGTCGCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                               <  1:190503‑M2/65‑1 (MQ=255)                                                                   ttacttttctacagacaaaacaaagccgagtaatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCaga                                                                           >  2:520529‑M2/75‑143 (MQ=255)                                                                   ttacttttctacaaaccaaagaactgcgtctcatcagtcgccttaaaaatcagtttgccagcgccgccttccgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                           <  2:750054‑M2/69‑1 (MQ=255)                                                                       ttttctacagacaaaaaaaaggcgactcatcagtcaccttaaaaatcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTGCGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                       >  1:944180‑M2/71‑143 (MQ=255)                                                                            tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCACTCCGGCAgagcctcccac                                                                  >  2:1037266‑M2/66‑138 (MQ=255)                                                                             acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                                 >  2:88258‑M2/65‑143 (MQ=255)                                                                                     aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg                                                         >  1:819857‑M2/57‑143 (MQ=255)                                                                                      aaaaagtcgccccatcagttgctttaaaaattagtttgccagctccgccttctgtCGCCTCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCTGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                        <  2:288377‑M2/88‑1 (MQ=255)                                                                                         aaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTCGGGCGGAGTCTTCACCGCACTCcagacgg                                                     >  2:425848‑M2/53‑140 (MQ=255)                                                                                            gcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCAGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaa                                                  >  1:188773‑M2/50‑143 (MQ=255)                                                                                                 tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCAGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTGCAACGTATTCCAGAACGCAAACCGc                                             >  1:126493‑M2/45‑143 (MQ=255)                                                                                                   atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTCCACTTAAATGATGCACCCGTCTTCGGTCCTCACTACCTTCATGTCCGTCTCTGCGGCCGACTCTTCAACGTCTTCCACATCGCACACCCCtt                                           >  2:146592‑M2/43‑143 (MQ=255)                                                                                                      agtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCAAAccgctcagc                                        >  2:1107668‑M2/40‑139 (MQ=255)                                                                                                       gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCTCCAACGTATTCCAGCTCGCAACCCGCCTCGcc                                       >  2:1104734‑M2/39‑143 (MQ=255)                                                                                                       gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGACTCTACGAAAGAGCATGCAACCTATTCCAGATCGCAAACCGCCTGGcc                                       >  2:238009‑M2/39‑143 (MQ=255)                                                                                                        tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTACTCCAGATCGCCACCCCCCTCGCCg                                      >  2:807081‑M2/38‑143 (MQ=255)                                                                                                        tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                      >  1:1047086‑M2/38‑143 (MQ=255)                                                                                                            cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCAAACCGCTTCGCCGTtcc                                  >  1:243360‑M2/34‑142 (MQ=255)                                                                                                            cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                                  <  2:576746‑M2/110‑1 (MQ=255)                                                                                                                     cagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                         <  1:734021‑M2/119‑1 (MQ=255)                                                                                                                                    gccgtctgcCGTCCCCTGCACATCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTa          <  1:440902‑M2/134‑1 (MQ=255)                                                                                                                                     ttacgcacCGACCCCCGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc         <  1:1062961‑M2/135‑1 (MQ=255)                                                                                                                                     ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCAATAAAAATACCGACAGAAACAGCGGCTAc         >  1:389566‑M2/9‑143 (MQ=255)                                                                                                                                       ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa       >  1:859800‑M2/7‑143 (MQ=255)                                                                                                                                            cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  <  1:1040612‑M2/142‑1 (MQ=255)                                                                                                                                            cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  <  1:1039588‑M2/142‑1 (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: ‑