Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,948,680 Δ1 bp 27.3% coding (213/435 nt) fliO ← flagellar protein FliO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,948,6800C.27.3% 77.7 / 30.1 33coding (213/435 nt)fliOflagellar protein FliO
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base . (5/4);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GGCACTTCCACGGGTTCGGCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCC  >  NC_002947/4948540‑4948815
                                                                                                                                            |                                                                                                                                       
ggCACTTCCACGGGTTCGGCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCa                                                                                                                                <  2:198915/150‑1 (MQ=255)
  cACTTCCACGGGTTCGGCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCAcc                                                                                                                              <  2:261994/150‑1 (MQ=255)
              ttCGGCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGt                                                                                                                  <  1:384689/150‑1 (MQ=255)
                  gCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGctct                                                                                                             >  2:19880/1‑150 (MQ=255)
                      gTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCa                                                                                                          >  1:514856/1‑150 (MQ=255)
                       tACGTGGACGGCCTCGCTGGTGCCGAGCCTATGGGCGATCAGGATCTCCCCCTTGCCAACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAg                                                                                                         <  2:113386/150‑1 (MQ=255)
                               gggCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCgca                                                                                                 <  2:514856/150‑1 (MQ=255)
                                    tCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACGATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAgcc                                                                                           >  2:4707/1‑150 (MQ=255)
                                                           gATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGcc                                                                     >  2:601236/1‑150 (MQ=255)
                                                                         ccTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagca                                                       >  2:542110/1‑150 (MQ=255)
                                                                          cTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagcag                                                      >  1:120846/1‑150 (MQ=255)
                                                                               ccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAgg                                                                       >  2:283024/1‑128 (MQ=255)
                                                                               ccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAgg                                                                       <  1:283024/128‑1 (MQ=255)
                                                                               ccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagcag                                                      <  2:5095/144‑1 (MQ=255)
                                                                               ccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcagcag                                                      >  1:5095/1‑144 (MQ=255)
                                                                                          agcagcagGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTgc                                      >  1:17294/1‑150 (MQ=255)
                                                                                                   cGGTCACGCGGGCTGATGGCACGGCGGCCGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGcc                            <  1:80811/150‑1 (MQ=255)
                                                                                                     gTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGccc                           <  1:542110/150‑1 (MQ=255)
                                                                                                     gTCACGCGGGCCGATGGCACGGCTCGGTGCAAGCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCg                          >  1:323531/1‑150 (MQ=255)
                                                                                                            gggCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCtt                                                                                                                           <  2:72168/47‑1 (MQ=255)
                                                                                                            gggCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCtt                                                                                                                           >  1:72168/1‑47 (MQ=255)
                                                                                                            gggCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGa                    <  2:120846/150‑1 (MQ=255)
                                                                                                                      cACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt          <  2:17294/150‑1 (MQ=255)
                                                                                                                       aCGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTg         >  1:73278/1‑150 (MQ=255)
                                                                                                                       aCGGCTGCCGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa        <  1:19880/150‑1 (MQ=255)
                                                                                                                           cTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCCGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAACATCAGGCCCAGCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgcc     <  2:183068/150‑1 (MQ=255)
                                                                                                                              ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgcc  >  1:516863/1‑150 (MQ=255)
                                                                                                                              ccGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt          <  2:596100/141‑1 (MQ=255)
                                                                                                                              ccGACAATCTCGATACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt          >  1:596100/1‑141 (MQ=255)
                                                                                                                               cGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt          >  2:241745/1‑141 (MQ=255)
                                                                                                                               cGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt          <  1:241745/141‑1 (MQ=255)
                                                                                                                                    aTCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCaa                                                >  2:113813/1‑98 (MQ=255)
                                                                                                                                    aTCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCaa                                                <  1:113813/98‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
GGCACTTCCACGGGTTCGGCCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAGGCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCC  >  NC_002947/4948540‑4948815

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