Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,932,635 C→T 8.9% R138H (CGC→CAC)  PP_4342 ← flagellar biosynthesis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,932,6350CT 8.9% 124.3 / 2.7 45R138H (CGC→CAC) PP_4342flagellar biosynthesis regulator
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/2);  total (24/21)
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.07e-01

GCCCTTCCTGCGGGCTCTGCGCCATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGTGCCCGAGG  >  NC_002947/4932489‑4932779
                                                                                                                                                  |                                                                                                                                                
gcCCTTCCTGCGGGCTCTGCGCCATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGcgcg                                                                                                                                                 >  6:356651/1‑148 (MQ=255)
   cTTCCTGCGGGCTCTGCGCCATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGcgaac                                                                                                                                              >  1:164141/1‑148 (MQ=255)
      ccTGCGGGCTCTGCGCCATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGcgaacgaa                                                                                                                                           <  2:164141/148‑1 (MQ=255)
                      cATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACaga                                                                                                                               >  6:198414/1‑144 (MQ=255)
                      cATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACaga                                                                                                                               <  5:198414/144‑1 (MQ=255)
                                     acGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGcgaa                                                                                                                                               <  1:1695/113‑1 (MQ=255)
                                     acGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGcgaa                                                                                                                                               >  2:1695/1‑113 (MQ=255)
                                      cGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCgg                                                                                                           >  3:81572/1‑148 (MQ=255)
                                       ggAAACGGTTCCTGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGGAGGCGTCGGTAATCGAGGTAGGTTCATCACACACGCCCAGCAGCACGTCCTGGGCGGCGCGCACGCAGCTGACCACAGCGGCACCGATAGCCGCAGCGGTa                                                                                                         >  6:357808/1‑149 (MQ=255)
                                       ggAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTa                                                                                                         <  4:81572/149‑1 (MQ=255)
                                        gAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTAt                                                                                                        <  5:356651/149‑1 (MQ=255)
                                        gAAACGGTTCATGCCATAGTAGCGATTGAGCAGCTTGATCAACGCGTAGGCGGCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACTAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTAt                                                                                                        <  8:273940/149‑1 (MQ=255)
                                                           tCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAg                                                                                     <  7:62917/149‑1 (MQ=255)
                                                           tCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAg                                                                                     >  3:209681/1‑149 (MQ=255)
                                                            cgcgATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGTGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAgg                                                                                    >  5:42772/1‑149 (MQ=255)
                                                            cgcgATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGTGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAgg                                                                                    >  5:48040/1‑149 (MQ=255)
                                                             gcgATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGAt                                                                                                    >  1:26797/1‑132 (MQ=255)
                                                             gcgATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGAt                                                                                                    <  2:26797/132‑1 (MQ=255)
                                                             gcgATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGt                                                                                   >  7:47023/1‑149 (MQ=255)
                                                                 ttGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTc                                                                               <  4:209681/149‑1 (MQ=255)
                                                                     gcagcTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCg                                                                           <  8:44849/149‑1 (MQ=255)
                                                                         cTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACAt                                                                                         <  3:210245/131‑1 (MQ=255)
                                                                         cTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACAt                                                                                         >  4:210245/1‑131 (MQ=255)
                                                                             aTCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGTGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCt                                                                   <  6:48040/149‑1 (MQ=255)
                                                                             aTCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGTGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCt                                                                   <  6:42772/149‑1 (MQ=255)
                                                                             aTCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCt                                                                   >  7:321750/1‑149 (MQ=255)
                                                                                aaCGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGaa                                                                >  8:150620/1‑149 (MQ=255)
                                                                                      tcggCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTg                                                          <  8:321750/147‑1 (MQ=255)
                                                                                            tCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCaa                                                    >  6:240567/1‑149 (MQ=255)
                                                                                                 aaTCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGa                                                        <  7:61951/140‑1 (MQ=255)
                                                                                                 aaTCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGa                                                        >  8:61951/1‑140 (MQ=255)
                                                                                                 aaTCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGAt                                                       >  7:114605/1‑141 (MQ=255)
                                                                                                 aaTCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGAt                                                       <  8:114605/141‑1 (MQ=255)
                                                                                                   tCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCa                                             >  5:297809/1‑149 (MQ=255)
                                                                                                      aggtaggtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTAt                                                                                                        >  4:220683/1‑87 (MQ=255)
                                                                                                        gtaggtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAAc                                                               >  4:324208/1‑126 (MQ=255)
                                                                                                        gtaggtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTg                                        >  1:326022/1‑149 (MQ=255)
                                                                                                         taggtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTgg                                       <  7:150620/149‑1 (MQ=255)
                                                                                                            gtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTg                                                          <  1:265130/127‑1 (MQ=255)
                                                                                                            gtTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTg                                                          >  2:265130/1‑127 (MQ=255)
                                                                                                                 tcaCACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGc                               >  4:263269/1‑149 (MQ=255)
                                                                                                                 tcaCACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGc                               >  3:156647/1‑149 (MQ=255)
                                                                                                                       accaccAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATg                         >  4:208610/1‑149 (MQ=255)
                                                                                                                       accaccAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATg                         >  3:158241/1‑149 (MQ=255)
                                                                                                                         caccaGCAGCACTTCCTGGGCGGCGCAAACGAAACTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCa                       >  8:220127/1‑149 (MQ=255)
                                                                                                                                    cTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTgg            <  3:208610/149‑1 (MQ=255)
                                                                                                                                       ccTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGt          <  6:297809/148‑1 (MQ=255)
                                                                                                                                            gcggcgCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGTGCCCGa    <  5:240567/149‑1 (MQ=255)
                                                                                                                                             cggcgCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGTGCCCGAg   <  2:326022/149‑1 (MQ=255)
                                                                                                                                              ggcgCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGTGCCCGAgg  <  5:357808/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCCCTTCCTGCGGGCTCTGCGCCATGTTGGCCAACACACGGAAACGGTTCATGCCATAGTCGCGATTGAGCAGCTTGATCAACGCGTAGGCGTCGGTAATCGAGGTAGGTTCATCACACACCACCAGCAGCACTTCCTGGGCGGCGCGAACGAAGCTGACCACAGAGTCACCGATACCCGCAGCGGTATCGATCACCAACACATCGAGGTTGTCGCCGATCTCGCTGAACGCCTGGATCAACCCGGCATGCTGGGCCGGGGCCAGATGCACCATGCTCTGGGTGCCCGAGG  >  NC_002947/4932489‑4932779

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