Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,963,758 A→G 7.5% D431D (GAT→GAC)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,963,7580AG 7.5% 157.8 / 4.5 53D431D (GAT→GAC) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base G (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

TGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAAC  >  NC_000913/1963613‑1963905
                                                                                                                                                 |                                                                                                                                                   
tgtACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATca                                                                                                                                                  <  2:396514/149‑1 (MQ=255)
  tACGCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATcacc                                                                                                                                                >  1:500273/1‑149 (MQ=255)
    ccccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATcaccac                                                                                                                                              <  2:43181/149‑1 (MQ=255)
      ccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTg                                                                                                                                            >  2:476120/1‑149 (MQ=255)
           aTCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAAt                                                                                                                                       <  1:176620/149‑1 (MQ=255)
               gcgcCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTcc                                                                                                                                   <  1:565374/149‑1 (MQ=255)
                 gcCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCAc                                                                                                                                 <  3:500407/149‑1 (MQ=255)
                     gTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCt                                                                                                                             >  1:367688/1‑149 (MQ=255)
                     gTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGATCCAGGGTTAAGCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCt                                                                                                                             >  3:425229/1‑149 (MQ=255)
                        cTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtc                                                                                                                          <  4:300135/149‑1 (MQ=255)
                        cTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtc                                                                                                                          <  1:488683/149‑1 (MQ=255)
                         ttttAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCGCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtca                                                                                                                         >  3:372826/1‑149 (MQ=255)
                         ttttAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtca                                                                                                                         >  1:12732/1‑149 (MQ=255)
                         ttttAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtca                                                                                                                         >  4:499947/1‑149 (MQ=255)
                               cgcACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCaaaa                                                                                                                   >  2:421355/1‑149 (MQ=255)
                                  aCTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGc                                                                                                                <  1:460618/149‑1 (MQ=255)
                                  aCTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGc                                                                                                                >  3:547692/1‑149 (MQ=255)
                                           ccAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCggg                                                                                                       <  3:197798/149‑1 (MQ=255)
                                              gATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAg                                                                                                    >  1:233345/1‑149 (MQ=255)
                                               aTAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAgg                                                                                                   >  4:97248/1‑149 (MQ=255)
                                                  gCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTg                                                                                                <  1:410440/149‑1 (MQ=255)
                                                    aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCGTCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCa                                                                                              >  2:340512/1‑149 (MQ=255)
                                                    aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCGTCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCa                                                                                              >  2:250934/1‑149 (MQ=255)
                                                    aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCa                                                                                              <  3:170039/149‑1 (MQ=255)
                                                    aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCa                                                                                              >  2:43094/1‑149 (MQ=255)
                                                     ttCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAg                                                                                             >  4:71875/1‑149 (MQ=255)
                                                      tCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGa                                                                                            <  1:314883/149‑1 (MQ=255)
                                                       cAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATATCGGGCAGGTTGCAGAt                                                                                           <  2:500273/149‑1 (MQ=255)
                                                       cAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGAt                                                                                           <  3:71875/149‑1 (MQ=255)
                                                       cAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGAGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGAt                                                                                           <  2:455156/149‑1 (MQ=255)
                                                                    ggATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTcgcg                                                                              >  4:135920/1‑149 (MQ=255)
                                                                      aTCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGaa                                                                            <  3:135920/149‑1 (MQ=255)
                                                                             gTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGGTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCTCGCCTTTCATCAAAATGCGATAGCGGGCAGGGTGCAGATCCATATTGTCGCGAATGTGcac                                                                     <  1:98579/149‑1 (MQ=255)
                                                                                     aCCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTggcg                                                             >  1:458122/1‑149 (MQ=255)
                                                                                            aaCGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGaaa                                                      >  2:450773/1‑149 (MQ=255)
                                                                                              cGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCGTCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATc                                                    <  1:340512/149‑1 (MQ=255)
                                                                                              cGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCGTCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATc                                                    <  1:250934/149‑1 (MQ=255)
                                                                                               gTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATcc                                                   >  3:270713/1‑149 (MQ=255)
                                                                                                       aGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCt                                           <  2:36743/149‑1 (MQ=255)
                                                                                                        gCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTg                                          >  2:448101/1‑149 (MQ=255)
                                                                                                                 gggTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAAttt                                 >  4:387320/1‑149 (MQ=255)
                                                                                                                    ttAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCtt                              <  3:97248/149‑1 (MQ=255)
                                                                                                                              gccagcGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTg                    <  1:12209/149‑1 (MQ=255)
                                                                                                                               ccagcGCCCCTGAAAAGCAACACCACCTGCAATCCCCCCGCCTTTCCTCAAAATGCGAAGGGGGGGAGGTTGCAGATCCATATTGTAGCTAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGc                   <  1:567538/149‑1 (MQ=255)
                                                                                                                               ccagcGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGc                   <  1:476120/149‑1 (MQ=255)
                                                                                                                               ccagcGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGc                   <  3:30907/149‑1 (MQ=255)
                                                                                                                                       ccGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACgg           >  2:306563/1‑149 (MQ=255)
                                                                                                                                         ggATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGcc         <  2:425005/149‑1 (MQ=255)
                                                                                                                                            tAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCcaa      <  1:306563/149‑1 (MQ=255)
                                                                                                                                               gCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCcaacaa   >  1:301364/1‑149 (MQ=255)
                                                                                                                                               gCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCcaacaa   <  3:30244/149‑1 (MQ=255)
                                                                                                                                                cATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCaacaac  <  2:367688/149‑1 (MQ=255)
                                                                                                                                                cATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCaacaac  >  2:272350/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGGCAAATTTCTTGCGGATACTGCGTATACGGCCCAACAAC  >  NC_000913/1963613‑1963905

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑