Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,001,157 G→A G211G (GGC→GGT)  fliA ← RNA polymerase, sigma 28 (sigma F) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,1570GA98.2% 180.1 / ‑4.2 57G211G (GGC→GGT) fliARNA polymerase, sigma 28 (sigma F) factor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/27);  total (29/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 = 7.30e-01

GAACTCCTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCGCCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCAGTAGTTGTTGTA  >  NC_000913/2001028‑2001281
                                                                                                                                 |                                                                                                                            
gAACTCCTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCa                                                                                                                           <  1:603607/133‑1 (MQ=255)
 aaCTCCTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCaa                                                                                                                          >  2:621978/1‑133 (MQ=255)
   cTCCTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAAtc                                                                                                                        >  1:487496/1‑133 (MQ=255)
     ccTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCACCACCACACCAAtctc                                                                                                                      >  2:220366/1‑133 (MQ=255)
             tCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTC‑TTGAGAtt                                                                                                             >  2:520829/1‑133 (MQ=255)
              cAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGAGATTCCCCTACCTCCATCACCGCACTAATCTCTTTGAGAtt                                                                                                             >  2:199394/1‑133 (MQ=255)
                aaGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGGTTCCCCTAGCTCTAGCACCGCACCAATCTCTTTGAGATTCa                                                                                                           >  2:204594/1‑133 (MQ=255)
                aaGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGACTGGCTGTGTAACTGACTGACCCGCGATTCCCCGTCCTCCAGCACCGCTCCAATCTCTTTGAGATTCa                                                                                                           <  2:617495/133‑1 (MQ=255)
                 aGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAg                                                                                                          <  1:551991/133‑1 (MQ=255)
                      aaGTGCGGTATTTACTGATGTTATAACTTATCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCtt                                                                                                     <  1:288628/133‑1 (MQ=255)
                      aaGTGCGGCATTTACTGACGTTTTAACTTATCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCTACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCtt                                                                                                     <  1:204594/133‑1 (MQ=255)
                       aGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTc                                                                                                    <  2:152094/133‑1 (MQ=255)
                           cGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCtt                                                                                                     >  2:344929/1‑128 (MQ=255)
                           cGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCtt                                                                                                     <  1:344929/128‑1 (MQ=255)
                               aTTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAAt                                                                                            <  2:487496/133‑1 (MQ=255)
                               aTTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTC‑TTGAGATTCAGCTCTTCCTGGTAATa                                                                                           <  1:520829/133‑1 (MQ=255)
                                tttACTGACGTTATACCTTACCCACTTTCGTGCGTAACCGTGTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATa                                                                                           <  2:187364/133‑1 (MQ=255)
                                tttACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATa                                                                                           <  2:564131/133‑1 (MQ=255)
                                  tACTGACGTTATAACATACCCACTTTAGTGAGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATaga                                                                                         <  1:242772/133‑1 (MQ=255)
                                      gACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGt                                                                                     <  2:121497/133‑1 (MQ=255)
                                        cGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGGTTCCGCTCTTCCTGGTAATAGAGGGTTa                                                                                   >  2:248387/1‑133 (MQ=255)
                                                   aCCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGtttt                                                                        >  1:470856/1‑133 (MQ=255)
                                                   aCCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGtttt                                                                        >  2:215573/1‑133 (MQ=255)
                                                   aCCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGtttt                                                                        >  2:63858/1‑133 (MQ=255)
                                                       aGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATACAGCTCTTCCTGATAACAGAGGGTTTATACCAGTTTTTcgc                                                                    >  2:256835/1‑133 (MQ=255)
                                                         tttAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGATAATAGAGGGTCAAGACCAGTGTTTCgcgc                                                                  >  1:293609/1‑133 (MQ=255)
                                                          ttAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCt                                                                 >  1:414212/1‑133 (MQ=255)
                                                                 gTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCaa                                                          >  1:553083/1‑133 (MQ=255)
                                                                  tAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGTGCTCCGGCAAc                                                         >  1:483708/1‑133 (MQ=255)
                                                                        tttAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCg                                                   >  1:135871/1‑133 (MQ=255)
                                                                         tgaaTAGCCCGGCTGCGTAACTGTGTTATCCGTGATTCCCCCACCTCCAGCACCGCACCAACCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGa                                                  <  2:293609/131‑1 (MQ=255)
                                                                          tAATAGCCTGGCTGTGTATCTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGAt                                                 <  1:215573/133‑1 (MQ=255)
                                                                                ccTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTc                                           >  1:602176/1‑133 (MQ=255)
                                                                                ccTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTc                                           >  1:416978/1‑133 (MQ=255)
                                                                                   ggCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCAt                                        >  1:99472/1‑133 (MQ=255)
                                                                                   ggCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCAt                                        >  1:514914/1‑133 (MQ=255)
                                                                                    gCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATc                                       >  2:387442/1‑133 (MQ=255)
                                                                                         gtAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCg                                  >  2:458364/1‑133 (MQ=255)
                                                                                          tAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCGCCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGc                                 >  2:155712/1‑133 (MQ=255)
                                                                                                 ctgacCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCa                          <  2:414212/133‑1 (MQ=255)
                                                                                                 ctgacCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCa                          <  2:470856/133‑1 (MQ=255)
                                                                                                 ctgacCCGAGATTCCTCGTCCTCCAGCACCTCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCa                          <  1:256835/133‑1 (MQ=255)
                                                                                                  tgacCCGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAg                         >  1:579529/1‑133 (MQ=255)
                                                                                                      ccGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTa                     <  1:621978/133‑1 (MQ=255)
                                                                                                      ccGCGATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTa                     <  2:602176/133‑1 (MQ=255)
                                                                                                          gATTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGt                 >  1:446401/1‑133 (MQ=255)
                                                                                                           aTTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTc                >  1:57394/1‑133 (MQ=255)
                                                                                                           aTTCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTc                <  2:57394/133‑1 (MQ=255)
                                                                                                            ttCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTcc               <  1:387442/133‑1 (MQ=255)
                                                                                                             tCCCCGACCTCCAGCGCCACACCAATCTCTATGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCa              <  1:248387/133‑1 (MQ=255)
                                                                                                             tCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCa              <  2:135871/133‑1 (MQ=255)
                                                                                                             tCCCCGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCa              <  1:390394/133‑1 (MQ=255)
                                                                                                               cccGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCagt            <  1:458364/133‑1 (MQ=255)
                                                                                                                ccGACCTCCAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCagta           >  1:388270/1‑133 (MQ=255)
                                                                                                                        cAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCAGTAGTtgttgt   <  2:514914/133‑1 (MQ=255)
                                                                                                                        cAGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCAGTAGTtgttgt   <  2:99472/133‑1 (MQ=255)
                                                                                                                         aGCACCGCACCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCAGTAGTTGTTGTa  >  2:384620/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                            
GAACTCCTGGTAGTCAAAGTTAAAGTGCGGCATTTACTGACGTTATAACTTACCCAGTTTAGTGCGTAACCGTTTAATAGCCTGGCTGTGTAACTGACTGACCCGCGATTCCCCGACCTCCAGCACCGCGCCAATCTCTTTGAGATTCAGCTCTTCCTGGTAATAGAGGGTTAATACCAGTTTTTCGCGCTCCGGCAACGTTTCGATGGCTTCCATCACCCGCTGGCGCAGATTACTGTCCAGTAGTTGTTGTA  >  NC_000913/2001028‑2001281

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