Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,113,649 T→G K330N (AAA→AAC glcA ← glycolate permease GlcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,113,6490TG100.0% 94.3 / NA 29K330N (AAA→AACglcAglycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/13);  total (16/13)

GCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113563‑3113731
                                                                                      |                                                                                  
gCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTa                                                                                 >  1:188350/1‑90 (MQ=255)
 ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                >  1:443516/1‑90 (MQ=255)
   gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAAc                                                                              <  1:288794/90‑1 (MQ=255)
      ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGc                                                                           <  1:434340/90‑1 (MQ=255)
           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                      <  2:443516/90‑1 (MQ=255)
           cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                      <  1:81942/90‑1 (MQ=255)
                 ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                <  1:1490/90‑1 (MQ=255)
                        aaaTGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTg                                                         >  2:477852/1‑90 (MQ=255)
                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTc                                                                        >  1:494208/1‑64 (MQ=255)
                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTc                                                                        <  2:494208/64‑1 (MQ=255)
                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGt                                              >  2:349024/1‑90 (MQ=255)
                                   cTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGt                                              >  1:172670/1‑90 (MQ=255)
                                     ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                            >  2:398479/1‑90 (MQ=255)
                                        aaTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGa                                         <  2:253203/90‑1 (MQ=255)
                                               cACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaaaaa                                  >  1:24217/1‑90 (MQ=255)
                                                    gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTg                             >  1:119157/1‑90 (MQ=255)
                                                         tAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCAc                        <  1:235070/90‑1 (MQ=255)
                                                              cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                >  1:108828/1‑45 (MQ=255)
                                                              cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                <  2:108828/45‑1 (MQ=255)
                                                              cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCg                   >  1:19721/1‑90 (MQ=255)
                                                                cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                 >  1:478491/1‑90 (MQ=255)
                                                                cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                 <  2:216869/90‑1 (MQ=255)
                                                                cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                 >  2:452456/1‑90 (MQ=255)
                                                                cgccgcCCGGAGCAAACAACGCGTTAAACGGCTTCAt                                                                      >  1:239385/1‑37 (MQ=25)
                                                                cgccgcCCGGAGCAAACAACGCGTTAAACGGCTTCAt                                                                      <  2:239385/37‑1 (MQ=25)
                                                                          aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  2:150954/90‑1 (MQ=255)
                                                                              aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa   >  2:412490/1‑90 (MQ=255)
                                                                              aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa   >  1:113211/1‑90 (MQ=255)
                                                                               aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc  <  2:488400/90‑1 (MQ=255)
                                                                                      |                                                                                  
GCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113563‑3113731

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