Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_bme_pgi	3,815,859	Δ82 bp		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	3815859	3815940	82	40 [0]	[0] 39	[rph]–[rph]	[rph],[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 3815858	0 (0.000)	36 (0.640)	31/278	0.8	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
	?	NC_000913_3_bme_pgi	3815941 =	0 (0.000)					pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

CAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAT  >  NC_000913_3_bme_pgi/3815796‑3816087                                                                                                                                                  |                                                                                                                                                   ttaacacctgcttgctaagcgcaaattcaataaactggcggtgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                                <  1:638163‑M2/5‑1 (MQ=255)      acctgctttctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGc                                                                                                                                           <  2:22043‑M2/10‑1 (MQ=255)                taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                                                                 <  2:672776‑M2/20‑1 (MQ=255)                  agcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagactaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttgtgcCGTCCCGGGCAGTTCTatgatg                                                                                                                               >  1:728793‑M2/129‑150 (MQ=255)                    cgctaattcaataacctggcgctgatatggttccctgccttcgctcctcatcttactttacaacagacaaaaaaaaggcgactcctcagtcgccttaaaaatcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcg                                                                                                                             <  2:805575‑M2/24‑1 (MQ=255)                      caaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaataagtttgccagcgccgcctactgcCGTCCCCTCCAATTCAATGATGCGcc                                                                                                                           >  1:71015‑M2/125‑150 (MQ=255)                       aaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                          >  2:98911‑M2/124‑150 (MQ=255)                       aaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                          >  2:676478‑M2/124‑150 (MQ=255)                            caataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGACCCCTGCACTTCAATGATGCGCCCGTCtt                                                                                                                     >  1:604020‑M2/119‑150 (MQ=255)                            caataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttcagcCGTCCCCTGCCCTTCAATGAAGCGCCCGTCtt                                                                                                                     >  2:904914‑M2/119‑150 (MQ=255)                                aaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatctgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                 <  2:276672‑M2/36‑1 (MQ=255)                                 aactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGGCCCCGGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                                >  2:587714‑M2/114‑150 (MQ=255)                                   ctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcagccttctgaCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtgat                                                                                                              >  1:118874‑M2/112‑150 (MQ=255)                                                  tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                               <  1:723263‑M2/54‑1 (MQ=255)                                                      atgccttcgctcctcatcttacttttctacacacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                                                                           <  1:435849‑M2/58‑1 (MQ=255)                                                                  ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCagag                                                                               >  1:859578‑M2/81‑150 (MQ=255)                                                                    atctaacttttctacagacaaaaaaaaaggcgactcatcagtcgcctaaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAAAGATGAGCCCGTCTTCGGTCATCAGTACGTTCAGTTCGGTCTCTGCGGCAGATTc                                                                             >  2:742955‑M2/79‑150 (MQ=255)                                                                        ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTGATGTCGGTCTCTGCGGCAGAGTCTTCa                                                                         >  1:766195‑M2/75‑150 (MQ=255)                                                                                   cagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg                                                              <  1:576447‑M2/87‑1 (MQ=255)                                                                                       caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCGGCAGCGTCTTCACCGTATTCACGATCg                                                          >  2:816789‑M2/60‑150 (MQ=255)                                                                                          aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaa                                                       >  2:943521‑M2/57‑150 (MQ=255)                                                                                             caaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                                    >  2:493428‑M2/54‑150 (MQ=255)                                                                                               aggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGc                                                  >  1:583123‑M2/52‑150 (MQ=255)                                                                                                     ctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTACAGATCGCAAACCGCTTCGcc                                            >  1:229194‑M2/46‑150 (MQ=255)                                                                                                               gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCCTGTCGGTCTCTGCGGCGGAGTCTTCAACGGATTCCAGATCGCAAACCGCTTCGCCGTTCCCAAtt                                  >  1:703762‑M2/36‑150 (MQ=255)                                                                                                                  ataaaaatcagtttgccagcgccgccttctgaCGTCCCGTGCACTTCATTGATGCTCACGTCTTCGGTCATCACAAGGTTAATGTCGGTGTCTGCAGCAGAGTCTTCAACGTATTCCAGATCGCAAACAGCTTCGGCGTTCACAATTCCg                               >  1:542560‑M2/33‑150 (MQ=255)                                                                                                                     aaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTGGCAAACCGCTTCGACGTTCACAATTCCGCCa                            >  1:564887‑M2/30‑150 (MQ=255)                                                                                                                      aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCACAAACCGCTTCGCCGTGCACAATTCCGACAg                           >  1:381750‑M2/29‑150 (MQ=255)                                                                                                                      aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAg                           >  2:98773‑M2/29‑150 (MQ=255)                                                                                                                       aatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTACACGTATTCCAGATCGCAAACCGCTTCGCCGGGCCCAATTCCGACAGa                          >  1:414124‑M2/28‑150 (MQ=255)                                                                                                                        atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCATCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCAAACCGCTTCGCCGTTCCCCATTCCGACAGaa                         >  2:812414‑M2/27‑150 (MQ=255)                                                                                                                          cagtttgccagcgcctccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGACTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCACTTCGCCGTTCACAATTCCGACAGAAAc                       <  2:71015‑M2/126‑1 (MQ=255)                                                                                                                               tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGGCGGTCTCTGCGGCAGAGTCTTCCACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                  <  2:798942‑M2/131‑1 (MQ=255)                                                                                                                                  cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTa               >  1:323610‑M2/17‑150 (MQ=255)                                                                                                                                    gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAGGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAcc             <  1:491259‑M2/136‑1 (MQ=255)                                                                                                                                      gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt           >  1:792936‑M2/13‑150 (MQ=255)                                                                                                                                         gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATccc        <  2:750045‑M2/141‑1 (MQ=255)                                                                                                                                               tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAt  <  1:98911‑M2/147‑1 (MQ=255)                                                                                                                                               tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAt  <  1:676478‑M2/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCAT  >  NC_000913_3_bme_pgi/3815796‑3816087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑