Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,964,698 A→T L118H (CTC→CAC)  flhA ← flagellar biosynthesis protein FlhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,6980AT89.3% 72.7 / NA 28L118H (CTC→CAC) flhAflagellar biosynthesis protein FlhA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (14/11);  minor base G (0/3);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTAAACAACAGGTCGAGTATGAATGCGGGCAGTGGCAGC  >  NC_000913/1964487‑1964943
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cACTTCGGAGCTGCGTGTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCGTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGGTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCGTCGACCACCTTCCCTgccgcc                                                                                                                                                                                                                 <  2:308326/250‑1 (MQ=255)
    gcGGAGCTGCGTTTTTTCGCCTCACCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGGGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgcc                                                                                                                                                                                                              <  2:114452/248‑1 (MQ=255)
          gggcgGTTTGTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGGCGGCGTCAATCTCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCTACTTCTGCGATACGCCCGGCCCCTTTGGTAATGACCATAAAGTTGATGATCACGAGACTGACAAACACCACGATACCGATAGCGAAATTGCCTCCAACGGGGAAGTGATCGATCGCGTCGACCACCTTCCCTGCCGCCGCTGCGCCGGt                                                                                                                                                                                                      <  2:311037/250‑1 (MQ=255)
                          cATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGTTATGCCCTTCCATTa                                                                                                                                                                                        <  1:76713/249‑1 (MQ=255)
                               tCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGATTCGACCACCTTCCCTGCCGCCGCCGCGACGGTATGCCCTTCCATTAAAATGAt                                                                                                                                                                                 >  1:66071/1‑251 (MQ=255)
                                                gTTAAGGTCGGCGTCAATCGCCATATGCTTACCCGGCCAACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACATTTGGAAATGACCATAAAGTTGATGATCACGAGAATGACAAACAACAAGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACATTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTGAAATGATACGGGTTGAAGCCACGt                                                                                                                                                                >  2:296532/1‑251 (MQ=255)
                                                        cGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                                                                                                                                                        >  1:176762/1‑251 (MQ=255)
                                                         ggCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACAGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                                                                                                                                                        >  2:80559/1‑250 (MQ=255)
                                                                               cccGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGacc                                                                                                                                                                                                                                >  2:48596/1‑156 (MQ=255)
                                                                               cccGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGacc                                                                                                                                                                                                                                <  1:48596/156‑1 (MQ=255)
                                                                                                       cgcgcACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGAATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACTACAGAATGGTCGGAAACGCAGCaa                                                                                                          >  2:217711/1‑250 (MQ=255)
                                                                                                                    cTGCGAAAAGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGCAGTAAACGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGGTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAACCTCAAGCGTgcac                                                                                            >  2:301324/1‑249 (MQ=255)
                                                                                                                              cccGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAAc                                                                                   >  1:136400/1‑250 (MQ=255)
                                                                                                                                              aTGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGCAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAaca                                                                   >  2:11206/1‑250 (MQ=255)
                                                                                                                                                         gTTGATGATCACGAGAATGACAAACACCACGATACCTATAGCGAAATTGCCACCACCGTGGAAGTGACCGAACGCTTCTTCCCCCTGCCCTGCCGACGCCGCGCCGGTATGCCCGGCCATAAAAATTATACGGGGTGAAGCCTCGTTAAGTGCCAGACGCAACAGCGTTGTAAACAACAGAATGGTCTGATACGCAGCGAACTCAAGCGTGCGCTTTGGAAACAGCGCCATCAGCACCGCCATTATCGACa                                                       <  2:29183/251‑1 (MQ=255)
                                                                                                                                                         gTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCTTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACa                                                       <  1:11206/251‑1 (MQ=255)
                                                                                                                                                               gatCACGAGAATGACAAACACCACGTTACCGCTCGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGAATGCCCATCCATTAAAATGATACGGATTGAAGCCACGTTAAGTACCAGACGCAACAGCGTGTTAACCAACAGAATTGTCGGAAACACAGCAAACTCAAGAGTGCGCTGGGTAAACAGCCACACCATAAACACCAAGATCGAAAAGGCa                                                  >  2:169958/1‑250 (MQ=255)
                                                                                                                                                               gatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCaa                                                 <  1:80559/251‑1 (MQ=255)
                                                                                                                                                                atCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAat                                                >  1:174677/1‑251 (MQ=255)
                                                                                                                                                                      agaATGACAAACACCACGTGACCGATATCTAAATTGCCACCATCGGGTAAGTGGCCGAACGCGTCTACCACCTGCCCTGCCGCCGCCGCTCTGTGATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGTCGCAACAGCGGGGTAAACAACAGAATGGTCGGAAACGCAGCAAGCTCAATCGTGCGCTGGGTAAACAGCGCCACCAGCAACACCATGATCTGCAAGGCAATATTg                                            <  2:285914/249‑1 (MQ=255)
                                                                                                                                                                           gACAAACACCACGATACCGATAGCGAAATTGCCTCCAACGTGGAAGTGACCGAACGCTCCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTTCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTaa                                     <  1:163265/251‑1 (MQ=255)
                                                                                                                                                                           gACAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTa                                      >  1:125475/1‑250 (MQ=255)
                                                                                                                                                                            aCAAACACCACGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTaa                                     >  1:43451/1‑250 (MQ=255)
                                                                                                                                                                             cAAACCCCACGATCCCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGAATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGACTGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGAGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTCaac                                   >  2:156933/1‑251 (MQ=255)
                                                                                                                                                                             cAAACACCACGAGACCGATAGCGAAATTCCCACTAACGTGAAAGTGACCGACCGCTTCGACCACCGTCCCGGCCGCCGCAGCGCCGGTATGCCCTGCCATGAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGGAGTGTCGGAAACGCAGCAAACTCAAGCGTGCGCTTGGTAAACATCACCACCAGCAACACCGTGATCGACAAGACAATATTGATGGTAaac                                   <  2:261863/251‑1 (MQ=255)
                                                                                                                                                                                  cccacGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTAAACAACAg                              <  1:211698/250‑1 (MQ=255)
                                                                                                                                                                                     acGATACCGATAGCGAAATTGCCACCAACGTGGAAGTGACCGAACTCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATGCGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAGCACCATGATCGACAAGGCAATATTGAAGGTAAACAACAGGTc                           <  1:296532/251‑1 (MQ=255)
                                                                                                                                                                                                              caACGTGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTAAACAACAGGTCGAGTATGAATGCGGGCAGTGgcagc  <  1:217711/251‑1 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
CACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAACGCAGCAAACTCAAGCGTGCGCTGGGTAAACATCGCCACCAGCAACACCATGATCGACAAGGCAATATTGAAGGTAAACAACAGGTCGAGTATGAATGCGGGCAGTGGCAGC  >  NC_000913/1964487‑1964943

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