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	65 [0]	[0] 66	[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)	60 (1.020)	47/178	0.2	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                                                                                                    |                                                                                                    tggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                     >  1:334721‑M2/100‑100 (MQ=255)   gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                   >  2:1094611‑M2/98‑100 (MQ=255)   gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                   >  2:408601‑M2/98‑100 (MQ=255)     ttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                 <  2:460944‑M2/5‑1 (MQ=255)     ttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                 <  2:327653‑M2/5‑1 (MQ=255)      tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                >  2:311482‑M2/95‑100 (MQ=255)          gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGc                                                                                            >  1:549426‑M2/91‑100 (MQ=255)            cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                          <  2:1302986‑M2/12‑1 (MQ=255)             ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagagccgccttctgcCGTCCCCTGCACt                                                                                         >  1:1466719‑M2/88‑100 (MQ=255)              tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                        <  2:1315549‑M2/14‑1 (MQ=255)                gctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTc                                                                                       <  1:1195047‑M2/15‑1 (MQ=255)                  tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                    <  2:994260‑M2/18‑1 (MQ=255)                   cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                   >  2:828731‑M2/82‑100 (MQ=255)                     tcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                 >  2:1581621‑M2/80‑100 (MQ=255)                            acttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCg                                                                          >  1:925647‑M2/73‑100 (MQ=255)                                  ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                     <  2:549430‑M2/33‑1 (MQ=255)                                  ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                    <  1:793634‑M2/34‑1 (MQ=255)                                    acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                  <  1:582786‑M2/36‑1 (MQ=255)                                       gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                 <  1:221951‑M2/37‑1 (MQ=255)                                         caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             >  1:1513908‑M2/60‑100 (MQ=255)                                         caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                             <  2:76176‑M2/41‑1 (MQ=255)                                          caaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                            >  1:1587132‑M2/59‑100 (MQ=255)                                          aaaaaaaaggcgactcattagttgccttaaaaatcaggttgccagcgacgccttctgcCGTCCCATGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                               >  1:61649‑M2/59‑97 (MQ=255)                                           aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                           <  1:251541‑M2/43‑1 (MQ=255)                                           aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                           >  2:1129601‑M2/58‑100 (MQ=255)                                            aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                          >  2:284273‑M2/57‑100 (MQ=255)                                            aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                          >  2:703906‑M2/57‑100 (MQ=255)                                             tttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                         >  2:1004528‑M2/56‑100 (MQ=255)                                             aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGTGCCCGTCTTCGGTCATCACTAc                                                         >  1:1368731‑M2/56‑100 (MQ=255)                                              ttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                        >  2:589576‑M2/55‑100 (MQ=255)                                              aaaaggcgactcatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                        <  1:1316467‑M2/46‑1 (MQ=255)                                               taaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                       >  1:485741‑M2/54‑100 (MQ=255)                                                  ggcgactcatcagtctccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                    <  2:1576484‑M2/50‑1 (MQ=255)                                                  ggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                    <  1:321074‑M2/50‑1 (MQ=255)                                                    cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                  >  1:94933‑M2/49‑100 (MQ=255)                                                        tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                              <  2:185653‑M2/56‑1 (MQ=255)                                                         catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCAGTCTTCGGTCATCACTACGTTCATGTCGGt                                             >  2:147120‑M2/44‑100 (MQ=255)                                                         catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTCCACTTCAATGATGTGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                             >  1:715809‑M2/44‑100 (MQ=255)                                                          atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                            <  2:424621‑M2/58‑1 (MQ=255)                                                          atcagtcgccttaaaaatcagtgtgccagcgccgccttctggCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                            <  1:105430‑M2/58‑1 (MQ=255)                                                            cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                          >  1:554403‑M2/41‑100 (MQ=255)                                                            cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                          <  1:449769‑M2/60‑1 (MQ=255)                                                               tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                       <  2:451947‑M2/63‑1 (MQ=255)                                                                  ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                    <  1:512419‑M2/66‑1 (MQ=255)                                                                  ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                    >  2:1456521‑M2/35‑100 (MQ=255)                                                                   cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                   >  2:1411800‑M2/34‑100 (MQ=255)                                                                     taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                 >  1:1132267‑M2/32‑100 (MQ=255)                                                                     taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGc                                    >  1:154087‑M2/32‑97 (MQ=255)                                                                      aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCagag                                >  2:414772‑M2/31‑100 (MQ=255)                                                                        aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                              >  2:274922‑M2/29‑100 (MQ=255)                                                                         aatcagtgtgcaagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCAGCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                             >  1:840194‑M2/28‑100 (MQ=255)                                                                               tttgccagctccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGTTCTCTGCGGCAGAGTCTTCAACg                       <  1:1349163‑M2/79‑1 (MQ=255)                                                                               tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                       <  1:239370‑M2/79‑1 (MQ=255)                                                                                 tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                     >  2:1145156‑M2/20‑100 (MQ=255)                                                                                   ccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTTGGTCTCTGCGGCAGAGTCTTCAACGTAtt                   <  2:1247529‑M2/83‑1 (MQ=255)                                                                                       cgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg               <  2:554407‑M2/87‑1 (MQ=255)                                                                                        gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa              <  1:974482‑M2/88‑1 (MQ=255)                                                                                        gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa              >  1:1203122‑M2/13‑100 (MQ=255)                                                                                         ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt             >  2:1446326‑M2/12‑100 (MQ=255)                                                                                         ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt             >  1:247850‑M2/12‑100 (MQ=255)                                                                                          cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc            <  2:334724‑M2/90‑1 (MQ=255)                                                                                           gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg           <  1:1456506‑M2/91‑1 (MQ=255)                                                                                             cttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa         >  1:1195662‑M2/8‑100 (MQ=255)                                                                                                ctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa       >  1:723687‑M2/5‑99 (MQ=255)                                                                                                   cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAaactg    >  2:495409‑M2/2‑97 (MQ=255)                                                                                                    cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt  <  2:1405052/100‑1 (MQ=255)                                                                                                    |                                                                                                    CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCT  >  NC_000913/3815842‑3816040

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