Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,057,873	C→G	W418C (TGG→TGC)	yihO ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,057,873	0	C	G	100.0%	138.4 / NA	39	W418C (TGG→TGC)	yihO	MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (23/16);  total (23/16)

TCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCA  >  NZ_CP009273/4057790‑4057952                                                                                    |                                                                                tCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATa                                                                           <  2:443010/90‑1 (MQ=255)      aatCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAAt                                                                      >  2:658780/1‑90 (MQ=255)          aaTGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGc                                                                  >  2:636886/1‑90 (MQ=255)            tGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTg                                                                <  1:620358/90‑1 (MQ=255)               aaaCCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACg                                                             >  1:210733/1‑90 (MQ=255)                aaCCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGc                                                            >  2:749762/1‑90 (MQ=255)                  ccGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCaa                                                          <  2:367533/90‑1 (MQ=255)                    gTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAAc                                                        <  2:628060/90‑1 (MQ=255)                     tttttCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAAcc                                                       >  1:725956/1‑90 (MQ=255)                              gAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGt                                              >  2:558391/1‑90 (MQ=255)                                  gtgtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCa                                          >  1:368019/1‑90 (MQ=255)                                  gtgtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCa                                          >  1:241359/1‑90 (MQ=255)                                  gtgtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCa                                          >  2:130617/1‑90 (MQ=255)                                   tgtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCAt                                         <  2:663399/90‑1 (MQ=255)                                    gtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATc                                        <  1:177019/90‑1 (MQ=255)                                    gtAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGCGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATc                                        >  2:145102/1‑90 (MQ=255)                                          gAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAAt                                                <  2:618914/76‑1 (MQ=255)                                          gAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAAt                                                >  1:618914/1‑76 (MQ=255)                                            aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                >  1:199105/1‑90 (MQ=255)                                            aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                >  2:607231/1‑90 (MQ=255)                                            aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                >  2:550569/1‑90 (MQ=255)                                            aTCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                >  1:638996/1‑90 (MQ=255)                                             tCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcgc                               <  2:241359/90‑1 (MQ=255)                                             tCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcgc                               <  2:641977/90‑1 (MQ=255)                                             tCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcgc                               <  2:638996/90‑1 (MQ=255)                                               ccATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGa                             <  1:648140/90‑1 (MQ=255)                                                  tCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGt                          >  1:557357/1‑90 (MQ=255)                                                  tCGTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGt                          >  1:532826/1‑90 (MQ=255)                                                    gTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTg                        <  1:658780/90‑1 (MQ=255)                                                    gTCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTg                        >  1:469644/1‑90 (MQ=255)                                                     tCAATGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTgg                       <  2:307689/90‑1 (MQ=255)                                                       aaTGCAGCAATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGa                     <  1:282607/90‑1 (MQ=255)                                                             gcaATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                <  2:132465/73‑1 (MQ=255)                                                             gcaATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTgcg                                >  1:132465/1‑73 (MQ=255)                                                             gcaATAATTGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATaa               >  1:607010/1‑90 (MQ=255)                                                             gcaATAATTGCCAGTACACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATaa               >  2:711027/1‑90 (MQ=255)                                                                     tGCCAGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTg       >  2:247753/1‑90 (MQ=255)                                                                         aGTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTc   <  1:130617/90‑1 (MQ=255)                                                                          gTGCACAAGGGCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCa  >  2:329055/1‑90 (MQ=255)                                                                                    |                                                                                TCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAAGAATCCCATCGTCAATGCAGCAATAATTGCCAGTGCACAAGGCCAGATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCA  >  NZ_CP009273/4057790‑4057952

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