Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,334,519	T→C	I445T (ATC→ACC)	melB →	melibiose:sodium transporter MelB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,334,519	0	T	C	100.0%	68.7 / NA	20	I445T (ATC→ACC)	melB	melibiose:sodium transporter MelB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/12);  total (8/12)

TGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCC  >  NZ_CP009273/4334444‑4334604                                                                            |                                                                                      tGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATa                                                                         >  2:154821/1‑90 (MQ=255)   tttttCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATaaa                                                                       >  2:212345/1‑90 (MQ=255)    ttttCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAat                                                                      <  2:260576/90‑1 (MQ=255)          tGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCa                                                                >  1:264741/1‑90 (MQ=255)           ggTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCaa                                                               <  1:139500/90‑1 (MQ=255)           ggTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCaa                                                               <  1:212345/90‑1 (MQ=255)                cGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTAc                                                          <  2:110393/90‑1 (MQ=255)                              ttCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGtt                                            >  2:454191/1‑90 (MQ=255)                               tCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTc                                           >  1:158000/1‑90 (MQ=255)                                 cGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTCCCGCCCGAGCCTGTTCAt                                         <  2:325691/90‑1 (MQ=255)                                 cGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCAt                                         <  1:440158/90‑1 (MQ=255)                                  gTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATg                                        >  1:483190/1‑90 (MQ=255)                                  gTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATg                                        >  1:119259/1‑90 (MQ=255)                                  gTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATg                                        <  2:264741/90‑1 (MQ=255)                                             ctaaaTGGTGACCCGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTCCCGCCCGAGCCTGTTCATGCTGATATTCcg                             <  2:319930/87‑1 (MQ=255)                                             cTCAATGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCcg                             <  1:207513/90‑1 (MQ=255)                                                  tGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCcggtcgg                        <  1:326500/90‑1 (MQ=255)                                                  tGGTGACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCcggtcgg                        <  1:238746/90‑1 (MQ=255)                                                      gACACGCTGCGCAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCa                    >  2:375754/1‑90 (MQ=255)                                                              gcgcAGGATCCAGACCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAgtgagtga            <  2:224675/90‑1 (MQ=255)                                                                        cagacccaTCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTc                               >  2:467954/1‑61 (MQ=255)                                                                        cagacccaTCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGcc  <  2:132911/90‑1 (MQ=255)                                                                        cagacccaTCTCCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTc                               <  1:467954/61‑1 (MQ=255)                                                                            |                                                                                      TGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCC  >  NZ_CP009273/4334444‑4334604

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