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

CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCT  >  NC_000913/3815842‑3816040                                                                                                    |                                                                                                    tggtttcatgccttcgctcctcatcttacttttatacagaaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                     <  1:1518426‑M2/1‑1 (MQ=255)  ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                    >  2:1653888‑M2/99‑100 (MQ=255)     ttcatgccttcgctcctcattttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                 <  2:312151‑M2/5‑1 (MQ=255)          gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCtgt                                                                                            <  2:266387‑M2/10‑2 (MQ=255)           ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                           <  1:1815916‑M2/11‑1 (MQ=255)           ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaataagtttgccagcgccgccttctgcCGTCCCctgaa                                                                                           >  1:736906‑M2/90‑98 (MQ=255)            cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                          >  2:2121820‑M2/89‑100 (MQ=255)            cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                          <  2:357949‑M2/12‑1 (MQ=255)            cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                          >  2:1563362‑M2/89‑100 (MQ=255)              tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                        >  1:2180924‑M2/87‑100 (MQ=255)                 ctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                     <  2:279248‑M2/17‑1 (MQ=255)                 ctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                     <  2:1738146‑M2/17‑1 (MQ=255)                  cctcatcttacttttctacagacaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCa                                                                                      >  1:207285‑M2/83‑98 (MQ=255)                   cgtcatcttacttgtctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                   >  2:1504901‑M2/82‑100 (MQ=255)                   cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                   >  1:259164‑M2/82‑100 (MQ=255)                          ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                            <  1:1543072‑M2/26‑1 (MQ=255)                          ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                            <  1:142739‑M2/26‑1 (MQ=255)                            acttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCg                                                                          <  2:1885261‑M2/28‑1 (MQ=255)                             cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                         >  2:1023212‑M2/72‑100 (MQ=255)                                ttctacaggcaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                                                                      <  2:812270‑M2/32‑1 (MQ=255)                                 tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                     >  1:1770664‑M2/68‑100 (MQ=255)                                  ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                    >  1:754369‑M2/67‑100 (MQ=255)                                   tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                   >  2:793570‑M2/66‑100 (MQ=255)                                     cagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                 >  2:1592768‑M2/64‑100 (MQ=255)                                      agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                <  1:1134868‑M2/38‑1 (MQ=255)                                      agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                <  2:1548309‑M2/38‑1 (MQ=255)                                       gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGACGCGCCCGTCTTCGGtcat                                                               >  1:2216512‑M2/62‑100 (MQ=255)                                         caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             >  2:239974‑M2/60‑100 (MQ=255)                                         caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             >  1:1259737‑M2/60‑100 (MQ=255)                                          aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                            <  1:1864760‑M2/42‑1 (MQ=255)                                          aaaaaaaaggcgactcatcagtcgcctgaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                            <  1:560022‑M2/42‑1 (MQ=255)                                              aaaaggcgactcatcagtcgccttaaaaatcagtttgccagtgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                        >  1:342701‑M2/55‑100 (MQ=255)                                               aaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAcgn                                                       <  1:1051555‑M2/47‑2 (MQ=255)                                               aaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                       <  1:852786‑M2/47‑1 (MQ=255)                                                 aggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                     <  1:1878997‑M2/49‑1 (MQ=255)                                                  ggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                    <  2:440495‑M2/50‑1 (MQ=255)                                                    cgactcatcagtcgccttaaaaatcagtttgcgagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                  >  1:2170503‑M2/49‑100 (MQ=255)                                                    cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                  <  2:878277‑M2/52‑1 (MQ=255)                                                        tcatcagtcgccttaaaaatcagtttgccagtgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                              <  2:342702‑M2/56‑1 (MQ=255)                                                         catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                             >  1:1891408‑M2/44‑100 (MQ=255)                                                          atcagtcgccttaaaaatcagtttgcgagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                            <  2:2170515‑M2/58‑1 (MQ=255)                                                          atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTCATCACTACGTTCATGTCGGtc                                            <  2:1732816‑M2/58‑1 (MQ=255)                                                          atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                            >  2:215737‑M2/43‑100 (MQ=255)                                                           tcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                           >  1:1821235‑M2/42‑100 (MQ=255)                                                            cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                          >  2:1833732‑M2/41‑100 (MQ=255)                                                               tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                       <  2:1821245‑M2/63‑1 (MQ=255)                                                               tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                       <  1:244268‑M2/63‑1 (MQ=255)                                                               tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                       <  2:577363‑M2/63‑1 (MQ=255)                                                                 gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                     >  2:468412‑M2/36‑100 (MQ=255)                                                                   cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                   >  2:58296‑M2/34‑100 (MQ=255)                                                                   cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                   <  1:272540‑M2/67‑1 (MQ=255)                                                                    ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                  >  1:1156096‑M2/33‑100 (MQ=255)                                                                     taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                 >  2:919755‑M2/32‑100 (MQ=255)                                                                      aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCagag                                >  1:343250‑M2/31‑100 (MQ=255)                                                                       aaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                               >  1:1444666‑M2/30‑100 (MQ=255)                                                                        aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                              >  2:18290‑M2/29‑100 (MQ=255)                                                                        aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                              <  2:499627‑M2/72‑1 (MQ=255)                                                                        aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGTCATCACTACGTTCATGTCCGTCTCTGCGGCACAGTc                              >  1:1886562‑M2/29‑100 (MQ=255)                                                                         tccgatcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                             <  2:1575473‑M2/73‑1 (MQ=255)                                                                           tcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCAGAGTCTTc                           <  1:1991779‑M2/75‑1 (MQ=255)                                                                            cagtttgcctgcgccgccttctgcCGTCCCCTGCACTTCAAGGGTGCGCCCGTCTTCGGTCATCACTACGTTCAGGTCGCTCTCTGCGGCAGAGTCTTc                           <  1:415368‑M2/75‑1 (MQ=255)                                                                             agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                         >  2:919996‑M2/24‑100 (MQ=255)                                                                              gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                        <  2:1925121‑M2/78‑1 (MQ=255)                                                                               tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                       <  2:1891418‑M2/79‑1 (MQ=255)                                                                               tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGACGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                       <  2:2216526‑M2/79‑1 (MQ=255)                                                                                 tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGgcagagtctagatcgga                     >  1:1575465‑M2/20‑92 (MQ=255)                                                                                     agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                 >  2:2043368‑M2/16‑100 (MQ=255)                                                                                     agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                 >  1:1786733‑M2/16‑100 (MQ=255)                                                                                     agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                 <  1:1833722‑M2/85‑1 (MQ=255)                                                                                        gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa              >  2:1338216‑M2/13‑100 (MQ=255)                                                                                         ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt             <  1:1023211‑M2/89‑1 (MQ=255)                                                                                         ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt             >  1:70378‑M2/12‑100 (MQ=255)                                                                                          cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc            <  1:1504894‑M2/90‑1 (MQ=255)                                                                                               tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa       >  2:1265889‑M2/6‑100 (MQ=255)                                                                                               tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa       >  2:1517124‑M2/6‑100 (MQ=255)                                                                                                ctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc      <  2:1331297‑M2/96‑1 (MQ=255)                                                                                                 tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc     <  2:1259739‑M2/97‑1 (MQ=255)                                                                                                  gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg    >  1:661789‑M2/3‑100 (MQ=255)                                                                                                  gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg    >  2:1340513‑M2/3‑100 (MQ=255)                                                                                                    cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt  >  2:1847871/1‑100 (MQ=255)                                                                                                    |                                                                                                    CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCT  >  NC_000913/3815842‑3816040

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