Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,408 A→G 20.0% intergenic (+774/+24) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,4080AG20.0% 96.8 / 24.0 55intergenic (+774/+24)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (34/10);  new base G (9/2);  total (43/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

TCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACA  >  NC_002947/4950284‑4950529
                                                                                                                            |                                                                                                                         
tCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTa                                                                                                                >  2:32690/1‑136 (MQ=255)
          gcCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACa                                                                                                      >  1:138255/1‑136 (MQ=255)
                 aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAAc                                                                                               >  1:86239/1‑136 (MQ=255)
                 aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAAc                                                                                               >  1:418413/1‑136 (MQ=255)
                 aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAAc                                                                                               >  2:317439/1‑136 (MQ=255)
                               cgcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCCCGCGGGCGGCGCGGGATCTGACAGGCACCCCCAAACCAACGCCGACCACAGCAAACGAGGCGGCGCCgg                                                                                 >  1:337184/1‑136 (MQ=12)
                                               gaataaCTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCa                                                                 >  2:80904/6‑136 (MQ=255)
                                                 ataaCTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgg                                                               >  2:303633/4‑136 (MQ=255)
                                                     cTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCAATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGgtcgt                                                           >  1:197909/1‑136 (MQ=255)
                                                          ccGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCtta                                                      >  1:175486/1‑135 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:157007/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCCGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:410382/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:85202/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:36527/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:280561/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:77684/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:21488/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCAATAGCCCAACAGATCAACAGTAGCGGCGCAGGTCCGGCGGGATAGGCAGGTCGTccgtaag                                                    >  2:356746/1‑131 (MQ=17)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGCTCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:220076/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTcc                                                                              <  2:168605/110‑1 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTcc                                                                              >  1:168605/1‑110 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:87626/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:24977/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:24480/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:73954/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  2:28055/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:85881/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:172331/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTACTTCat                                                    >  2:204397/1‑135 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCCGGATCGGCAGGTCGTCCTTCaa                                                    >  2:107/1‑135 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCAGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:259399/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACACGCGACAAAAACCAATCGCCGAACAGAACAACTGTAGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAg                                                    >  1:172947/1‑136 (MQ=21)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGACAAAAACCAATACCCCAACAGATCAACAGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTAGTACGTCaa                                                    >  1:246368/1‑135 (MQ=17)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCGCGATCACACAGGCAACGCAAACCTATCGCCGAACACATCAACTGTCGCGGCGCCGGTGCCGCGGGATCCGCAGGTCCTCCATCac                                                    >  1:191297/1‑135 (MQ=21)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCCGCGCTCGATCTCACAGGCAACCCAAACCTATCGCCGAACAGCTAAACTGACGCGGCGCAGGTCCGCCGGGATCGGGAGGTCGTCCTTCAg                                                    >  2:247660/1‑136 (MQ=21)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCCCTCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATAGCCGCACAGATAAACAGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGGGCTTCAg                                                    >  2:211360/1‑136 (MQ=18)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCCCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACAGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTGCTTCAg                                                    >  1:409826/1‑136 (MQ=255)
                                                            gAAACAAATGTAGGAGCGAGCTTGCTCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCCCCGCACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGCTCGGCAGGTCGTCCTTCAg                                                    >  1:385646/1‑136 (MQ=255)
                                                             aaaCAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGt                                                   >  1:43575/1‑136 (MQ=255)
                                                             aaaCAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAAACTATCGCCCAACACATAAACTGTCGCGGCGCACGTCCCGCGGGATCGGCAGGTCGTCCTTAAGt                                                   >  1:299428/1‑136 (MQ=25)
                                                              aaCAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTGGTCCTTCAGTg                                                  >  2:44830/1‑136 (MQ=255)
                                                                     gTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTgg                                           <  1:73954/136‑1 (MQ=255)
                                                                                 ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                               >  2:284587/1‑136 (MQ=255)
                                                                                 ttGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                               >  2:397281/1‑136 (MQ=255)
                                                                                   gcgcgcgATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                             <  2:191297/133‑1 (MQ=255)
                                                                                   gCTCGCGATGCGCCGCTCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                             <  1:211360/136‑1 (MQ=255)
                                                                                   gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                             <  2:21488/136‑1 (MQ=255)
                                                                                   gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                             <  1:216715/136‑1 (MQ=255)
                                                                                   gCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgcccttgc                             <  1:24480/136‑1 (MQ=255)
                                                                                     tCGCGATGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGTATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCg                           <  2:86239/136‑1 (MQ=255)
                                                                                     tCGCGATGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGTATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCg                           <  2:418413/136‑1 (MQ=255)
                                                                                       gcgATGCGCCGCGCGGGCGGCGCTCGATATCAAGGGCAACGCAAACCTATCGCCGAACAGATCAACTGTAGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCATGCCCTTGCCGgc                         >  2:238026/1‑136 (MQ=255)
                                                                                                             cTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTTTTGGCGAATCTGGAac   <  1:115714/136‑1 (MQ=255)
                                                                                                              tCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca  <  1:28055/136‑1 (MQ=255)
                                                                                                              tCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAaca  <  1:121012/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                         
TCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACA  >  NC_002947/4950284‑4950529

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