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	44 [0]	[0] 43	[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)	39 (0.830)	36/262	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

ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTC  >  NC_000913/3815719‑3815998                                                                                                                                             |                                                                                                                                             acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTGCcgt                                                                                                                                           >  2:280442‑M1/1‑140 (MQ=255) acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGCTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCACCTTAAACATCAGTTTCCACGCGCCGCCTACTGCagc                                                                                                                                           >  2:122539‑M1/1‑140 (MQ=255)  ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                                                          <  2:684938‑M1/143‑5 (MQ=255)                         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                   <  1:829654‑M1/143‑28 (MQ=255)                         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                   >  1:741915‑M1/1‑116 (MQ=255)                          tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGAGACAGATCAGTAGCCTTAAAAATAAGTTTGCCAGGGCAGCCTTCTGCagtcccctgcacttcaatgatgcgaccg                                                                                                                  >  1:88677‑M1/1‑115 (MQ=255)                                ggtgctgaTATGCTTTCGTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                            <  2:215817‑M1/140‑35 (MQ=255)                                 gcgcTGATATGGTTTCATGCCTTAGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgatgcgccagccttcgg                                                                                                           >  2:924230‑M1/1‑108 (MQ=255)                                   gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacatacatgctgcgcccgccctcggtg                                                                                                         >  2:267292‑M1/1‑106 (MQ=255)                                              atactgcctCCTCTCCTCATCTTCCTTTTCGACACAC‑AAAAAAAACGCGACTCATCATTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                              <  2:799388‑M1/139‑49 (MQ=255)                                               tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                             <  1:99171‑M1/143‑50 (MQ=255)                                                  tGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttccgtcctccctacgttcacg                                                                                          >  2:923525‑M1/1‑91 (MQ=255)                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactccaatgatgcgcccgtcttcggtcatcactccgttcatgtc                                                                                        >  1:684869‑M1/1‑89 (MQ=255)                                                     cTTCTCTCCTCATCTGACTTCGCTACAGAC‑GAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                                                                       <  2:872412‑M1/143‑56 (MQ=255)                                                        cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                                                    >  1:840227‑M1/1‑85 (MQ=255)                                                            cctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGTggtgtgccagacgccaatgatgcgcccgtcttcggtaatcactacgttaatgtcgggctatg                                                                                >  1:428436‑M1/1‑81 (MQ=255)                                                             ctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatctcggtctctgc                                                                               >  2:203927‑M1/1‑80 (MQ=255)                                                                 tCTTACTTTTAAACGGAA‑AGAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCGGTTTGCCAGCGCCGCCTTCTGCcgtcccccgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                           <  1:326653‑M1/143‑68 (MQ=255)                                                                  ctcccttttCTACAGAC‑AAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcag                                                                          <  1:293500‑M1/139‑69 (MQ=255)                                                                    tACTGTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                         <  1:571508‑M1/143‑70 (MQ=255)                                                                     aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctggacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggccgagt                                                                       >  2:108322‑M1/1‑72 (MQ=255)                                                                          tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca                                                                  <  2:840227‑M1/143‑77 (MQ=255)                                                                               aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaacgta                                                             >  1:20251‑M1/1‑62 (MQ=255)                                                                               aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCCTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggccatcactacgttcgtggcggccacggccgcagggtcttcaccgga                                                             >  2:144810‑M1/1‑62 (MQ=255)                                                                                         aaaGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCTTTCTGTcgatccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                                                    <  1:834130‑M1/143‑91 (MQ=255)                                                                                                     tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                        >  1:69304‑M1/1‑41 (MQ=255)                                                                                                        gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgcc                                     >  1:75232‑M1/1‑38 (MQ=255)                                                                                                              ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                               >  1:917196‑M1/1‑32 (MQ=255)                                                                                                              tgaaaaaccAGGTTGCCGTCGCCGCCTTCTGCcgtcccctgcacttccatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                               <  2:774733‑M1/141‑112 (MQ=255)                                                                                                              cgaaaaatcCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggccgagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                               <  2:592785‑M1/141‑112 (MQ=255)                                                                                                               tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctgagccgttccca                              >  2:140278‑M1/1‑31 (MQ=255)                                                                                                                  aaaTCCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgctgcagagttttcaacgtattccagatcgcaaaccgcttcgccgttcacaatt                           <  1:535296‑M1/143‑116 (MQ=255)                                                                                                                   aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                          <  1:828453‑M1/143‑117 (MQ=255)                                                                                                                      cAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccga                       <  1:564674‑M1/143‑120 (MQ=255)                                                                                                                       aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatagcaaaccgcttcgccgttctcaattacgac                      >  1:47039‑M1/1‑23 (MQ=255)                                                                                                                       aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgcggtcatcactacgttcaggtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                      <  1:374499‑M1/143‑121 (MQ=255)                                                                                                                        gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatcccagatctcaaaccgcttcgccgttcacaattccgtca                     >  1:368057‑M1/1‑22 (MQ=255)                                                                                                                         tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacag                    >  1:736888‑M1/1‑21 (MQ=255)                                                                                                                         tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgcattccagatcgcaaaccgcctcgccgttcccaactccgaccg                    >  2:234604‑M1/1‑21 (MQ=255)                                                                                                                                gcgcCGACTTCTGCcgtcccatgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgctaaccgcttcgccgttcacaattccgacagaaactgc             <  2:238404‑M1/143‑130 (MQ=255)                                                                                                                                    cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcaccattccgacagaacctgcggct         >  2:587790‑M1/1‑10 (MQ=255)                                                                                                                                      ccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgctcatgtcggtctctgcggcagagtcttcaacgcattccagatcgcaaaccgcttcgccgttgccaattccgccagaacatgcggctgc       >  1:302752‑M1/1‑8 (MQ=255)                                                                                                                                        ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca     >  1:453660‑M1/1‑6 (MQ=255)                                                                                                                                           tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatcc  <  2:741915‑M1/143‑141 (MQ=255)                                                                                                                                             |                                                                                                                                             ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTC  >  NC_000913/3815719‑3815998

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