Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,970,038 T→G 11.4% T84P (ACT→CCT)  fliC ← flagellin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,970,0380TG11.4% 105.2 / 8.0 44T84P (ACT→CCT) fliCflagellin
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base G (3/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

GTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAG  >  NC_002947/4969899‑4970182
                                                                                                                                           |                                                                                                                                                 
gtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgcg                                                                                                                                         <  1:1103139/150‑1 (MQ=255)
gtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgcg                                                                                                                                         >  2:769428/1‑150 (MQ=255)
   gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCt                                                                                                                                      >  2:142367/1‑150 (MQ=255)
      ttAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAACGCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCACCTTca                                                                                                                                   <  1:1050522/150‑1 (MQ=255)
         gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATCGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTAACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcagta                                                                                                                                >  2:573504/1‑148 (MQ=255)
             tACGGGTGAGCTCGCCCGACAGCACGATTACCTCCTGCTGCCTGGAGGGGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCTATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGttt                                                                                                                           <  1:76/150‑1 (MQ=255)
                gggTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTgcg                                                                                                                        <  1:378245/150‑1 (MQ=255)
                  gTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCa                                                                                                                      >  1:1129364/1‑150 (MQ=255)
                         cGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAgaga                                                                                                               <  1:1107549/150‑1 (MQ=255)
                             cGACAGCGCGGTAAACTCCTGCTGCATCGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATAcc                                                                                                           >  1:476541/1‑150 (MQ=255)
                                      ggTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGc                                                                                                  >  1:404907/1‑150 (MQ=255)
                                      ggTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTACCCCAGCCGATCGTTGATTGCAGGGCCAGTTCACGCATACGCTGCATGATATTGCCTGGAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGc                                                                                                  >  1:570668/1‑150 (MQ=255)
                                               ctgctgCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAc                                                                                         >  2:693218/1‑150 (MQ=255)
                                                      aTGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACg                                                                                  >  1:421399/1‑150 (MQ=255)
                                                       tGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGAGAATAATTAGATAAGCACAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGc                                                                                 >  2:29357/1‑150 (MQ=255)
                                                         gAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTc                                                                               <  1:693218/150‑1 (MQ=255)
                                                               cgGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCt                                                                         <  2:476541/150‑1 (MQ=255)
                                                                 gtcgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCttt                                                                       <  2:1031298/150‑1 (MQ=255)
                                                                   cgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                     <  1:67318/150‑1 (MQ=255)
                                                                          gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGt                                                              <  1:572737/150‑1 (MQ=255)
                                                                             gTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                           <  2:637297/150‑1 (MQ=255)
                                                                                         cGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                               <  2:1129364/150‑1 (MQ=255)
                                                                                          gCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcag                                                                                                                                  <  1:840193/67‑1 (MQ=255)
                                                                                          gCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcag                                                                                                                                  >  2:840193/1‑67 (MQ=255)
                                                                                           ccGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGgaga                                             >  2:207302/1‑150 (MQ=255)
                                                                                                  cAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAg                                      <  2:421399/150‑1 (MQ=255)
                                                                                                      gCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttgatttga                                                                >  2:312412/1‑120 (MQ=255)
                                                                                                      gCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttgatttga                                                                <  1:312412/120‑1 (MQ=255)
                                                                                                           ttCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTgcggc                             >  1:675448/1‑150 (MQ=255)
                                                                                                              aCGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGTCACTGGGAGATCTGCAGGCCTGCGGcgtc                          <  1:29357/150‑1 (MQ=255)
                                                                                                                 cATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCgtcgtc                       <  2:404907/150‑1 (MQ=255)
                                                                                                                  aTACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAAc                                                                                   <  2:184002/89‑1 (MQ=255)
                                                                                                                  aTACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAAc                                                                                   >  1:184002/1‑89 (MQ=255)
                                                                                                                                gTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAt        >  2:631857/1‑150 (MQ=255)
                                                                                                                                gTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAt        >  2:1127706/1‑150 (MQ=255)
                                                                                                                                     cTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCa   >  1:310602/1‑150 (MQ=255)
                                                                                                                                     cTGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                           <  2:307062/94‑1 (MQ=255)
                                                                                                                                     cTGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                           >  1:307062/1‑94 (MQ=255)
                                                                                                                                      tGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                               >  2:204389/1‑105 (MQ=255)
                                                                                                                                      tGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                               <  1:204389/105‑1 (MQ=255)
                                                                                                                                      tGAAGTCCCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGCTGGCGTTCTTGACCGCAACGCTCAGACCTTGGATTTGAGTGGTCAGACGGTTAGAGATCTGCAGGCCTGCGGCGTCGTCTTTGTCACTGTTGATTTTCAg  <  2:570668/150‑1 (MQ=255)
                                                                                                                                      tGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg  >  2:872675/1‑150 (MQ=255)
                                                                                                                                       gAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                          >  1:850696/1‑125 (MQ=255)
                                                                                                                                       gAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                          <  2:850696/125‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAG  >  NC_002947/4969899‑4970182

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