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	70 [1]	[0] 70	[rph]–[rph]	[rph], [rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	1 (0.010)	65 (1.030)	52/178	0.2	99.3%	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

AGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCT  >  NC_000913/3815844‑3816040                                                                                                  |                                                                                                    gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                   >  2:163601‑M2/98‑100 (MQ=255)     tctgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                               <  2:162385‑M2/7‑1 (MQ=255)     catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                               >  1:593222‑M2/94‑100 (MQ=255)     catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                               <  2:378602‑M2/7‑1 (MQ=255)     catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                               <  1:775610‑M2/7‑1 (MQ=255)       tgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTg                                                                                             >  2:1455537‑M2/92‑100 (MQ=255)        gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgccgtccccaga                                                                                            >  1:162385‑M2/91‑97 (MQ=255)         ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                           <  2:218103‑M2/11‑1 (MQ=255)         ccttcgctcctcatcttactttaaaaaaaaaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                           <  1:1248606‑M2/11‑1 (MQ=255)          cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                          >  1:938336‑M2/89‑100 (MQ=255)            tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                        >  2:1143982‑M2/87‑100 (MQ=255)            tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACt                                                                                         >  1:29601‑M2/87‑99 (MQ=255)                tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                    >  2:7065‑M2/83‑100 (MQ=255)                 cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                   <  1:940341‑M2/19‑1 (MQ=255)                 cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                   >  1:28054‑M2/82‑100 (MQ=255)                  ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                  >  1:773974‑M2/81‑100 (MQ=255)                  ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                  <  1:201029‑M2/20‑1 (MQ=255)                    catcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                >  1:641126‑M2/79‑100 (MQ=255)                     atcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgc                                                                               >  1:1691068‑M2/78‑100 (MQ=255)                     aacttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgc                                                                               >  2:1148281‑M2/78‑100 (MQ=255)                      tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgc                                                                               <  1:447630‑M2/23‑1 (MQ=255)                        ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                            <  1:116598‑M2/26‑1 (MQ=255)                         tacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                           <  2:593224‑M2/27‑1 (MQ=255)                           cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGcgcccat                                                                         >  1:1036553‑M2/72‑98 (MQ=255)                           cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCg                                                                          >  1:632922‑M2/72‑99 (MQ=255)                           cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctccCGTCCCCTTCACTTCAATGATGCTCCCg                                                                          >  1:1090137‑M2/72‑99 (MQ=255)                               tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                     <  1:1231546‑M2/33‑1 (MQ=255)                                ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                    >  1:626111‑M2/67‑100 (MQ=255)                                 tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                   <  2:1623010‑M2/35‑1 (MQ=255)                                  acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                  <  2:201492‑M2/36‑1 (MQ=255)                                   cagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                 <  1:1514174‑M2/37‑1 (MQ=255)                                    agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                <  2:1029749‑M2/38‑1 (MQ=255)                                    agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                 >  1:430851‑M2/63‑99 (MQ=255)                                    agacaaaaaaaaggcgactcatcagtcgccgtaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                >  1:1117022‑M2/63‑100 (MQ=255)                                     ctgatttttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                               >  1:1197866‑M2/62‑100 (MQ=255)                                       caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             >  2:1367870‑M2/60‑100 (MQ=255)                                       caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             >  2:1044170‑M2/60‑100 (MQ=255)                                        aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                            >  1:1464640‑M2/59‑100 (MQ=255)                                         tttttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                           >  2:514304‑M2/58‑100 (MQ=255)                                         aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                           >  1:1269813‑M2/58‑100 (MQ=255)                                            aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                        <  1:982819‑M2/46‑1 (MQ=255)                                                ggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                    >  2:1662142‑M2/51‑100 (MQ=255)                                                  cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                  >  1:409020‑M2/49‑100 (MQ=255)                                                   gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                 <  2:362910‑M2/53‑1 (MQ=255)                                                    actcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                >  1:1360357‑M2/47‑100 (MQ=255)                                                     ctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                               >  1:1390096‑M2/46‑100 (MQ=255)                                                     ctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                               >  2:1461922‑M2/46‑100 (MQ=255)                                                           agtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                         >  1:740409‑M2/40‑100 (MQ=255)                                                            gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                        >  2:764942‑M2/39‑100 (MQ=255)                                                            gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                        <  2:824203‑M2/62‑1 (MQ=255)                                                             tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                       <  2:1538545‑M2/63‑1 (MQ=255)                                                              cgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                      <  1:781148‑M2/64‑1 (MQ=255)                                                                 cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCAGTCTCTGCGGc                                    >  1:817057‑M2/34‑99 (MQ=255)                                                                 cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                   >  2:712785‑M2/34‑100 (MQ=255)                                                                   taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                 >  1:392186‑M2/32‑100 (MQ=255)                                                                    aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCagag                                <  2:441526‑M2/70‑1 (MQ=255)                                                                     aaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                               >  1:519886‑M2/30‑100 (MQ=255)                                                                      aaatcagtttgccagcgccgccttctgcCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTAATCACTACGTTCATGTCGGTATCTCCGgcagacc                               >  1:1385229‑M2/29‑97 (MQ=255)                                                                          cagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                          <  1:1691079‑M2/76‑1 (MQ=255)                                                                           agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                         >  2:288782‑M2/24‑100 (MQ=255)                                                                           agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                         <  2:626113‑M2/77‑1 (MQ=255)                                                                             tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                       <  2:1691075‑M2/79‑1 (MQ=255)                                                                                gccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                    >  2:1676218‑M2/19‑100 (MQ=255)                                                                                    gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                <  1:255453‑M2/86‑1 (MQ=255)                                                                                     cgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg               >  2:1454114‑M2/14‑100 (MQ=255)                                                                                        cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACATATTCCAGATc            <  2:1692999‑M2/90‑1 (MQ=255)                                                                                         gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg           <  2:1464647‑M2/91‑1 (MQ=255)                                                                                            ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACATATTCCAGATCGCaa        >  1:607192‑M2/7‑100 (MQ=255)                                                                                              ctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc      <  2:1390103‑M2/96‑1 (MQ=255)                                                                                                  cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt  <  1:621920/100‑1 (MQ=255)                                                                                                  |                                                                                                    AGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCT  >  NC_000913/3815844‑3816040

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