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	45 [0]	[0] 45	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)	41 (0.820)	34/184	0.6	100%	coding (576/891 nt)	malG	maltose transporter subunit
	?	NC_000913	4243008 =	0 (0.000)					coding (509/891 nt)	malG	maltose transporter subunit

GCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCGCCG  >  NC_000913/4242843‑4243040                                                                                                   |                                                                                                    gCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAa                                                                                                    >  1:605965‑M1/1‑99 (MQ=255)  cGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAat                                                                                                   >  1:1129284‑M1/1‑98 (MQ=255)    aTAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcc                                                                                                 <  1:730120‑M1/100‑5 (MQ=255)    aTAAACGACAGGACGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcc                                                                                                 >  1:875374‑M1/1‑96 (MQ=255)       aaCGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccac                                                                                              >  2:1073411‑M1/1‑93 (MQ=255)         cGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccaccc                                                                                            <  2:654589‑M1/100‑10 (MQ=255)           aCAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccag                                                                                          <  1:1084543‑M1/100‑12 (MQ=255)             aGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagat                                                                                        <  1:781816‑M1/100‑14 (MQ=255)             aGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagat                                                                                        <  2:605961‑M1/100‑14 (MQ=255)              ggATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagata                                                                                       <  1:526398‑M1/100‑15 (MQ=255)               gATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatac                                                                                      <  1:518584‑M1/100‑16 (MQ=255)                   aaTACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcga                                                                                  >  2:1258255‑M1/1‑81 (MQ=255)                   aaTACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcga                                                                                  >  2:682606‑M1/1‑81 (MQ=255)                     tACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcga                                                                                  >  1:553718‑M1/1‑79 (MQ=255)                       ccaccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaat                                                                              >  1:1014119‑M1/1‑77 (MQ=255)                       ccaccGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaat                                                                              >  2:151277‑M1/1‑77 (MQ=255)                                     ggTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgag                                                                <  2:237876‑M1/100‑38 (MQ=255)                                              aGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggc                                                       <  1:317892‑M1/100‑47 (MQ=255)                                                  gCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaat                                                   >  1:263792‑M1/1‑50 (MQ=255)                                                   cAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatg                                                  <  2:193267‑M1/100‑52 (MQ=255)                                                   cAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatg                                                  <  1:210763‑M1/100‑52 (MQ=255)                                                        ggACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatgg                                             >  1:842214‑M1/1‑44 (MQ=255)                                                           cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaat                                          <  1:614396‑M1/100‑60 (MQ=255)                                                           cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaat                                          >  2:1202263‑M1/1‑41 (MQ=255)                                                           cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaa                                           >  1:1261516‑M1/1‑41 (MQ=255)                                                           cAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgacgtgagtattcaggccaatgaatggaat                                          >  2:41610‑M1/1‑41 (MQ=255)                                                            aaGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatg                                         >  1:1090743‑M1/1‑40 (MQ=255)                                                              ggCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgt                                        <  2:479728‑M1/99‑62 (MQ=255)                                                                cGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtact                                     >  1:1026181‑M1/1‑36 (MQ=255)                                                                 ggAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactc                                    >  1:351056‑M1/1‑35 (MQ=255)                                                                  gAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaaggtactca                                   <  1:685982‑M1/100‑67 (MQ=255)                                                                    aGGCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacc                                 <  1:791399‑M1/100‑69 (MQ=255)                                                                      gCCTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcaccca                               >  2:1249803‑M1/1‑30 (MQ=255)                                                                        cTGCCACGGTGTCGCACCATCCAGTGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccaga                             >  2:738247‑M1/1‑28 (MQ=255)                                                                        cTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccaga                             <  2:446407‑M1/100‑73 (MQ=255)                                                                        cTGCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccaga                             >  1:1217318‑M1/1‑28 (MQ=255)                                                                          gCCACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacg                           >  1:422647‑M1/1‑26 (MQ=255)                                                                            cACGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgat                         >  1:607661‑M1/1‑24 (MQ=255)                                                                             aCGGTGTCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatc                        >  2:1181943‑M1/1‑23 (MQ=255)                                                                                 tgtCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaac                    <  1:1181954‑M1/100‑82 (MQ=255)                                                                                  gtCGCACCATCCAGCGCAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaaca                   <  2:1014110‑M1/100‑83 (MQ=255)                                                                                               gcgcAatcccacccagatacgctaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaacg       <  1:317137‑M1/99‑95 (MQ=255)                                                                                               gcgcAatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaacgc      <  2:722521‑M1/100‑96 (MQ=255)                                                                                                cgcAatcccacccagatacgcgaaaattacgccgccgtgaggattcaggccactgactggaatgtactcacacagacgcacaaacaacgcatacaacgcg     >  2:172543‑M1/1‑4 (MQ=255)                                                                                                   aatcccacccagatacgcgaaaattacgccgccgtgagtattcaggccaatgaatggaatgtactcacccagacgatcaaacaacgcatacaacgcgacc  <  2:570540‑M1/100‑100 (MQ=255)                                                                                                   |                                                                                                    GCGATAAACGACAGGATGAATACCACCGCCAGAATCGGTACTGACAGCGGCAACAGGACAAGGCGGAAGGCCTGCCACGGTGTCGCACCATCCAGCGCAGCAGCTTCTTCCAGCGAACTGTCGATGGTTTCGAAATAGCCTTTGATGGTCCAGACATGCAGCGCAATCCCACCCAGATACGCGAAAATTACGCCGCCG  >  NC_000913/4242843‑4243040

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