Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,951,360 C→T 7.4% G219D (GGC→GAC)  fliG ← Flagellar motor switch protein FliG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,951,3600CT 7.4% 178.0 / 4.1 54G219D (GGC→GAC) fliGFlagellar motor switch protein FliG
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base T (2/2);  total (26/28)
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.97e-01

GAGGACACTTCGCGCAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCTTCAGCGCCGCCGGCT  >  ACB143/4951228‑4951495
                                                                                                                                    |                                                                                                                                       
gaggaCACTTCGCGCAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTATTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCAATAAGTGAGCCt                                                                                                                                       >  4:273896/1‑135 (MQ=255)
gaggaCACTGCGCGCAACAGCGCCTGGATGGCCCGGTCGTCGACGGCGGCCAGGTGGTTGAAGACAAACTTGAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGGGCGCCtt                                                                                                                                      <  6:117273/136‑1 (MQ=255)
         tCGCGCAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGTTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCg                                                                                                                                <  7:225013/133‑1 (MQ=255)
         tCGCGCAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGTTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCg                                                                                                                                >  8:225013/1‑133 (MQ=255)
           gcgcAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTg                                                                                                                           >  1:179965/1‑136 (MQ=255)
           gcgcAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCATTCCACCGAGcag                                                                                                                           >  3:102663/1‑134 (MQ=255)
             gcAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTc                                                                                                                         >  8:119618/1‑136 (MQ=255)
              cAACAGCGCCTGGATGCCCCGTTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCg                                                                                                                        <  7:225040/136‑1 (MQ=255)
               aaCAGCGCCTGGATGCCCCGGGCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGa                                                                                                                       <  3:15474/136‑1 (MQ=255)
                aCAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGTCTTCCACCGAGCTGTCGAg                                                                                                                      >  4:207336/1‑136 (MQ=255)
                aCAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGTCTTCCACCGAGCTGTCGAg                                                                                                                      >  4:207333/1‑136 (MQ=255)
                    cgcCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAg                                                                                                                  <  7:144466/136‑1 (MQ=255)
                     gcCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGt                                                                                                                 >  8:95284/1‑136 (MQ=255)
                       cTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTAACGGATCGCGTCCATCAGTGCGCCTTCCACCCAGCTGTCGAGGAAGtt                                                                                                                >  4:289094/1‑135 (MQ=255)
                       cTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCACGTCTTCGATCTGCTCCCACAGGGCGCTGTGGATTTGACCGATAGCGGGCAACAGGGCGCCATACCCCGCCCGGTCGAGGAAggt                                                                                                                >  6:178021/1‑133 (MQ=255)
                         ggATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCat                                                                                                             >  4:278782/1‑136 (MQ=255)
                                  ggTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGTTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGgc                                                                                                    <  5:91557/136‑1 (MQ=255)
                                  ggTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGgc                                                                                                    <  7:9147/136‑1 (MQ=255)
                                   gtcgtcGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGgca                                                                                                   <  3:76772/136‑1 (MQ=255)
                                       tcgACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCAc                                                                                               <  4:219209/136‑1 (MQ=255)
                                       tcgACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCAc                                                                                               <  6:138603/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACTGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  2:94274/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  1:23966/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  3:265107/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  5:218781/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGGCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  3:221339/136‑1 (MQ=255)
                                            gTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGGCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  8:135996/136‑1 (MQ=255)
                                                gCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATac                                                                                      >  5:276828/1‑136 (MQ=255)
                                                 ccAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATacc                                                                                     <  2:179965/136‑1 (MQ=255)
                                                              gACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGgtggtg                                                                        >  8:71811/1‑136 (MQ=255)
                                                               aCAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTgc                                                                       <  3:207066/136‑1 (MQ=255)
                                                               aCAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTgc                                                                       <  6:192442/136‑1 (MQ=255)
                                                                    cATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGcgcgcg                                                                  >  3:199307/1‑136 (MQ=255)
                                                                    cATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATAAAACCCAGTGTGGTGcgcgcg                                                                  <  7:80311/136‑1 (MQ=255)
                                                                                  tCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGcc                                                    >  4:123352/1‑136 (MQ=255)
                                                                                  tCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGcc                                                    >  4:254736/1‑136 (MQ=255)
                                                                                  tCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGcc                                                    >  8:200340/1‑136 (MQ=255)
                                                                                      cTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTgg                                                           <  8:10683/125‑1 (MQ=255)
                                                                                       tCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGCGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATTCCACCCAGGGTGGTGCGCGCGGCGGTGGAGTTGCCCgag                                                >  7:238238/1‑135 (MQ=255)
                                                                                        ccGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGaa                                              >  8:127212/1‑136 (MQ=255)
                                                                                         cGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAAc                                             >  4:70235/1‑136 (MQ=255)
                                                                                              ggTCGCTGTCGATTTCACGGATCGCTTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGGTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCtt                                        >  2:30870/1‑136 (MQ=255)
                                                                                              ggTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCtt                                        >  8:249292/1‑136 (MQ=255)
                                                                                                   cTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTtctc                                     >  8:27373/1‑134 (MQ=255)
                                                                                                      tCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAGCTGCTTCTCGAgg                                 <  7:25103/135‑1 (MQ=255)
                                                                                                      tCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAgg                                 <  1:153264/135‑1 (MQ=255)
                                                                                                             cACGGATCGCGTCCATCAGTGCGTCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGt                          <  3:207336/135‑1 (MQ=255)
                                                                                                             cACGGATCGCGTCCATCAGTGCGTCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGt                          <  3:207333/135‑1 (MQ=255)
                                                                                                              aCGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          <  6:247915/70‑1 (MQ=255)
                                                                                                              aCGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTg                                                                                          >  5:247915/1‑70 (MQ=255)
                                                                                                              aCGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTc                        <  6:185444/136‑1 (MQ=255)
                                                                                                                 gATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGc                     >  4:84866/1‑136 (MQ=255)
                                                                                                                     gcgTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCt                 <  8:144341/136‑1 (MQ=255)
                                                                                                                       gTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCtt                >  3:31796/1‑135 (MQ=255)
                                                                                                                                   gcCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCTTCAGCGCCGCCGGc   >  1:158762/1‑136 (MQ=255)
                                                                                                                                    ccTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCTTCAGCGCCGCCGGCt  <  5:11587/136‑1 (MQ=255)
                                                                                                                                    ccTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCTTCAGCGCCGCCGGCt  <  2:161126/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GAGGACACTTCGCGCAACAGCGCCTGGATGCCCCGGTCGTCGACGTCGGCCAGGTTGTTGAAGACAAACATCAGGTCTTCGATCTGCTCCGACAGGTCGCTGTCGATTTCACGGATCGCGTCCATCAGTGCGCCTTCCACCGAGCTGTCGAGGAAGTTCATGATGTCGGCAGCACGTTTGATACCACCCAGGGTGGTGCGCGCGGCGTTGGAGTTGCCCGAGAACTGCTTCTCGAGGATCTGGTTCAGCTCCTTCAGCGCCGCCGGCT  >  ACB143/4951228‑4951495

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: ‑