Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,983,125 T→C 7.7% S5G (AGT→GGT)  flgC ← flagellar basal‑body rod protein FlgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,983,1250TC 7.7% 144.5 / 2.5 52S5G (AGT→GGT) flgCflagellar basal‑body rod protein FlgC
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base C (2/2);  total (23/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 = 7.71e-01

GAACACCGGGTGGCGGGCGCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTGGTCGATC  >  NC_002947/4982991‑4983272
                                                                                                                                       |                                                                                                                                                   
gAACACCGGGTGGCGGGCGCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAgg                                                                                                                                              >  2:45717/1‑142 (MQ=255)
gAACACCGGGTGGCGGGCGCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAgg                                                                                                                                              <  1:45717/142‑1 (MQ=255)
gAACACCGGGTGGCGGGCGCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATgg                                                                                                                                       <  2:285534/149‑1 (MQ=255)
              gggcgCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCg                                                                                                                           >  1:215699/1‑147 (MQ=255)
                  ccgGTACGTCTGGTCGATGCTCTAAGTCACTGTGTCGGCGTTGGCGATGTTGGACGCCACTGTGTTCAGACGCGTGTTTTGTGGGCTCATGCCGCTACGGGCAATG‑TTGTAGACTGTGGAAAGGGACATGGTCTTTCTCCCCTCAgggg                                                                                                                     <  8:169498/148‑2 (MQ=255)
                      tACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGa                                                                                                                  >  8:299713/1‑148 (MQ=255)
                         gTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAAcc                                                                                                              <  1:129520/149‑1 (MQ=255)
                          tCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAAc                                                                                                               >  2:79733/1‑147 (MQ=255)
                             ggTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGcc                                                                                                          >  8:232457/1‑149 (MQ=255)
                              gTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCAGGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGACAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAg                                                                                                            <  8:181841/146‑1 (MQ=255)
                              gTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGTAAAGGGACATGGTCATTCTCCCCAAATTGAGGAAACCAGccc                                                                                                         >  8:81617/1‑149 (MQ=255)
                              gTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGcc                                                                                                          <  7:299713/148‑1 (MQ=255)
                                    gCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACCGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGTCAAGCATCCCTTTGaa                                                                                                   >  8:239241/1‑149 (MQ=255)
                                    gCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACCGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGaa                                                                                                   >  8:130941/1‑149 (MQ=255)
                                     cTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATGCTTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGa                                                                                                    <  1:259841/148‑1 (MQ=255)
                                          aGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTc                                                                                                                                     >  7:185110/1‑109 (MQ=255)
                                          aGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTc                                                                                                                                     <  8:185110/109‑1 (MQ=255)
                                           gACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTg                                                                                            <  1:79733/149‑1 (MQ=255)
                                             cacGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGtt                                                                                          <  4:293495/149‑1 (MQ=255)
                                               cGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGa                                                                                        <  8:82175/149‑1 (MQ=255)
                                                      ggCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAg                                                                                   <  2:92464/147‑1 (MQ=255)
                                                      ggCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAg                                                                                   <  4:161069/147‑1 (MQ=255)
                                                                    tGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGACCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTgg                                                                    >  3:291219/1‑148 (MQ=255)
                                                                     ccaCGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGaa                                                                  <  1:117489/147‑1 (MQ=255)
                                                                                    cAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGTAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTc                                                   >  6:261429/1‑149 (MQ=255)
                                                                                    cAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTc                                                   >  4:177221/1‑149 (MQ=255)
                                                                                     aGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCCCAGGGCGGAAACAAGCCCTTTGAATTTACAGTTGAGCAAGGAGAAGCTGGCCTGGAATCCAATAGCGATTTcc                                                  >  2:246166/1‑149 (MQ=255)
                                                                                       aCGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGaa                                                                  >  8:276986/1‑131 (MQ=255)
                                                                                       aCGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGaa                                                                  <  7:276986/131‑1 (MQ=255)
                                                                                       aCGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCgt                                                >  7:18063/1‑149 (MQ=255)
                                                                                         gcgTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCgtg                                               >  3:122556/1‑148 (MQ=255)
                                                                                          cgTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTa                                             <  3:177221/149‑1 (MQ=255)
                                                                                          cgTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTa                                             >  1:177846/1‑149 (MQ=255)
                                                                                              ttCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGTAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGtt                                          <  5:261429/148‑1 (MQ=255)
                                                                                                   tgCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACCGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGAtt                                     <  7:130941/148‑1 (MQ=255)
                                                                                                   tgCGCTAATGCCGCTCCCGGCAATG‑TTGAAGACACCGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGAtt                                     <  7:239241/148‑1 (MQ=255)
                                                                                                    tcaacattGCCGGTAGCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGt                                   >  1:59515/8‑149 (MQ=255)
                                                                                                     cgcTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGtt                                  <  7:232457/149‑1 (MQ=255)
                                                                                                     cgcTCATGCAGCTACCGGCAATG‑TTGAATACACTTGAAAGGGAGATGGTCATTCTCCCCGCAGGTCTGAAACCAGCCCTTTTAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGtt                                  <  1:246166/149‑1 (MQ=255)
                                                                                                           tGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCt                            <  4:122556/149‑1 (MQ=255)
                                                                                                             ccGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTg                           <  2:215699/148‑1 (MQ=255)
                                                                                                                cTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCa                       >  8:23131/1‑149 (MQ=255)
                                                                                                                             ttGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAAt                                                            >  6:165020/1‑100 (MQ=255)
                                                                                                                             ttGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAAt                                                            <  5:165020/100‑1 (MQ=255)
                                                                                                                             ttGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTc           <  2:177846/149‑1 (MQ=255)
                                                                                                                                aaGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTgg        >  5:104521/1‑149 (MQ=255)
                                                                                                                                aaGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTgg        >  5:208476/1‑149 (MQ=255)
                                                                                                                                aaGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTg         <  6:92359/148‑1 (MQ=255)
                                                                                                                                 aGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTGGt       <  8:18063/149‑1 (MQ=255)
                                                                                                                                 aGACACTGGAAAGGGACATGGTCATTCTCCCCACAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTGGt       <  7:81617/149‑1 (MQ=255)
                                                                                                                                     acTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTgc                          <  1:72318/126‑1 (MQ=255)
                                                                                                                                     acTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTgc                          >  2:72318/1‑126 (MQ=255)
                                                                                                                                       tGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTGGtcgatc  >  3:287725/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
GAACACCGGGTGGCGGGCGCGGTACGTCTGGTCGATGCTCGAAGACACGGTCTCGGCGTTGGCGATGTTGGACGCCACGGTGTTCAGACGCGTGTTCTGTGCGCTCATGCCGCTACCGGCAATG‑TTGAAGACACTGGAAAGGGACATGGTCATTCTCCCCGCAGGGCGGAAACCAGCCCTTTGAATTTACTGTTGAGCAAGGTGAAGCTGGCCTGGAAGCCAATCGCGTTTTCCGTGTAGTTCGATTGTTCGATCTGCGCATCCACGGTGTTCTGGTCGATC  >  NC_002947/4982991‑4983272

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