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

TGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCG  >  NC_000913/366025‑366510                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             gagtaacaacccgtcggattctccgtgggaacaaacggcggattgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccac                                                                                                                                                                                                                                             >  1:222464‑M2/250‑251 (MQ=255)    taacaacccgtcggattctccgtgggaacaaacggcggattgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccacgac                                                                                                                                                                                                                                          >  1:198521‑M2/247‑251 (MQ=255)            cgtcggattctccgggggaacaaccgggggattgaccgtaaagggataggtcacgttggtgtagatgggcgcatcgtacccttgcatcggccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTg                                                                                                                                                                                                                                  <  2:198521‑M2/13‑1 (MQ=255)                ggattctccgtgggaacaaacggcggattgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGGGAATc                                                                                                                                                                                                                              >  1:165156‑M2/235‑251 (MQ=255)                                                           aggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                                                                                                                                   <  1:124669‑M2/60‑1 (MQ=255)                                                                                    atcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                                                                                                                                           >  1:249463‑M2/167‑250 (MQ=255)                                                                                                      ccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCCTACAgtgt                                                                                                                                        >  2:232364‑M2/149‑251 (MQ=255)                                                                                                             gaggggacgacgacagtagcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTTACAATTCCACCCAACATACGAGCCGGAAGCATAAAGTGTAAAGcc                                                                                                                                  <  1:107805‑M2/109‑1 (MQ=255)                                                                                                              aggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttacggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggccgcgcaactgttgggaagggcgatcagtgcgggccACGACGGCCAGTGAATCCTGAATCATGGTCATAGCTGTATCCTGTGTGCAATTGTTATCCGCTCACCATTCCACACAAAATAAGAGCCGGAAGAATAAAGCGTAAAGCCTg                                                                                                                                >  1:197035‑M2/141‑251 (MQ=255)                                                                                                                    cgacgacagaatcggcctcaggaagatcgcactccagccagctttccggcacctcttctggtgccggaaacgaggcaaagcgccattcgcctttcatgcggcgcagctgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCTTAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATCCGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGc                                                                                                                          <  1:41428‑M2/117‑1 (MQ=255)                                                                                                                                                  actccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAACGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAAttgc                                                                                            >  1:210769‑M2/105‑251 (MQ=255)                                                                                                                                                                        cgcttctcgagccgaaaaccaggcaaagcgccattcgccattcaggctgcgcgactgttgggaagggcgatcggtgcgggcaACGACGGCCAGTGAATCCGCAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCAGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGTGTGCCTAATGAGTGAGCTAAATCACATTAATTTCGTTGTGCTCACTGCCCGCATTc                                                                      >  2:124202‑M2/83‑251 (MQ=255)                                                                                                                                                                                        aaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGt                                                      <  1:224945‑M2/185‑1 (MQ=255)                                                                                                                                                                                                   gctcctttcgccatcccggcggcgcacctgttgggaagggcgatcgtttctggccACGTCGGCCAGTGAATCCGTAATCATTGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACTCAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGATACCTGTCGTGCCAGCTg                                           <  2:90329‑M2/196‑1 (MQ=255)                                                                                                                                                                                                            cccattcaggctgcgcaactgttgggaagggcgatcggggcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGa                                  <  2:222464‑M2/205‑1 (MQ=255)                                                                                                                                                                                                                       ggcgcaactggggggaagggcgatcggggcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAcg                       <  2:57519‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                       cgcgatactggtgggaagggcgatcggtgcggtccACTACGGCCAGTGAATCCGTAATCATTGTCATAGCTGGTTCCTGTGTGAAATTGTTTACCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTATCTCACATTAATTGCGTTGCGCGCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAGCGGCCAAcg                       <  2:151611‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                               tgttgggaaggggggttggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCgggg                <  2:76994‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                                            gatcggtgcgggccacgacgGCCAGTGACACCGTAATCATGGTCATAGCTGATTCCTGTGTGAAACTGGTATACGCACACAATTCCAAACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAAATCACATTAATGGCGATGCGCTCACTGCCGGCGTTCCAGTCGGGAAAACGGTCGTGCAAGCTGCCTTAATGAATCGGCAAACGCGGCGGGAGAggagcagttcg  >  2:135239/15‑242 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             TGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCG  >  NC_000913/366025‑366510

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