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,6530AG97.8% 125.6 / ‑3.5 46G358G (GGT→GGCfrvBPTS fructose‑like transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/24);  minor base C (1/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084575‑4084738
                                                                              |                                                                                     
ggCCCCGACCACAATGGCACTATTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                            <  1:846849/90‑1 (MQ=255)
ggCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                            <  1:416765/90‑1 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                      >  1:926968/1‑90 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                      >  2:1005420/1‑90 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGCCCCGCCGCGACGggcgg                                                                      >  2:1514051/1‑90 (MQ=255)
          acaATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATc                                                                  <  1:531762/90‑1 (MQ=255)
          acaATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATc                                                                  <  1:701457/90‑1 (MQ=255)
             aTGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAGTCGCAAGGCCCGCCGCGACGGGCGGCATCAgcg                                                              >  1:855380/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  1:77924/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  1:1052603/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  2:1250871/1‑90 (MQ=255)
                     actGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAc                                                       >  1:538089/1‑90 (MQ=255)
                        ggCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGcc                                                    >  2:1255634/1‑90 (MQ=255)
                            gCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgacg                                                <  2:302930/90‑1 (MQ=255)
                             ccTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgacga                                               >  1:657429/1‑90 (MQ=255)
                              ctctTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                              <  1:152908/90‑1 (MQ=255)
                              ctctTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                              <  2:235654/90‑1 (MQ=255)
                                ctTTTTCCTGTGGGGTAAACAGTTTGGGCCCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATgg                                            <  1:1268004/90‑1 (MQ=255)
                                 tttttCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGc                                           >  1:1234072/1‑90 (MQ=255)
                                     tCCGGTCGGTTAAACAGTTTGGGCGCAACGACCGCCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAc                                       <  1:1221284/90‑1 (MQ=255)
                                     cccTGTGGGGTAAACAGTTTGGGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAc                                       <  1:927243/89‑1 (MQ=255)
                                       ccgccGGGTAACAAGTTAGGCCGCAATGATCGGCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCa                                     <  1:810563/86‑1 (MQ=255)
                                       cTGTCGGGTAAACAGTTTGGGGGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     <  1:1437279/90‑1 (MQ=255)
                                        tGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAg                                    >  2:1069378/1‑90 (MQ=255)
                                        tGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAg                                    >  1:1306823/1‑90 (MQ=255)
                                         gTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGt                                   >  2:320742/1‑90 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGAc                                                          <  2:1349382/66‑1 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGAc                                                          >  1:1349382/1‑66 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTc                                  >  1:1220537/1‑90 (MQ=255)
                                                     gTTTGGGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtt                       <  2:885950/90‑1 (MQ=255)
                                                     gTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtt                       >  1:103474/1‑90 (MQ=255)
                                                      tttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGttt                      >  1:554166/1‑90 (MQ=255)
                                                      tttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGttt                      >  1:1052076/1‑90 (MQ=255)
                                                       ttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     >  1:695000/1‑74 (MQ=255)
                                                       ttGGGCGCAAGGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     <  2:695000/74‑1 (MQ=255)
                                                          ggCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCAGTCATAAATGTGTTTTTcc                  <  1:1224763/90‑1 (MQ=255)
                                                               aaTGAAGGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             <  1:543247/90‑1 (MQ=255)
                                                               aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             >  1:116932/1‑90 (MQ=255)
                                                               aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             <  2:80821/90‑1 (MQ=255)
                                                                 tGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTg           <  2:538089/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  1:23827/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  1:1209637/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  2:659180/90‑1 (MQ=255)
                                                                        cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGTGTACCAGTCATAAATGTGTTTTTCcagcag             <  1:1471550/81‑1 (MQ=255)
                                                                        cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGTGTACCAGTCATAAATGTGTTTTTCcagcag             >  2:1471550/1‑81 (MQ=255)
                                                                          cAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCaga  <  2:116932/90‑1 (MQ=255)
                                                                              |                                                                                     
GGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084575‑4084738

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