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

TCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGA  >  NC_000913/366029‑366497                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                aacatcccgtcggattctccgtgggaacaaacggcggattgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcagtctgcgcaactgttgggaagggcgatcggtgcgggccacgacg                                                                                                                                                                                                                            <  1:110170‑M2/6‑1 (MQ=255)             gattctccgtgggaacaaacggcggattgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATcc                                                                                                                                                                                                                <  1:35084‑M2/18‑1 (MQ=255)                          aacaaacggcggattgaccgtaatgggataggtcacgttggggtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttcctgcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTc                                                                                                                                                                                                    >  1:161705‑M2/221‑250 (MQ=255)                               cctggggatttcccgtcatgggataggtcatgttggtgtggatggtcgcatcgtaaccgtgcatctgccagttggatgggaagacgacagtatcggcctcaggacgatcgcactccagccagctttccggcaccgcttctggttccggaaaccaggcaaagcgccattcgccattcaggcttcgcaactgttttgaaggtcgatcgttgctggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCt                                                                                                                                                                                              <  2:10473‑M2/36‑1 (MQ=255)                                      attgaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCt                                                                                                                                                                                       >  1:234398‑M2/209‑251 (MQ=255)                                         gaccgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCtgtg                                                                                                                                                                                    >  1:175240‑M2/206‑251 (MQ=255)                                            cgtaatgggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCtgtgtg                                                                                                                                                                                  >  1:197642‑M2/203‑250 (MQ=255)                                                  gggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTg                                                                                                                                                                            <  2:189300‑M2/54‑1 (MQ=255)                                                   ggataggtcacgttggtgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGtt                                                                                                                                                                          <  1:151854‑M2/56‑1 (MQ=255)                                                    gataggtcacgttggcgtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGtt                                                                                                                                                                          >  1:37601‑M2/195‑250 (MQ=255)                                                                    gtagatgggcgcatcgtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                                                                                                                         <  2:51919‑M2/73‑1 (MQ=255)                                                                                   gtaaccgtgcatctgccagtttgaggggacgacgacagtatcggcctcaggaagatcgccctccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGGGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                                                                                                                          <  2:118782‑M2/88‑1 (MQ=255)                                                                                                                cgacgacagtatcggcctcaggaagatcgcactccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGc                                                                                                             >  1:91922‑M2/135‑251 (MQ=255)                                                                                                                                       agatcgccctccagccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCTAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACAt                                                                                      <  1:212891‑M2/140‑1 (MQ=255)                                                                                                                                         atcgcactccagccagctttccggcaccgcttctggtaccggaaaccaggcaaagcgccattcgccattcaggctgcgccactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATACAGAGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACAtt                                                                                     >  2:121114‑M2/110‑250 (MQ=255)                                                                                                                                              cctccagccatctttccggaaccgcttcgggggccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGTCGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTGCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAAttgc                                                                               <  2:124477‑M2/147‑1 (MQ=255)                                                                                                                                                     ccagctttccggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTgcgc                                                                        >  1:153422‑M2/98‑251 (MQ=255)                                                                                                                                                               ggcaccgcttctggtgccggaaaccaggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCg                                                              <  2:90808‑M2/164‑1 (MQ=255)                                                                                                                                                                                       caggcaaagcgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGt                                      <  1:109145‑M2/188‑1 (MQ=255)                                                                                                                                                                                                cgccattcgccattcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGATTTCCAGTCGGGAAACCTGTCGTGCCAGCTGc                             >  1:224359‑M2/55‑251 (MQ=255)                                                                                                                                                                                                           attcaggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATc                  >  2:198867‑M2/44‑251 (MQ=255)                                                                                                                                                                                                           attcaggctgcgcaactgttgggaagggcgatcggcgcgggccACGACGGCCAGTAAATCCGTCCTCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATCAAGTGCAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCCCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATc                  >  2:235804‑M2/44‑251 (MQ=255)                                                                                                                                                                                                              caggctgcgcaactgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCTTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGc               <  1:88007‑M2/211‑1 (MQ=255)                                                                                                                                                                                                                           tgttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCgggg   >  1:227654‑M2/28‑250 (MQ=255)                                                                                                                                                                                                                            gttgggcagggcgaccggtgcgggccACGCCGGCCAGCGAATCCGTAATCATGGTCATAGATGTTTACTGTGTGAAATTGTTATCTGCTCACAATTCCACACAACACACGAGCCGGACGCATAAAGTGTAAAGCCTGGGGAGCCTAATGAGTGAGCTAACTCACATTACTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTCCATTAATGAATCGGCCAACGCGCgggg   >  2:76498‑M2/27‑249 (MQ=255)                                                                                                                                                                                                                            gttgggaagggcgatcggtgcgggccACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGga  <  2:91149‑M2/224‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                TCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGA  >  NC_000913/366029‑366497

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: ‑