Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,014,548 A→G 9.3% P440P (CCA→CCG)  fliF → flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,014,5480AG 9.3% 119.7 / 3.7 43P440P (CCA→CCG) fliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (2/2);  total (24/19)
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.95e-01

AACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAGCAA  >  NC_000913/2014403‑2014692
                                                                                                                                                 |                                                                                                                                                
aaCCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTc                                                                                                                                               >  6:142451/1‑149 (MQ=255)
        ccTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCaa                                                                                                                                         >  1:43891/1‑147 (MQ=255)
         ctctCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAg                                                                                                                                      >  6:121132/1‑149 (MQ=255)
            tcAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGCACTGCCATTCTGGCAACAGc                                                                                                                                     >  1:40526/1‑147 (MQ=255)
                cAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGt                                                                                                                               <  4:293577/149‑1 (MQ=255)
                cAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGt                                                                                                                               <  2:49844/149‑1 (MQ=255)
                   cGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTa                                                                                                                            <  6:58063/149‑1 (MQ=255)
                    gAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTAt                                                                                                                           <  7:73146/149‑1 (MQ=255)
                          aTGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGAt                                                                                                                       >  7:248278/1‑147 (MQ=255)
                          aTGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCa                                                                                                                     <  1:149351/149‑1 (MQ=255)
                                  aaTTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCGTTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTtg                                                                                                             >  1:213271/1‑149 (MQ=255)
                                  aaTTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCGTTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTtg                                                                                                             >  5:141144/1‑149 (MQ=255)
                                             tGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGc                                                                                                  >  8:178429/1‑149 (MQ=255)
                                               aCCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCGTTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGc                                                                                                  <  2:213271/147‑1 (MQ=255)
                                               aCCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCGTTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGc                                                                                                  <  6:141144/147‑1 (MQ=255)
                                                cccGCGAGGCGCTGGGCTTTTCTGAAATACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGGAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTgg                                                                                               <  4:159590/149‑1 (MQ=255)
                                                 ccGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGt                                                                                              >  3:150432/1‑149 (MQ=255)
                                                   gcgAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTg                                                                                            <  3:132340/149‑1 (MQ=255)
                                                     gAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCt                                                                                          >  8:93361/1‑149 (MQ=255)
                                                      aGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTg                                                                                         <  7:98064/149‑1 (MQ=255)
                                                       ggCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTgg                                                                                        <  7:93361/149‑1 (MQ=255)
                                                          gATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTAc                                                                                     >  7:32302/1‑149 (MQ=255)
                                                           aTGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCTGAGCACTGCCATTCTGGCCACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTAct                                                                                    >  4:233790/1‑149 (MQ=255)
                                                            tGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTActg                                                                                   <  5:142451/149‑1 (MQ=255)
                                                                      tGAATAACGCGGTGACTCGCGCAATGTCGTTAACTCTACGTTGAATAGCATTGACGAAATCTGCGGAGAACTGCCATTCTGGCTACAGCAAGCGTTTATCGATCAGTTAGTTGCTGCCgg                                                                                                      <  2:171675/120‑1 (MQ=255)
                                                                      tGAAAAACGCGGTGACTCGCTCAATGTCCATAACTCGCCGATCAATACCAGTGACGAAAGCGGCGGATAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGATACATGCTGCCgg                                                                                                      >  1:171675/1‑120 (MQ=255)
                                                                         aaaaCGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGgtggcgtgg                                                                       <  2:43891/148‑1 (MQ=255)
                                                                                  tGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGc                                                             >  3:292805/1‑149 (MQ=255)
                                                                                  tGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGc                                                             <  8:248278/149‑1 (MQ=255)
                                                                                      tCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGa                                                          >  2:136430/1‑148 (MQ=255)
                                                                                                tCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGgctgc                                                                   <  2:241235/129‑1 (MQ=255)
                                                                                                tCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGgctgc                                                                   >  1:241235/1‑129 (MQ=255)
                                                                                                     aaCTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGc                                           <  4:233932/148‑1 (MQ=255)
                                                                                                        tCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGc                                        >  6:186791/1‑148 (MQ=255)
                                                                                                             gTTCAATAGCAGTGACCAAAGCGGCGGAGAACAGCCATTCTGGCACCAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCACCTAacac                                  >  5:149066/1‑149 (MQ=255)
                                                                                                               tCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGt                                >  2:327466/1‑149 (MQ=255)
                                                                                                                       aGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCt                          >  8:38433/1‑147 (MQ=255)
                                                                                                                        gTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAg                       >  1:164836/1‑149 (MQ=255)
                                                                                                                                         gaACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACa      <  7:17815/149‑1 (MQ=255)
                                                                                                                                          aaCTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAg     >  8:45738/1‑149 (MQ=255)
                                                                                                                                               ccATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAGCaa  >  5:322708/1‑147 (MQ=255)
                                                                                                                                               ccATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACAGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAGCaa  >  3:200331/1‑147 (MQ=255)
                                                                                                                                               ccAGTCTGGCAACCGCAAGCGTGTATCGATCAGGTACTTGCTGCCGGTCGCTGGTGGCGGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAGCaa  <  4:276159/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGGCTTTTCTGAAAAACGCGGTGACTCGCTCAATGTCGTTAACTCGCCGTTCAATAGCAGTGACGAAAGCGGCGGAGAACTGCCATTCTGGCAACAGCAAGCGTTTATCGATCAGTTACTTGCTGCCGGTCGCTGGTTGCTGGTACTGCTGGTGGCGTGGCTGCTGTGGCGGAAAGCGGTACGTCCGCAGCTAACACGTCGCGCTGAGGCGATGAAAGCTGTACAGCAA  >  NC_000913/2014403‑2014692

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