Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,947,557 G→C 60.5% P27R (CCG→CGG)  fliQ ← flagellar export apparatus protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,947,5570GC60.5% 24.1 / 46.1 43P27R (CCG→CGG) fliQflagellar export apparatus protein
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base C (13/12);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ACTCCATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGCCCATGATCAGCGCC  >  NC_002947/4947411‑4947703
                                                                                                                                                   |                                                                                                                                                   
aCTCCATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCgg                                                                                                                                                    >  2:167685/1‑148 (MQ=255)
 cTCCATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCAGca                                                                                                                                                  >  3:216974/1‑149 (MQ=255)
    cATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGcacca                                                                                                                                               >  6:239038/1‑149 (MQ=255)
    cATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGcacca                                                                                                                                               >  5:162074/1‑149 (MQ=255)
     ttGAACTTCTGCACCAGCCACGGCCCT‑GAGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGcaccac                                                                                                                                              <  4:216974/148‑1 (MQ=255)
     aTGAACTTCTGCACCAGCCACGGCCCGAGCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTAAGGGTCTGTTCGTTGATATGCGTGCCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGATCCccacca                                                                                                                                               >  1:15261/1‑149 (MQ=255)
                cACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACgg                                                                                                                                      <  2:310497/146‑1 (MQ=255)
                 aCCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGAc                                                                                                                                  >  4:299084/1‑149 (MQ=255)
                      ccACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGa                                                                                                                             <  2:15261/149‑1 (MQ=255)
                         cGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAggg                                                                                                                          <  5:324461/149‑1 (MQ=255)
                          ggCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGggt                                                                                                                         <  3:299084/149‑1 (MQ=255)
                            cccG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGggtgg                                                                                                                       <  2:255971/149‑1 (MQ=255)
                                       atCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAg                                                                                                             >  3:303084/1‑149 (MQ=255)
                                          aGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGc                                                                                                          <  5:239038/149‑1 (MQ=255)
                                            cGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCAt                                                                                                        <  4:303084/149‑1 (MQ=255)
                                                 tCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCa                                                                                                                                >  7:156415/1‑120 (MQ=255)
                                                 tCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCa                                                                                                                                <  8:156415/120‑1 (MQ=255)
                                                        cACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGt                                                                                            >  2:13083/1‑149 (MQ=255)
                                                           cagcagGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGa                                                                                        >  4:200099/1‑149 (MQ=255)
                                                                        aaaaaGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTc                                                                              >  6:60451/1‑146 (MQ=255)
                                                                        aaaaaGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTc                                                                              <  5:60451/146‑1 (MQ=255)
                                                                        aaaaaGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAgg                                                                           >  1:257099/1‑149 (MQ=255)
                                                                                    gTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGTCCGACAGCGACTTCAGGTGTCATcggcg                                                                <  1:167685/149‑1 (MQ=255)
                                                                                    gTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATcggcgg                                                               <  2:257099/149‑1 (MQ=255)
                                                                                     tCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGggtg                                                                                                                        >  6:62927/1‑92 (MQ=255)
                                                                                     tCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGggtg                                                                                                                        <  5:62927/92‑1 (MQ=255)
                                                                                           cGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGc                                                                                    <  3:235211/121‑1 (MQ=255)
                                                                                           cGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGc                                                                                    >  4:235211/1‑121 (MQ=255)
                                                                                           cGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTcc                                                         <  3:200099/148‑1 (MQ=255)
                                                                                               gATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGa                                                     <  1:13083/148‑1 (MQ=255)
                                                                                                aTCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGCGCATCACGCACCAGGT‑CGACAGCGACTTCAGGTGTCATAGGCGGCTCTCCTTGACg                                                   >  8:228358/1‑149 (MQ=255)
                                                                                                     cGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑c                                                                                          <  8:151974/105‑1 (MQ=255)
                                                                                                          ccgccTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGAcg                                         <  6:162074/149‑1 (MQ=255)
                                                                                                              cTGGCACTTCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGGGGTCAGCCACAGTGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACCGCGGCAGAcgccgc                                     <  3:277147/149‑1 (MQ=255)
                                                                                                                  aaCATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCCCCTTGACGGCGTCAGACGCCGCCGa                                  >  3:231753/1‑148 (MQ=255)
                                                                                                                        gCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCGTCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAActgct                           >  1:67680/1‑149 (MQ=255)
                                                                                                                        gCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAActgct                           >  8:209985/1‑149 (MQ=255)
                                                                                                                         ccaccaccaGGCCCACCAACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAActgctg                          >  6:192509/1‑149 (MQ=255)
                                                                                                                            ccaccaGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGCCGCCGCCGAAACTGCTGGc                        >  3:104956/1‑148 (MQ=255)
                                                                                                                               ccaGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCCCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAgg                    <  4:231753/149‑1 (MQ=255)
                                                                                                                                     ccaccaACAGGCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGcc               >  1:255448/1‑148 (MQ=255)
                                                                                                                                              ggCTCCGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGCCCATGATCAgc     >  2:8037/1‑149 (MQ=255)
                                                                                                                                                cTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGTACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGCCCAAGATCAgcgc   >  3:193791/1‑149 (MQ=255)
                                                                                                                                                 tCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGCCCATGATCAGCGcc  >  4:122306/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTCCATGAACTTCTGCACCAGCCACGGCCCG‑GCGATGATCAGCGTGATCAGCATCACCAGCAGGCGCGGCAAAAAGCTCAGGGTCTGTTCGTTGATCTGCGTGGCCGCCTGGAACATCGCCACCACCAGGCCCACCAACAGGCTCGGCACCACCAGCACGGCGACCATGAGGGTGGTCAGCCACAGGGCATCACGGAACAGGT‑CGACAGCGACTTCAGGTGTCATCGGCGGCTCTCCTTGACGGCGTCAGACGCCGCCGAAACTGCTGGCCAGGGTGCCCATGATCAGCGCC  >  NC_002947/4947411‑4947703

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