Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,242,942	Δ66 bp	coding (510‑575/891 nt)	malG ←	maltose transporter subunit

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4242942	4243007	66	49 [0]	[1] 49	malG	maltose transporter subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4242941	0 (0.000)	44 (0.740)	38/184	0.8	98.9%	coding (576/891 nt)	malG	maltose transporter subunit
	?	NC_000913	4243008 =	1 (0.020)					coding (509/891 nt)	malG	maltose transporter subunit

AGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCG  >  NC_000913/4242842‑4243037                                                                                                    |                                                                                                 agggataAACGACAGGAGGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACATCACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCA                                                                                                  <  2:292287‑M1/97‑1 (MQ=255) aGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCA                                                                                                  <  1:932768‑M1/100‑1 (MQ=255)      tAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatccc                                                                                             <  1:1460896‑M1/100‑6 (MQ=255)       aaaCGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatccca                                                                                            >  1:97200‑M1/1‑94 (MQ=255)       aaaCGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatccca                                                                                            >  1:97205‑M1/1‑94 (MQ=255)        aaCGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGCCAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccac                                                                                           >  2:1034392‑M1/1‑93 (MQ=255)          cGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccaccc                                                                                         <  2:1235403‑M1/100‑10 (MQ=255)               gaatgAATACAACCGCCAGAATCGGTACTGACAGCGGCAACAGAACAAGCCGGAAGACCTGCCACGGGGTCGAACCATCCAGCGCAatcccacccagata                                                                                    <  2:213588‑M1/98‑15 (MQ=255)                 aTGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacg                                                                                  >  2:1563889‑M1/1‑84 (MQ=255)                      tACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaa                                                                             >  1:117298‑M1/1‑79 (MQ=255)                       accaccGCCAGAATCCGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaa                                                                             <  2:920293‑M1/99‑22 (MQ=255)                         caccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcgcacccagatacgcgaaaatt                                                                          >  2:562763‑M1/1‑76 (MQ=255)                         caccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaatt                                                                          >  1:767128‑M1/1‑76 (MQ=255)                           ccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattac                                                                        >  1:1538309‑M1/1‑74 (MQ=255)                           ccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattac                                                                        <  1:804074‑M1/100‑27 (MQ=255)                            cGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacg                                                                       <  2:1429297‑M1/100‑28 (MQ=255)                                    tCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtg                                                               <  1:477662‑M1/100‑36 (MQ=255)                                     cGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtga                                                              >  2:187514‑M1/1‑64 (MQ=255)                                          cTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtatt                                                         >  1:1504980‑M1/1‑59 (MQ=255)                                           tGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtcttc                                                        >  2:59147‑M1/1‑58 (MQ=255)                                           tGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattc                                                        <  1:508269‑M1/100‑43 (MQ=255)                                            gACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattca                                                       <  1:1225615‑M1/100‑44 (MQ=255)                                             aCAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcag                                                      >  2:362053‑M1/1‑56 (MQ=255)                                             aCAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcag                                                      <  2:98115‑M1/100‑45 (MQ=255)                                              cAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcagg                                                     <  1:1439640‑M1/100‑46 (MQ=255)                                                    cAACAGGACAAGGCGGAAGGCCTTCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatg                                               <  2:1191574‑M1/100‑52 (MQ=255)                                                     aaCAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatga                                              >  1:108430‑M1/1‑48 (MQ=255)                                                       cAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgcgtattcaggccaatga                                              >  2:1450‑M1/1‑46 (MQ=255)                                                        aGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatg                                           <  2:978059‑M1/100‑56 (MQ=255)                                                            cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaat                                       >  1:373121‑M1/1‑41 (MQ=255)                                                            cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaat                                       >  1:373078‑M1/1‑41 (MQ=255)                                                             aaGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatg                                      <  2:1452050‑M1/100‑61 (MQ=255)                                                                   gAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactca                                >  1:1213154‑M1/1‑34 (MQ=255)                                                                       gCCTGCCACGGTGTCGCACCATCCAGCGCGatcccacccatatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcaccca                            <  2:680127‑M1/100‑71 (MQ=255)                                                                        ccTGCCACGGTGTCGCACCATCCAGCGCAatcccacccggatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccag                           >  2:969688‑M1/1‑29 (MQ=255)                                                                        ccTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccag                           >  2:1110412‑M1/1‑29 (MQ=255)                                                                         cTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccaga                          >  1:1050028‑M1/1‑28 (MQ=255)                                                                         cTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccaga                          <  1:250128‑M1/100‑73 (MQ=255)                                                                          tGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagac                         <  1:1055090‑M1/100‑74 (MQ=255)                                                                            ccACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacga                       >  2:591340‑M1/1‑25 (MQ=255)                                                                              aCGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatc                     <  2:1024529‑M1/100‑78 (MQ=255)                                                                              aCGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatc                     <  1:963809‑M1/100‑78 (MQ=255)                                                                                  tgtCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaac                 <  1:927874‑M1/100‑82 (MQ=255)                                                                                      gCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacg             >  1:1431348‑M1/1‑15 (MQ=255)                                                                                      gCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacg             >  1:1216369‑M1/1‑15 (MQ=255)                                                                                         ccATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcat          >  1:1465284‑M1/1‑12 (MQ=255)                                                                                              cAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaac     >  1:949141‑M1/1‑7 (MQ=255)                                                                                                gcgcAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaacgc   >  2:774424‑M1/1‑5 (MQ=255)                                                                                                 cgcAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaacgcg  <  2:708600‑M1/100‑97 (MQ=255)                                                                                                    |                                                                                                 AGCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCG  >  NC_000913/4242842‑4243037

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