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

AAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGG  >  NC_000913_3_bme_pgi/3815709‑3816005                                                                                                                                                       |                                                                                                                                                    aaaCTTTAACACCTGCTTGCTAAGCGCAAATTAAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTATACACAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                                                                     <  1:304412‑M1/150‑1 (MQ=255)      ttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGGTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                                                                <  1:638163‑M1/150‑6 (MQ=255)           acCTGCTTTCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                                                                           <  2:22043‑M1/150‑11 (MQ=255)                     tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                                                                                 <  2:672776‑M1/150‑21 (MQ=255)                       aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑TAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTGTGCcgtcccgggcagttctatgatg                                                                                                                               >  1:728793‑M1/1‑128 (MQ=255)                         cgcTAATTCAATAACCTGGCGCTGATATGGTTCCCTGCCTTCGCTCCTCATCTTACTTTACAACAGAC‑AAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcg                                                                                                                             <  2:805575‑M1/150‑25 (MQ=255)                           cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTACTGCcgtcccctccaattcaatgatgcgcc                                                                                                                           >  1:71015‑M1/1‑124 (MQ=255)                            aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                          >  2:98911‑M1/1‑123 (MQ=255)                            aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                          >  2:676478‑M1/1‑123 (MQ=255)                                 cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgacccctgcacttcaatgatgcgcccgtctt                                                                                                                     >  1:604020‑M1/1‑118 (MQ=255)                                 cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCAGCcgtcccctgcccttcaatgaagcgcccgtctt                                                                                                                     >  2:904914‑M1/1‑118 (MQ=255)                                     aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                                 <  2:276672‑M1/150‑37 (MQ=255)                                      aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggccccggcacttcaatgatgcgcccgtcttcggtc                                                                                                                >  2:587714‑M1/1‑113 (MQ=255)                                        cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGAcgtcccctgcacttcaatgatgcgcccgtcttcggtgat                                                                                                              >  1:118874‑M1/1‑111 (MQ=255)                                                       tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                                                                               <  1:723263‑M1/150‑55 (MQ=255)                                                           aTGCCTTCGCTCCTCATCTTACTTTTCTACACAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                                                           <  1:435849‑M1/150‑59 (MQ=255)                                                                       ctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                                                                               >  1:859578‑M1/1‑80 (MQ=255)                                                                          aTCTAACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTAAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcaaagatgagcccgtcttcggtcatcagtacgttcagttcggtctctgcggcagattc                                                                             >  2:742955‑M1/1‑78 (MQ=255)                                                                             ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacggtgatgtcggtctctgcggcagagtcttca                                                                         >  1:766195‑M1/1‑74 (MQ=255)                                                                                        cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                                                              <  1:576447‑M1/150‑88 (MQ=255)                                                                                            c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaattatgcgcccgtcttcggtcatcactacgttcatgtcggtctcggcggcagcgtcttcaccgtattcacgatcg                                                          >  2:816789‑M1/1‑59 (MQ=255)                                                                                                aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaa                                                       >  2:943521‑M1/1‑56 (MQ=255)                                                                                                   caagGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaacc                                                    >  2:493428‑M1/2‑53 (MQ=255)                                                                                                     aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgc                                                  >  1:583123‑M1/1‑51 (MQ=255)                                                                                                           cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttcaacgtattacagatcgcaaaccgcttcgcc                                            >  1:229194‑M1/1‑45 (MQ=255)                                                                                                                     gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttcctgtcggtctctgcggcggagtcttcaacggattccagatcgcaaaccgcttcgccgttcccaatt                                  >  1:703762‑M1/1‑35 (MQ=255)                                                                                                                        ataAAAATCAGTTTGCCAGCGCCGCCTTCTGAcgtcccgtgcacttcattgatgctcacgtcttcggtcatcacaaggttaatgtcggtgtctgcagcagagtcttcaacgtattccagatcgcaaacagcttcggcgttcacaattccg                               >  1:542560‑M1/2‑32 (MQ=255)                                                                                                                           aaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccggcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagctggcaaaccgcttcgacgttcacaattccgcca                            >  1:564887‑M1/1‑29 (MQ=255)                                                                                                                            aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgtattccagatcacaaaccgcttcgccgtgcacaattccgacag                           >  1:381750‑M1/1‑28 (MQ=255)                                                                                                                            aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacag                           >  2:98773‑M1/1‑28 (MQ=255)                                                                                                                             aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttacacgtattccagatcgcaaaccgcttcgccgggcccaattccgacaga                          >  1:414124‑M1/1‑27 (MQ=255)                                                                                                                              aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcatcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagctcgcaaaccgcttcgccgttccccattccgacagaa                         >  2:812414‑M1/1‑26 (MQ=255)                                                                                                                                cAGTTTGCCAGCGCCTCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggactctgcggcagagtcttcaacgtattccagatcgcaaaccacttcgccgttcacaattccgacagaaac                       <  2:71015‑M1/150‑127 (MQ=255)                                                                                                                                     tGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcatggcggtctctgcggcagagtcttccacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg                  <  2:798942‑M1/150‑132 (MQ=255)                                                                                                                                        cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggcta               >  1:323610‑M1/1‑16 (MQ=255)                                                                                                                                          gcgcCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaggttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacc             <  1:491259‑M1/150‑137 (MQ=255)                                                                                                                                            gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccat           >  1:792936‑M1/1‑12 (MQ=255)                                                                                                                                               gccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccc        <  2:750045‑M1/150‑142 (MQ=255)                                                                                                                                                     tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcat  <  1:98911‑M1/150‑148 (MQ=255)                                                                                                                                                     tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcat  <  1:676478‑M1/150‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                    AAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGG  >  NC_000913_3_bme_pgi/3815709‑3816005

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