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	26 [0]	[0] 26	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)	26 (1.100)	25/496	0.1	100%	coding (126/3075 nt)	lacZ	beta‑galactosidase
	?	NC_000913	366274 =	0 (0.000)					coding (32/3075 nt)	lacZ	beta‑galactosidase

AACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCT  >  NC_000913/365936‑366404                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 aacatcCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGTCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacg                                                                                                                                                                                                                            <  1:110170‑M1/251‑7 (MQ=255)             gATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatcc                                                                                                                                                                                                                <  1:35084‑M1/251‑19 (MQ=255)                          aaCAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGGGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCTGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtc                                                                                                                                                                                                    >  1:161705‑M1/1‑220 (MQ=255)                               cctggggatttcccgtcatgggatagGTCATGTTGGTGTGGATGGTCGCATCGTAACCGTGCATCTGCCAGTTGGATGGGAAGACGACAGTATCGGCCTCAGGACGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTTCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTTCGCAACTGTTTTGAAGGTCGATCGTTGCTGGCCacgacggccagtgaatccgtaatcatggtcatagct                                                                                                                                                                                              <  2:10473‑M1/239‑37 (MQ=255)                                      aTTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcct                                                                                                                                                                                       >  1:234398‑M1/1‑208 (MQ=255)                                         gACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtg                                                                                                                                                                                    >  1:175240‑M1/1‑205 (MQ=255)                                            cGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtg                                                                                                                                                                                  >  1:197642‑M1/1‑202 (MQ=255)                                                  gggATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattg                                                                                                                                                                            <  2:189300‑M1/250‑55 (MQ=255)                                                   ggATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgtt                                                                                                                                                                          <  1:151854‑M1/251‑57 (MQ=255)                                                    gATAGGTCACGTTGGCGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgtt                                                                                                                                                                          >  1:37601‑M1/1‑194 (MQ=255)                                                                    gtAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                                                                                                                                         <  2:51919‑M1/251‑74 (MQ=255)                                                                                   gTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCCCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagggaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagcc                                                                                                                                          <  2:118782‑M1/251‑89 (MQ=255)                                                                                                                cgacgacAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgc                                                                                                             >  1:91922‑M1/1‑134 (MQ=255)                                                                                                                                       aGATCGCCCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggctagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacat                                                                                      <  1:212891‑M1/251‑141 (MQ=255)                                                                                                                                         aTCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTACCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCCACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcatacagagtaaagcctggggtgcctaatgagtgagctaactcacatt                                                                                     >  2:121114‑M1/1‑109 (MQ=255)                                                                                                                                              cctCCAGCCATCTTTCCGGAACCGCTTCGGGGGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgtcggccagtgaatccgtaatcatggtcatagctgttgcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgc                                                                               <  2:124477‑M1/250‑148 (MQ=255)                                                                                                                                                     ccagcTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgc                                                                        >  1:153422‑M1/1‑97 (MQ=255)                                                                                                                                                               ggCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccg                                                              <  2:90808‑M1/251‑165 (MQ=255)                                                                                                                                                                                       cAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgt                                      <  1:109145‑M1/251‑189 (MQ=255)                                                                                                                                                                                                cgccattcgccattCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgatttccagtcgggaaacctgtcgtgccagctgc                             >  1:224359‑M1/1‑54 (MQ=255)                                                                                                                                                                                                           attcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatc                  >  2:198867‑M1/1‑43 (MQ=255)                                                                                                                                                                                                           attcAGGCTGCGCAACTGTTGGGAAGGGCGATCGGCGCGGGCCacgacggccagtaaatccgtcctcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcatcaagtgcaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctccctgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatc                  >  2:235804‑M1/1‑43 (MQ=255)                                                                                                                                                                                                              cAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgctttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggc               <  1:88007‑M1/251‑212 (MQ=255)                                                                                                                                                                                                                           tGTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcgggg   >  1:227654‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                            gTTGGGCAGGGCGACCGGTGCGGGCCacgccggccagcgaatccgtaatcatggtcatagatgtttactgtgtgaaattgttatctgctcacaattccacacaacacacgagccggacgcataaagtgtaaagcctggggagcctaatgagtgagctaactcacattacttgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctccattaatgaatcggccaacgcgcgggg   >  2:76498‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                            gTTGGGAAGGGCGATCGGTGCGGGCCacgacggccagtgaatccgtaatcatggtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaagtgtaaagcctggggtgcctaatgagtgagctaactcacattaattgcgttgcgctcactgcccgctttccagtcgggaaacctgtcgtgccagctgcattaatgaatcggccaacgcgcgggga  <  2:91149‑M1/250‑225 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 AACAACCCGTCGGATTCTCCGTGGGAACAAACGGCGGATTGACCGTAATGGGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCT  >  NC_000913/365936‑366404

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