Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,015,409 A→G E177G (GAA→GGA)  fliG → flagellar motor switch protein FliG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,015,4090AG96.5% 98.5 / NA 29E177G (GAA→GGA) fliGflagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/13);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGAAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAATTA  >  NC_000913/2015332‑2015486
                                                                             |                                                                             
aCGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGt                                                                            <  1:353374/81‑1 (MQ=255)
 cGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTa                                                                           <  2:353375/81‑1 (MQ=255)
  gTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTAc                                                                          >  1:353423/1‑81 (MQ=255)
       cgcCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAAt                                                                     >  2:353415/1‑81 (MQ=255)
        gcCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATg                                                                    >  1:353435/1‑81 (MQ=255)
               cgTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCt                                                             >  1:353437/1‑81 (MQ=255)
                       tGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGt                                                      >  2:353434/1‑80 (MQ=255)
                         cGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAg                                                   <  2:353376/81‑1 (MQ=255)
                          gTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGa                                                  <  2:353378/81‑1 (MQ=255)
                          gTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGa                                                  <  2:353377/81‑1 (MQ=255)
                            aTCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAt                                                >  1:353421/1‑81 (MQ=255)
                               gCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAAtctc                                             <  1:353379/81‑1 (MQ=255)
                                  aCCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAg                                          >  2:353419/1‑81 (MQ=255)
                                       tGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCgcag                                     >  1:353420/1‑81 (MQ=255)
                                           ggcgTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCaaa                                 <  1:353380/81‑1 (MQ=255)
                                           ggcgTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCAACGGTCAGAATCTCAAGCGCAGCaa                                  >  1:353427/1‑80 (MQ=255)
                                              gTGCAGCCAGCCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATg                              >  2:353418/1‑81 (MQ=255)
                                                 cagccagcCGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATggg                            <  1:353381/80‑1 (MQ=255)
                                                 cagccagcCGCGCTGGCGGAGCTGACCGAAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGggc                           <  1:353382/81‑1 (MQ=255)
                                                  agccagccGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGggcg                          >  2:353422/1‑81 (MQ=255)
                                                    ccagccGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGggcggc                        <  1:353383/81‑1 (MQ=255)
                                                     cagccGCGCTGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGgcggcg                       >  2:353433/1‑81 (MQ=255)
                                                              tGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACgg              >  1:353436/1‑81 (MQ=255)
                                                              tGGCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACgg              >  2:353430/1‑81 (MQ=255)
                                                                gCGGAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCa            <  1:353384/81‑1 (MQ=255)
                                                                  ggAGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGc          >  1:353431/1‑81 (MQ=255)
                                                                    aGCTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCg        <  2:353385/81‑1 (MQ=255)
                                                                      cTGACCGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGaa      <  2:353386/81‑1 (MQ=255)
                                                                          ccGGAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAATTa  <  1:353387/81‑1 (MQ=255)
                                                                             |                                                                             
ACGTCTGCGCCACGACGTGATGTTGCGTATCGCCACCTTTGGCGGCGTGCAGCCAGCCGCGCTGGCGGAGCTGACCGAAGTACTGAATGGCTTGCTCGACGGTCAGAATCTCAAGCGCAGCAAAATGGGCGGCGTGAGAACGGCAGCCGAAATTA  >  NC_000913/2015332‑2015486

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