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

TTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGC  >  NC_000913/2877927‑2878405                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         atatctcaccttcgttgccggtttatccctgctggcgcggggaactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacgttcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGt                                                                                                                                                                                                                                      <  2:123542‑M2/21‑1 (MQ=255)  aatctctgcttcgttgccggtttatccctgctggggcggtgaaacctcggttcaggcgttgcaaacctggcttccggactgtgtatccccgctcacgcggggaacgcatagtacatcattccgcctatgtctgaactcgcggtttatccccgatggcgcggggaactcccgggggataatgtttacggtcatgagccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTa                                                                                                                                                                                                                                     <  2:210554‑M2/22‑1 (MQ=255)         ccgtcgttgccggtttatccctgctggcgcggggaactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGa                                                                                                                                                                                                                              <  1:109071‑M2/29‑1 (MQ=255)                         atccctgctggcgcggggaactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATccc                                                                                                                                                                                                               <  1:122234‑M2/44‑1 (MQ=255)                         atccctgctggcgcggggaactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATccc                                                                                                                                                                                                               <  1:497167‑M2/44‑1 (MQ=255)                            cctgcgggcgcggggaactctcggttcaggcgtggcaatccgggctaccgggcggtttatccccgctaacgcggggaactcggagtccatcattccaccgatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggagaatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGTGGTTTATCCCCGCt                                                                                                                                                                                                           <  2:67630‑M2/48‑1 (MQ=255)                                          gaactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                                                                                                             <  1:459643‑M2/62‑1 (MQ=255)                                           aactctcggttcaggcgttgcaaacctggctaccgggcggtttatccccgctaacgcggggaactcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTcg                                                                                                                                                                                            >  2:333644‑M2/189‑251 (MQ=255)                                                      caggcgttgcaaacctggctaccgagcggtttatccccgctaacgcggggaactcgtcgtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgagccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTc                                                                                                                                                                                   >  1:197901‑M2/178‑249 (MQ=255)                                                        ggcgtggcaaacctggctcccgggcaggttatccccggtgacgcggggaactcgtagtccatcattcctcctatgtctgcactcccggtttatccccgctggcggggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGaa                                                                                                                                                                               <  2:115165‑M2/76‑1 (MQ=255)                                                                             gggcggtgtatccccgctaacgcggggaactcgtagtccatcattccacctatgtcggaactcccggtttatccccgcgggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACg                                                                                                                                                          <  2:61959‑M2/97‑1 (MQ=255)                                                                                 ggtttatccccgctaacgcggggaactcgtantccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGttt                                                                                                                                                      <  2:225415‑M2/101‑1 (MQ=255)                                                                                    ttatccccgctaacgcggggaattcgtagtccatcattccacctatgtctgaacgcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTAt                                                                                                                                                    <  2:108215‑M2/103‑1 (MQ=255)                                                                                                      ggaacgcgtagtccatcattccacctatgtctgaactcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAAc                                                                                                                                 <  1:294144‑M2/122‑1 (MQ=255)                                                                                                                                        actcccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGcc                                                                                               <  1:317953‑M2/156‑1 (MQ=255)                                                                                                                                           cccggtttatccccgctggcgcggggaactcccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGCATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGccc                                                                                              <  1:417777‑M2/157‑1 (MQ=255)                                                                                                                                                   atccccgctggcgcggggaactcccgggggataatgtttacggtcaagcgccccccggtttatccccgctggcgcggggaacgcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCTCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGTAGCCCGGTTTATcc                                                                                     <  1:78646‑M2/166‑1 (MQ=255)                                                                                                                                                               gcggggaactcctggggtataatgttgacggtcattcgccccccgtttgatccccgctggcgcggggaactcGCTGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGTCTCGGGGAACTCGCGACCGCTCAGAAATTCCAGTCCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATGATCAATTACAACCGACAGGGAGCCCTGGTTATCACCGCTGGCGCggg                                                                        <  2:38290‑M2/179‑1 (MQ=255)                                                                                                                                                                    gacctcccgggggctaatgtgtacggtcaggcgccccccggtttatccccgctggcgcggggaacgcGCAGGCGGCGACGGGCAGGGTATGCTGGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATGATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACt                                                                   <  2:433934‑M2/184‑1 (MQ=255)                                                                                                                                                                         cccgggggataatgtttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAg                                                                <  1:119596‑M2/187‑1 (MQ=255)                                                                                                                                                                                     ggtttacggtcatgcgccccccggtttagccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCa                                                  <  2:320935‑M2/201‑1 (MQ=255)                                                                                                                                                                                        ttacggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCAcc                                                <  1:15966‑M2/203‑1 (MQ=255)                                                                                                                                                                                           cggtcatgcgccccccggtttatccccgctggcgcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTCTCCCCGCTGGCGCGGGGCACTCGCGACCGCTCAGAAATTACAGACCCGAACCAAACGGTTTATCCCCGCGGACGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGGACCCGGTCTATCCCCCCCGGCGCGGGGAACTCAGCGTGTTCTGCAACACctttt                                            >  2:275889‑M2/45‑250 (MQ=255)                                                                                                                                                                                                  gcgccccccggttttgccccggtggggcggggaactcGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCgg                                     <  2:164663‑M2/214‑1 (MQ=255)                                                                                                                                                                                                              ttttccctgcgggcggcgggtaCTCGCCGGCGGCGAGGCGTAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGGGCGGGGAACGCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCGGGCGCGGGGACCCCTCAACATTAGCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTAt                          <  2:8559/229‑1 (MQ=255)                                                                                                                                                                                                                     cGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGCACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCACTTAAAACCGACAGGGAGCCCGGGTTATCCCCGCTGGCGCGGGGAACTCCGCGTGTTCCGCATCCCCTTTGACTTCGACTGCGGTTTATCCCCGCTg                  >  2:376565/1‑251 (MQ=255)                                                                                                                                                                                                                      gCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCCGGCGCGGGGAACTCTCAACATTATCAATTACACCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTATTCGGCATAACTTTTGGCTCCGGCTCCGGTTCATCCCCGccgg                 >  2:115994/1‑248 (MQ=255)                                                                                                                                                                                                                                      cGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGc  <  1:153658/250‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         TTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGC  >  NC_000913/2877927‑2878405

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