Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,877,974	Δ183 bp	intergenic (+356/+413)	iap → / ← cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2877974	2878156	183	24 [0]	[4] 28	iap/cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2877973	0 (0.000)	24 (0.420)	23/420	0.8	100%	intergenic (+355/+596)	iap/cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2
	?	NC_000913	2878157 =	0 (0.000)					intergenic (+539/+412)	iap/cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

AAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATT  >  NC_000913/2877744‑2878187                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       atatcTCACCTTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGTTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggt                                                                                                                                                                                                   <  2:123542‑M1/249‑22 (MQ=255)  aatctctgcttcgttgccggtttaTCCCTGCTGGGGCGGTGAAACCTCGGTTCAGGCGTTGCAAACCTGGCTTCCGGACTGTGTATCCCCGCTCACGCGGGGAACGCATAGTACATCATTCCGCCTATGTCTGAACTCGCGGTTTATCCCCGATGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGAGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggta                                                                                                                                                                                                  <  2:210554‑M1/241‑23 (MQ=255)         ccGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcga                                                                                                                                                                                           <  1:109071‑M1/251‑30 (MQ=255)                         aTCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccc                                                                                                                                                                            <  1:122234‑M1/250‑45 (MQ=255)                         aTCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccc                                                                                                                                                                            <  1:497167‑M1/250‑45 (MQ=255)                            ccTGCGGGCGCGGGGAACTCTCGGTTCAGGCGTGGCAATCCGGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGGAGTCCATCATTCCACCGATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGAGAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgtggtttatccccgct                                                                                                                                                                        <  2:67630‑M1/251‑49 (MQ=255)                                          gAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactc                                                                                                                                                          <  1:459643‑M1/251‑63 (MQ=255)                                           aaCTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcg                                                                                                                                                         >  2:333644‑M1/1‑188 (MQ=255)                                                      cAGGCGTTGCAAACCTGGCTACCGAGCGGTTTATCCCCGCTAACGCGGGGAACTCGTCGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGAGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctc                                                                                                                                                >  1:197901‑M1/1‑177 (MQ=255)                                                        ggCGTGGCAAACCTGGCTCCCGGGCAGGTTATCCCCGGTGACGCGGGGAACTCGTAGTCCATCATTCCTCCTATGTCTGCACTCCCGGTTTATCCCCGCTGGCGGGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaa                                                                                                                                            <  2:115165‑M1/251‑77 (MQ=255)                                                                             gggCGGTGTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCGGAACTCCCGGTTTATCCCCGCGGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacg                                                                                                                       <  2:61959‑M1/251‑98 (MQ=255)                                                                                 ggTTTATCCCCGCTAACGCGGGGAACTCGTANTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggttt                                                                                                                   <  2:225415‑M1/251‑102 (MQ=255)                                                                                    ttATCCCCGCTAACGCGGGGAATTCGTAGTCCATCATTCCACCTATGTCTGAACGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttat                                                                                                                 <  2:108215‑M1/250‑104 (MQ=255)                                                                                                      ggAACGCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaac                                                                                              <  1:294144‑M1/251‑123 (MQ=255)                                                                                                                                        aCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagcc                                                            <  1:317953‑M1/251‑157 (MQ=255)                                                                                                                                           cccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggcatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagccc                                                           <  1:417777‑M1/249‑158 (MQ=255)                                                                                                                                                   aTCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCAAGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACGCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatctccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacaggtagcccggtttatcc                                                  <  1:78646‑M1/250‑167 (MQ=255)                                                                                                                                                               gcgGGGAACTCCTGGGGTATAATGTTGACGGTCATTCGCCCCCCGTTTGATCCCCGCTGGCGCGGGGAACTCgctggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctgtctcggggaactcgcgaccgctcagaaattccagtcccgatccaaacggtttatccccgctggcgcggggaactctcaacatgatcaattacaaccgacagggagccctggttatcaccgctggcgcggg                                     <  2:38290‑M1/251‑180 (MQ=255)                                                                                                                                                                    gacctccCGGGGGCTAATGTGTACGGTCAGGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACGCgcaggcggcgacgggcagggtatgctggattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacatgatcaattacaaccgacagggagcccggtttatccccgctggcgcggggaact                                <  2:433934‑M1/248‑185 (MQ=255)                                                                                                                                                                         cccGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagcccggtttatccccgctggcgcggggaactcag                             <  1:119596‑M1/249‑188 (MQ=255)                                                                                                                                                                                     ggtTTACGGTCATGCGCCCCCCGGTTTAGCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagcccggtttatccccgctggcgcggggaactcagcgtgttcggcatca               <  2:320935‑M1/250‑202 (MQ=255)                                                                                                                                                                                        ttACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagcccggtttatccccgctggcgcggggaactcagcgtgttcggcatcacc             <  1:15966‑M1/250‑204 (MQ=255)                                                                                                                                                                                           cGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttctccccgctggcgcggggcactcgcgaccgctcagaaattacagacccgaaccaaacggtttatccccgcggacgcggggaactctcaacattatcaattacaaccgacaggggacccggtctatcccccccggcgcggggaactcagcgtgttctgcaacacctttt         >  2:275889‑M1/1‑44 (MQ=255)                                                                                                                                                                                                  gcgcCCCCCGGTTTTGCCCCGGTGGGGCGGGGAACTCgcaggcggcgacgcgcagggtatgcgcgattcgcggtttatccccgctggcgcggggaactcgcgaccgctcagaaattccagacccgatccaaacggtttatccccgctggcgcggggaactctcaacattatcaattacaaccgacagggagcccggtttatccccgctggcgcggggaactcagcgtgttcggcatcacctttggcttcgg  <  2:164663‑M1/251‑215 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       AAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATT  >  NC_000913/2877744‑2878187

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