Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,132,502 G→A 7.1% G217S (GGC→AGC)  flgD → flagellar hook assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,132,5020GA 7.1% 175.4 / 2.7 56G217S (GGC→AGC) flgDflagellar hook assembly protein
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base A (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGCAACAATATCGCCA  >  NC_000913/1132361‑1132649
                                                                                                                                             |                                                                                                                                                   
gATGGCCGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACt                                                                                                                                                >  4:154080/1‑147 (MQ=255)
  tGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACCCACCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGGAATACGCTGGCTCTCGGCCCTTAc                                                                                                                                             >  8:331142/1‑148 (MQ=255)
    gcagcaCTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACg                                                                                                                                            <  1:355018/147‑1 (MQ=255)
            gcgcCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGcaccacca                                                                                                                                   >  5:16730/1‑148 (MQ=255)
            gcgcCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGcaccacca                                                                                                                                   >  2:183480/1‑148 (MQ=255)
                cGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCt                                                                                                                               >  4:216868/1‑148 (MQ=255)
                   aCGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTcga                                                                                                                            <  3:216868/148‑1 (MQ=255)
                         cTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCAGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAg                                                                                                                       >  2:103280/1‑147 (MQ=255)
                            aCAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACtt                                                                                                                                               <  1:124714/120‑1 (MQ=255)
                            aCAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACtt                                                                                                                                               >  2:124714/1‑120 (MQ=255)
                            aCAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACg                                                                                                                   >  4:221537/1‑148 (MQ=255)
                                                  tAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACtt                                                                                                                                               >  7:302760/1‑98 (MQ=255)
                                                  tAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACtt                                                                                                                                               <  8:302760/98‑1 (MQ=255)
                                                   aaCGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACAt                                                                                            >  1:348650/1‑148 (MQ=255)
                                                       gtggtACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCa                                                                                        <  4:229726/148‑1 (MQ=255)
                                                         ggtACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAgg                                                                                      <  7:100498/148‑1 (MQ=255)
                                                             cacaACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGAAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGagt                                                                                   >  1:84672/1‑147 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGagtcagtc                                                                              <  3:221537/148‑1 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGagtcagt                                                                               >  1:270831/1‑147 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCAGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCa                                                                                        >  3:179242/1‑138 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCAGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCa                                                                                        >  3:173429/1‑138 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCAGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCa                                                                                        <  4:173429/138‑1 (MQ=255)
                                                                 aCTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCAGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCa                                                                                        <  4:179242/138‑1 (MQ=255)
                                                                   tGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCAt                                                                            >  5:106156/1‑148 (MQ=255)
                                                                      ttGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGc                                                                         >  3:269758/1‑148 (MQ=255)
                                                                      ttGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGc                                                                         >  3:331224/1‑148 (MQ=255)
                                                                          ccAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCttt                                                                     >  3:253927/1‑148 (MQ=255)
                                                                            aGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGCCGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTtc                                                                   >  1:243682/1‑148 (MQ=255)
                                                                            aGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTtc                                                                   >  5:147440/1‑148 (MQ=255)
                                                                              ccGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTtctc                                                                 <  3:344981/148‑1 (MQ=255)
                                                                              ccGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGAAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTtctc                                                                 <  2:84672/148‑1 (MQ=255)
                                                                                        ttGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAg                                                       >  2:284504/1‑148 (MQ=255)
                                                                                        ttGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAg                                                       >  8:202643/1‑148 (MQ=255)
                                                                                           ctctGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCgg                                                    <  1:112632/148‑1 (MQ=255)
                                                                                             ctGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTa                                                                                                       >  3:335572/1‑95 (MQ=255)
                                                                                             ctGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTa                                                                                                       <  4:335572/95‑1 (MQ=255)
                                                                                                gTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGAt                                                  <  1:183480/145‑1 (MQ=255)
                                                                                                gTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTaa                                               <  6:16730/148‑1 (MQ=255)
                                                                                                 tGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTaaa                                              <  1:103280/148‑1 (MQ=255)
                                                                                                 tGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATATCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTaaa                                              <  6:281260/148‑1 (MQ=255)
                                                                                                    aGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGc                                           <  1:209257/148‑1 (MQ=255)
                                                                                                     gggTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTACACGCt                                          >  8:341241/1‑148 (MQ=255)
                                                                                                      ggTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTg                                         >  2:98685/1‑148 (MQ=255)
                                                                                                           gATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTgccgcc                                    <  1:284504/148‑1 (MQ=255)
                                                                                                             tCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGccac                                  >  4:99340/1‑148 (MQ=255)
                                                                                                                      aCAGCGGTAATACGCTGGATCTCGTCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGa                         >  3:306937/1‑148 (MQ=255)
                                                                                                                       cAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAAGTTAAGCCTTCACATGTCAGGAGTCCGTCCTGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGag                        >  3:208480/1‑147 (MQ=255)
                                                                                                                           ggTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGtt                     <  4:306937/147‑1 (MQ=255)
                                                                                                                           ggTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTa                    <  2:348650/148‑1 (MQ=255)
                                                                                                                            gTAATACGCTGTATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATATAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGAGTAATCGCTGCCGCCACCAACCTCGATGTTAt                   <  6:106156/148‑1 (MQ=255)
                                                                                                                               aTACGCTGGATCTCGGCACTTACGGCACCACCTCCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCATGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTTAACGCTGCCGCCACCAACCTCGATGTTATTgg                <  3:96681/148‑1 (MQ=255)
                                                                                                                                 aCGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGc               >  5:221784/1‑147 (MQ=255)
                                                                                                                                 aCGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGc               >  5:136700/1‑147 (MQ=255)
                                                                                                                                  cGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGa                         >  5:310126/1‑136 (MQ=255)
                                                                                                                                  cGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGa                         <  6:310126/136‑1 (MQ=255)
                                                                                                                                          ctcGGCACTTACGGCACCACCACCCACGACGAAGTACGGCAGATAATTTAAGCCATCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGCAACAATATCg     >  6:262255/1‑148 (MQ=255)
                                                                                                                                             ggCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGCAACAATATCGCCa  <  7:202643/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GATGGCAGCACTGCGCCGAACGGTTCTTACAATGTAGCGATTAGCGCCAGTAACGGTGGTACACAACTGGTTGCCCAGCCGCTGCAGTTTGCTCTGGTGCAGGGTGTGATCCGCGGCAACAGCGGTAATACGCTGGATCTCGGCACTTACGGCACCACCACCCTCGACGAAGTACGGCAGATAATTTAAGCCTTCACATTTCAGGAGTCAGTCATGGCCTTTTCTCAAGCGGTTAGCGGATTAAACGCTGCCGCCACCAACCTCGATGTTATTGGCAACAATATCGCCA  >  NC_000913/1132361‑1132649

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