Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,815,859	Δ82 bp		[rph]	[rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3815859	3815940	82	28 [4]	[0] 24	[rph]	[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	3 (0.060)	21 (0.490)	15/114	0.6	94.3%	intergenic (‑90/+5)	pyrE/rph	orotate phosphoribosyltransferase/truncated RNase PH
	?	NC_000913	3815941 =	0 (0.000)					coding (609/687 nt)	rph	truncated RNase PH

TTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGA  >  NC_000913/3815791‑3815924                                                                    |                                                                   ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                    <  1:1289962‑M1/68‑1 (MQ=255)   cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                  >  1:1164550‑M1/1‑66 (MQ=255)    tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                 <  1:463910‑M1/68‑4 (MQ=255)    tACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                 <  1:341100‑M1/68‑4 (MQ=255)      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                               <  1:3713009‑M1/68‑6 (MQ=255)      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                               <  1:1128112‑M1/68‑6 (MQ=255)      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                               <  1:2214496‑M1/68‑6 (MQ=255)      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                               <  1:600918/68‑1 (MQ=255)      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                               <  1:3558303/68‑1 (MQ=255)       aGAAAAAAAAAAGGCGACTCATCAGTGGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCt                                                              <  1:2990580/68‑1 (MQ=255)         aCAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                            <  1:955234‑M1/68‑9 (MQ=255)            aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                         <  1:715585‑M1/68‑12 (MQ=255)              aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                       <  1:1450345‑M1/68‑14 (MQ=255)              aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                       <  1:1622016‑M1/68‑14 (MQ=255)               aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                      >  1:2927642‑M1/1‑54 (MQ=255)                aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                     <  1:490560‑M1/68‑16 (MQ=255)                             cggtcGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                        <  1:2656905‑M1/66‑29 (MQ=255)                              aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                       <  1:1587305‑M1/68‑30 (MQ=255)                              aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                       <  1:3308745‑M1/68‑30 (MQ=255)                                tCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCTGCcgtcccctgcacttcaatgatgcgcccgtct                                     <  1:2659722‑M1/68‑32 (MQ=255)                                         aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                            <  1:3052525‑M1/68‑41 (MQ=255)                                          aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                           >  1:430384‑M1/1‑27 (MQ=255)                                               gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                      <  1:1334530‑M1/68‑47 (MQ=255)                                                tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                     <  1:744319‑M1/68‑48 (MQ=255)                                                         gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg            <  1:2798245‑M1/68‑57 (MQ=255)                                                           cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc          <  1:2067717‑M1/68‑59 (MQ=255)                                                            gccCTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct         <  1:506182‑M1/68‑60 (MQ=255)                                                                   gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggc  <  1:1713148‑M1/68‑67 (MQ=255)                                                                    |                                                                   TTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGA  >  NC_000913/3815791‑3815924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑