Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	366,181	Δ93 bp	coding (33‑125/3075 nt)	lacZ ←	beta‑galactosidase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	366181	366273	93	18 [0]	[0] 19	lacZ	beta‑galactosidase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 366180	0 (0.000)	18 (0.700)	17/486	0.5	100%	coding (126/3075 nt)	lacZ	beta‑galactosidase
	?	NC_000913	366274 =	0 (0.000)					coding (32/3075 nt)	lacZ	beta‑galactosidase

GAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGC  >  NC_000913/365932‑366403                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                gAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCac                                                                                                                                                                                                                               >  1:222464‑M1/1‑249 (MQ=255)    tAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgac                                                                                                                                                                                                                            >  1:198521‑M1/1‑246 (MQ=255)            cGTCGGATTCTCCGGGGGAACAACCGGGGGATTGACCGTAAAGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTACCCTTGCATCGGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtg                                                                                                                                                                                                                    <  2:198521‑M1/251‑14 (MQ=255)                ggATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagggaatc                                                                                                                                                                                                                >  1:165156‑M1/1‑234 (MQ=255)                                                           aGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatcc                                                                                                                                                                     <  1:124669‑M1/251‑61 (MQ=255)                                                                                    aTCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                                                                                                                                             >  1:249463‑M1/1‑166 (MQ=255)                                                                                                      ccAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcctacagtgt                                                                                                                          >  2:232364‑M1/1‑148 (MQ=255)                                                                                                             gAGGGGACGACGACAGTAGCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgcttacaattccacccaacatacgagccggaagcataaagtgtaaagcc                                                                                                                    <  1:107805‑M1/250‑110 (MQ=255)                                                                                                              aGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTACGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCCGCGCAACTGTTGGGAAGGGCGATCAGTGCGGGCCacgacggccagtgaatcctgaatcatggtcatagctgtatcctgtgtgcaattgttatccgctcaccattccacacaaaataagagccggaagaataaagcgtaaagcctg                                                                                                                  >  1:197035‑M1/1‑140 (MQ=255)                                                                                                                    cgacgacAGAATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCTCTTCTGGTGCCGGAAACGAGGCAAAGCGCCATTCGCCTTTCATGCGGCGCAGCTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcttagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatccgagccggaagcataaagtgtaaagcctggggtgc                                                                                                            <  1:41428‑M1/251‑118 (MQ=255)                                                                                                                                                  aCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaacgcctggggtgcctaatgagtgagctaactcacattaattgc                                                                              >  1:210769‑M1/1‑104 (MQ=255)                                                                                                                                                                        cGCTTCTCGAGCCGAAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCGACTGTTGGGAAGGGCGATCGGTGCGGGCAacgacggccagtgaatccgcaatcatggtcatagctgtttcctgtgtgaaattgttatcagctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggtgtgcctaatgagtgagctaaatcacattaatttcgttgtgctcactgcccgcattc                                                        >  2:124202‑M1/1‑82 (MQ=255)                                                                                                                                                                                        aaCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgt                                        <  1:224945‑M1/251‑186 (MQ=255)                                                                                                                                                                                                   gctcctttcgccatcccggcggCGCACCTGTTGGGAAGGGCGATCGTTTCTGGCCacgtcggccagtgaatccgtaatcattgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccactcaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggatacctgtcgtgccagctg                             <  2:90329‑M1/245‑197 (MQ=255)                                                                                                                                                                                                            cccattcaGGCTGCGCAACTGTTGGGAAGGGCGATCGGGGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatga                    <  2:222464‑M1/250‑206 (MQ=255)                                                                                                                                                                                                                       ggcgcaACTGGGGGGAAGGGCGATCGGGGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacg         <  2:57519‑M1/250‑217 (MQ=255)                                                                                                                                                                                                                       cgcgatactggtgGGAAGGGCGATCGGTGCGGTCCactacggccagtgaatccgtaatcattgtcatagctggttcctgtgtgaaattgtttaccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctatctcacattaattgcgttgcgcgcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaagcggccaacg         <  2:151611‑M1/245‑217 (MQ=255)                                                                                                                                                                                                                               tGTTGGGAAGGGGGGTTGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcgggg  <  2:76994‑M1/250‑224 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                GAGTAACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGC  >  NC_000913/365932‑366403

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