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	80 [0]	[0] 80	[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)	72 (0.870)	53/178	0.5	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

ATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCA  >  NC_000913/3815759‑3815956                                                                                                     |                                                                                                   atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                    >  1:577362‑M1/1‑100 (MQ=255)  tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGAA‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                   <  1:1518426‑M1/100‑2 (MQ=255)   ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                  >  2:1653888‑M1/1‑98 (MQ=255)      ttCATGCCTTCGCTCCTCATTTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                               <  2:312151‑M1/100‑6 (MQ=255)           gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgt                                                                                          <  2:266387‑M1/100‑11 (MQ=255)            ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                         <  1:1815916‑M1/100‑12 (MQ=255)            ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgaa                                                                                         >  1:736906‑M1/1‑89 (MQ=255)             cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                                        >  2:2121820‑M1/1‑88 (MQ=255)             cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                                        <  2:357949‑M1/100‑13 (MQ=255)             cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                                        >  2:1563362‑M1/1‑88 (MQ=255)               tCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                                      >  1:2180924‑M1/1‑86 (MQ=255)                  ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                   <  2:279248‑M1/100‑18 (MQ=255)                  ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                   <  2:1738146‑M1/100‑18 (MQ=255)                    cgtcatCTTACTTGTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                 >  2:1504901‑M1/3‑81 (MQ=255)                    cctcATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                                                                    >  1:207285‑M1/1‑82 (MQ=255)                    cctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                 >  1:259164‑M1/1‑81 (MQ=255)                           ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                          <  1:142739‑M1/100‑27 (MQ=255)                           ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                          <  1:1543072‑M1/100‑27 (MQ=255)                             aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                                                        <  2:1885261‑M1/100‑29 (MQ=255)                              cTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                                                       >  2:1023212‑M1/1‑71 (MQ=255)                                 ttCTACAGGC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                    <  2:812270‑M1/100‑33 (MQ=255)                                  tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                   >  1:1770664‑M1/1‑67 (MQ=255)                                   cTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                  >  1:754369‑M1/1‑66 (MQ=255)                                    tACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                 >  2:793570‑M1/1‑65 (MQ=255)                                      cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                               >  2:1592768‑M1/1‑63 (MQ=255)                                       aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                                                              <  1:1134868‑M1/100‑39 (MQ=255)                                       aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                                                              <  2:1548309‑M1/100‑39 (MQ=255)                                        gAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgacgcgcccgtcttcggtcat                                                             >  1:2216512‑M1/1‑61 (MQ=255)                                          c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  1:1259737‑M1/1‑59 (MQ=255)                                          c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  2:239974‑M1/1‑59 (MQ=255)                                            aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                                                          <  1:1864760‑M1/100‑43 (MQ=255)                                            aaaaaaaaGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                                                          <  1:560022‑M1/100‑43 (MQ=255)                                                aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGTGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                      >  1:342701‑M1/1‑54 (MQ=255)                                                 aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                     <  1:852786‑M1/100‑48 (MQ=255)                                                 aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgn                                                     <  1:1051555‑M1/100‑48 (MQ=255)                                                   aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                   <  1:1878997‑M1/100‑50 (MQ=255)                                                    ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca                                                  <  2:440495‑M1/100‑51 (MQ=255)                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCGAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                >  1:2170503‑M1/1‑48 (MQ=255)                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                <  2:878277‑M1/100‑53 (MQ=255)                                                          tcatcaGTCGCCTTAAAAATCAGTTTGCCAGTGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                            <  2:342702‑M1/100‑57 (MQ=255)                                                           catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt                                           >  1:1891408‑M1/1‑43 (MQ=255)                                                            atcaGTCGCCTTAAAAATCAGTTTGCGAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                          <  2:2170515‑M1/100‑59 (MQ=255)                                                            atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttctgtcatcactacgttcatgtcggtc                                          <  2:1732816‑M1/100‑59 (MQ=255)                                                            atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                          >  2:215737‑M1/1‑42 (MQ=255)                                                             tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                         >  1:1821235‑M1/1‑41 (MQ=255)                                                              caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                        >  2:1833732‑M1/1‑40 (MQ=255)                                                                 tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                     <  2:577363‑M1/100‑64 (MQ=255)                                                                 tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                     <  1:244268‑M1/100‑64 (MQ=255)                                                                 tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                     <  2:1821245‑M1/100‑64 (MQ=255)                                                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgg                                   >  2:468412‑M1/1‑35 (MQ=255)                                                                     cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                 >  2:58296‑M1/1‑33 (MQ=255)                                                                     cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                 <  1:272540‑M1/100‑68 (MQ=255)                                                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcag                                >  1:1156096‑M1/1‑32 (MQ=255)                                                                       tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                               >  2:919755‑M1/1‑31 (MQ=255)                                                                        aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                              >  1:343250‑M1/1‑30 (MQ=255)                                                                         aaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                             >  1:1444666‑M1/1‑29 (MQ=255)                                                                          aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                            >  2:18290‑M1/1‑28 (MQ=255)                                                                          aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                            <  2:499627‑M1/100‑73 (MQ=255)                                                                          aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcctcggtcatcactacgttcatgtccgtctctgcggcacagtc                            >  1:1886562‑M1/1‑28 (MQ=255)                                                                           tccgatcttgccagcgCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                           <  2:1575473‑M1/93‑74 (MQ=255)                                                                             tCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgctgcagagtcttc                         <  1:1991779‑M1/100‑76 (MQ=255)                                                                              cAGTTTGCCTGCGCCGCCTTCTGCcgtcccctgcacttcaagggtgcgcccgtcttcggtcatcactacgttcaggtcgctctctgcggcagagtcttc                         <  1:415368‑M1/99‑76 (MQ=255)                                                                               aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                       >  2:919996‑M1/1‑23 (MQ=255)                                                                                gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                      <  2:1925121‑M1/100‑79 (MQ=255)                                                                                 tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                     <  2:1891418‑M1/100‑80 (MQ=255)                                                                                 tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgacgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                     <  2:2216526‑M1/100‑80 (MQ=255)                                                                                   tGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctagatcgga                   >  1:1575465‑M1/1‑19 (MQ=255)                                                                                       aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc               >  2:2043368‑M1/1‑15 (MQ=255)                                                                                       aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc               >  1:1786733‑M1/1‑15 (MQ=255)                                                                                       aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc               <  1:1833722‑M1/100‑86 (MQ=255)                                                                                          gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga            >  2:1338216‑M1/1‑12 (MQ=255)                                                                                           ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat           <  1:1023211‑M1/100‑90 (MQ=255)                                                                                           ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat           >  1:70378‑M1/1‑11 (MQ=255)                                                                                            cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc          <  1:1504894‑M1/100‑91 (MQ=255)                                                                                                 tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaa     >  2:1265889‑M1/1‑5 (MQ=255)                                                                                                 tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaa     >  2:1517124‑M1/1‑5 (MQ=255)                                                                                                  cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac    <  2:1331297‑M1/100‑97 (MQ=255)                                                                                                   tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaacc   <  2:1259739‑M1/100‑98 (MQ=255)                                                                                                    gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg  >  1:661789‑M1/1‑2 (MQ=255)                                                                                                    gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg  >  2:1340513‑M1/1‑2 (MQ=255)                                                                                                     |                                                                                                   ATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCA  >  NC_000913/3815759‑3815956

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