Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	3,912,149	Δ15 bp	coding (28‑42/99 nt)	ilvL →	ilvG operon leader peptide

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	3912149	3912163	15	33 [0]	[0] 33	ilvL	ilvG operon leader peptide

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3912148	0 (0.000)	28 (1.080)	25/132	0.1	100%	coding (27/99 nt)	ilvL	ilvG operon leader peptide
	?	REL606	3912164 =	0 (0.000)					coding (43/99 nt)	ilvL	ilvG operon leader peptide

CCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAA  >  REL606/3912093‑3912238                                                                        |                                                                           aactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcgtggt                                                                        <  1:20276‑M2/5‑1 (MQ=255)  actctttaggccttcttacgaacaagatgccagaacagacaanangacaggcnttctacgagtgattagcgtggtg                                                                       <  1:897443‑M2/6‑1 (MQ=255)      tttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGa                                                                   >  2:555976‑M2/67‑76 (MQ=255)      tttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGa                                                                   >  1:899862‑M2/67‑76 (MQ=255)       ttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGat                                                                  <  2:719378‑M2/11‑1 (MQ=255)       ttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGat                                                                  >  1:608079‑M2/66‑76 (MQ=255)           gcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttcgacgagtgattagcGTGGTGGTGattatt                                                              >  2:435724‑M2/62‑76 (MQ=255)            cattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGAttatta                                                             >  1:410695‑M2/61‑76 (MQ=255)             attccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATtattat                                                            >  1:244618‑M2/60‑76 (MQ=255)                 cttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCa                                                        <  1:179548‑M2/21‑1 (MQ=255)                 cttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCa                                                        <  2:637596‑M2/21‑1 (MQ=255)                  ttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCAc                                                       <  1:513669‑M2/22‑1 (MQ=255)                    cgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCg                                                     <  1:169657‑M2/24‑1 (MQ=255)                      aacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTg                                                   <  1:555599‑M2/26‑1 (MQ=255)                           gatgcaagaaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCgggg                                              <  1:518249‑M2/31‑1 (MQ=255)                                   aaaagacaaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACtt                                      <  1:555976‑M2/39‑1 (MQ=255)                                         caaaatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGa                                <  2:114739‑M2/45‑1 (MQ=255)                                           aaatgacagcccttctacgagtgattagcGTGGCGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAgg                              >  2:322586‑M2/30‑76 (MQ=255)                                            aatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGa                             >  2:4500‑M2/29‑76 (MQ=255)                                            aatgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGa                             <  2:17195‑M2/48‑1 (MQ=255)                                              tgacagcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGaaa                           >  1:927547‑M2/27‑76 (MQ=255)                                                  agcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCt                       >  2:45247‑M2/23‑76 (MQ=255)                                                  agcccttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCt                       >  1:485496‑M2/23‑76 (MQ=255)                                                    cccttctacgagtgattagcGTGGTGGTGATTATTATGCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTa                     <  2:69635‑M2/56‑1 (MQ=255)                                                       ttctacgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTagag                  <  1:367987‑M2/59‑1 (MQ=255)                                                          tacgagtgattagcGTGGTGGTGATTACTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATc               >  2:93373‑M2/15‑76 (MQ=255)                                                           acgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCa              <  1:877789‑M2/63‑1 (MQ=255)                                                            cgagtgattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCaa             <  2:480298‑M2/64‑1 (MQ=255)                                                                 gattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCtt        >  2:271292‑M2/8‑76 (MQ=255)                                                                 gattagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCtt        >  2:259114‑M2/8‑76 (MQ=255)                                                                   ttagcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGGCGAGGAAAGGCTTAGAGATCAAGCCTTaa      <  2:740730‑M2/71‑1 (MQ=255)                                                                     agcGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACg    <  1:45247‑M2/73‑1 (MQ=255)                                                                       cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGaa  >  2:427965‑M2/2‑76 (MQ=255)                                                                        |                                                                           CCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAA  >  REL606/3912093‑3912238

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