Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,333 A→G 12.2% intergenic (+693/+105) 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,945,3330AG12.2% 98.0 / 6.3 41intergenic (+693/+105)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base G (3/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

CCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCG  >  NC_002947/4945198‑4945477
                                                                                                                                       |                                                                                                                                                
ccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGCCGGTCGCGCGCGCCGGGCACAGGAGGGACTgg                                                                                                                                     <  5:378364/149‑1 (MQ=255)
      ggcgCATCGCGAGAAAGCTCGCGCCTACATTTGTTGCGGGCAATTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCCATATTGCGCCGTAAAACCACGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaa                                                                                                                               <  2:90517/149‑1 (MQ=255)
                  gCATGCGCGCTCCTAATTTTGTTTCGGGCCTGTTATTCCGGCGGTTTTTACGCGCGACCGCCTGCGCTTACACCTCCATATGGCGCCTTACAACCAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                                                                                   <  6:83031/149‑1 (MQ=21)
                     aGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagcg                                                                                                                >  3:87411/1‑149 (MQ=255)
                         cgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACACGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTa                                                                                                                      <  6:167101/139‑1 (MQ=255)
                         cgctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTa                                                                                                                      >  5:167101/1‑139 (MQ=255)
                          gctcCTACATTTGTTTAGGGCCAGTTATTCATGCTGGATTTACGCGCGAACGCCTTCGCGTACAACTAGATATTGCGCCGTACAAACAAGGCGGTCGCGAGCGCCGGGCACAGGAGGGACTGGCACGAAACAAATGTAGGAGCGAGCtt                                                                                                           >  8:235073/1‑149 (MQ=25)
                               tACATTTGTTTGGGGCATTTTTTTCGTGCGGTTTTAACGCGCGACCGCCTGCGCTACCCCCTCCATATTGCGCTGTACAAACACGGCGGTGGGGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTcg                                                                                                      <  4:81032/149‑1 (MQ=21)
                                      gTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATgcg                                                                                               >  8:277619/1‑149 (MQ=255)
                                       tttCGGGCCATTAATTCCGCCGGGATTTACGCGCGACCCCCTTCGCTTACACCTCGATATTGCCCCGTACAAACATGGCGGGCGGGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATgcg                                                                                               <  4:87411/148‑1 (MQ=37)
                                         tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCAATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgccg                                                                                            >  1:98027/1‑149 (MQ=255)
                                           gggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGc                                                                                                                                    >  7:41152/1‑107 (MQ=255)
                                           gggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCcgcg                                                                                          >  1:301326/1‑149 (MQ=255)
                                            ggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCcgcgc                                                                                         >  1:172580/1‑149 (MQ=255)
                                                        cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGTAAATCCCGCAGGAATAACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTc                                                                              >  2:257326/1‑148 (MQ=17)
                                                        cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCcgcgc                                                                                         >  1:349022/1‑137 (MQ=255)
                                                        cTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTc                                                                              >  3:374984/1‑148 (MQ=255)
                                                         tGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGa                                                                            >  5:34673/1‑149 (MQ=255)
                                                                  tACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                    >  3:285195/1‑148 (MQ=255)
                                                                  accgcgcgACCCCCTTCGCGTACACCTCGATATTGCGCCTTACAACCATGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                    <  4:285195/146‑1 (MQ=255)
                                                                       gcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACg                                                              >  7:224779/1‑149 (MQ=255)
                                                                            cGCCTTCGCTTACACCTCTATTTTGCGCCGTAAAAACACGGCGGTGGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGACAAAAAc                                                         <  2:98027/149‑1 (MQ=18)
                                                                            cGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAAc                                                         >  1:154341/1‑149 (MQ=255)
                                                                               cTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGTCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTa                                                      >  2:287258/1‑149 (MQ=255)
                                                                                           cTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGaa                                              >  5:391399/1‑145 (MQ=255)
                                                                                           cTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGa                                          <  1:287258/149‑1 (MQ=255)
                                                                                                 atTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTca                                                                       <  2:290864/114‑1 (MQ=255)
                                                                                                 atTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTca                                                                       >  1:290864/1‑114 (MQ=255)
                                                                                                  ttGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTg                                   >  6:7802/1‑149 (MQ=255)
                                                                                                   tGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGt                                  >  3:68911/1‑149 (MQ=255)
                                                                                                           acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaaacctccggcaaacacctcaaaattccaacgccg                         <  5:215920‑M2/149‑83 (MQ=255)
                                                                                                            caaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcggcg                           >  4:91089/1‑147 (MQ=255)
                                                                                                             aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCa                                                            <  5:357056/113‑1 (MQ=255)
                                                                                                             aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCa                                                            >  6:357056/1‑113 (MQ=255)
                                                                                                                     gggggccgcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCcggcg                <  6:34673/143‑1 (MQ=255)
                                                                                                                      ggTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctcc                                                    <  1:262891‑M2/111‑56 (MQ=255)
                                                                                                                      ggTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctcc                                                    >  2:262891‑M2/1‑56 (MQ=255)
                                                                                                                        tCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaaacctccggcaaacacctcaaaattccaacgccgtccactcaaaatt            >  6:215920‑M2/1‑54 (MQ=255)
                                                                                                                            gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCgg         >  5:389233/1‑149 (MQ=255)
                                                                                                                                gcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgg     <  8:224779/149‑1 (MQ=255)
                                                                                                                                   gggCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGgtcg  <  2:154341/149‑1 (MQ=255)
                                                                                                                                       acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTcg                                >  2:53221/1‑115 (MQ=255)
                                                                                                                                       acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTcg                                >  2:58692/1‑115 (MQ=255)
                                                                                                                                       acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTcg                                <  1:58692/115‑1 (MQ=255)
                                                                                                                                       acaGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTcg                                <  1:53221/115‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCG  >  NC_002947/4945198‑4945477

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: