Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,084,653 A→G G358G (GGT→GGC frvB ← PTS fructose‑like transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,084,6530AG90.0% 46.6 / ‑3.5 21G358G (GGT→GGCfrvBPTS fructose‑like transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/7);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084569‑4084738
                                                                                    |                                                                                     
gACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCccgc                                                                                  >  2:700122/1‑90 (MQ=255)
gACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCccgc                                                                                  >  1:265632/1‑90 (MQ=255)
    gtgGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCcgcg                                                                              >  1:231552/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                      >  1:47528/1‑90 (MQ=255)
            gACCACAATGCCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGCCCCGCCGCCACGggcgg                                                                      >  2:519691/1‑90 (MQ=255)
                      gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  2:434438/1‑90 (MQ=255)
                           actGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAAc                                                       >  1:137234/1‑90 (MQ=255)
                                 aGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATAAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgac                                                 <  1:508596/90‑1 (MQ=255)
                                    ctctTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                              <  2:244275/90‑1 (MQ=255)
                                           tCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAc                                       <  1:39056/90‑1 (MQ=255)
                                             cTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     <  2:209262/90‑1 (MQ=255)
                                                       aaCAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCa                                                                    >  1:47899/1‑49 (MQ=255)
                                                       aaCAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCa                                                                    <  2:47899/49‑1 (MQ=255)
                                                         cAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAAtgtg                         <  2:601193/90‑1 (MQ=255)
                                                           gTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtt                       <  1:208474/90‑1 (MQ=255)
                                                             ttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtttt                     >  1:712406/1‑90 (MQ=255)
                                                             ttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGGGtttt                     >  1:80846/1‑90 (MQ=255)
                                                                     aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             >  1:250607/1‑90 (MQ=255)
                                                                          aCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTgaga        >  1:20738/1‑90 (MQ=255)
                                                                              cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCa    >  2:136822/1‑90 (MQ=255)
                                                                                cAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCaga  <  2:137234/90‑1 (MQ=255)
                                                                                    |                                                                                     
GACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084569‑4084738

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