Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,518,147 C→T 8.1% R30H (CGC→CAC)  fecI ← RNA polymerase sigma‑19 factor, fec operon‑specific; ECF sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,1470CT 8.1% 147.5 / 4.2 50R30H (CGC→CAC) fecIRNA polymerase sigma‑19 factor, fec operon‑specific; ECF sigma factor
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (2/2);  total (20/30)
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.42e-01

GGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCAT  >  NC_000913/4518001‑4518294
                                                                                                                                                  |                                                                                                                                                   
ggCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGTAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTgcgc                                                                                                                                                   <  3:299711/149‑1 (MQ=255)
    aaaCAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTca                                                                                                                                               <  1:112213/149‑1 (MQ=255)
       cAGGCCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTcagcc                                                                                                                                            <  2:282311/149‑1 (MQ=255)
             caccaTCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCttt                                                                                                                                      <  3:215594/149‑1 (MQ=255)
              accaTCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCtttt                                                                                                                                     <  2:267928/149‑1 (MQ=255)
                  tCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAAc                                                                                                                                 <  4:187535/149‑1 (MQ=255)
                    aCGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCa                                                                                                                               <  3:188316/149‑1 (MQ=255)
                     cgcgTTTGGCGGTAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAg                                                                                                                              >  3:279219/1‑149 (MQ=255)
                           tGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTCTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGa                                                                                                                        <  4:287474/149‑1 (MQ=255)
                            ggCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGAt                                                                                                                       <  1:9427/149‑1 (MQ=255)
                                aTAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGAtgtgt                                                                                                                   <  1:261565/149‑1 (MQ=255)
                                         aggaaggaGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAg                                                                                                          <  3:32108/149‑1 (MQ=255)
                                                 gCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGTGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGa                                                                                                  >  2:203444/1‑149 (MQ=255)
                                                 gCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGa                                                                                                  >  1:87142/1‑149 (MQ=255)
                                                   gAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAAc                                                                                                <  1:303408/149‑1 (MQ=255)
                                                   gAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAAc                                                                                                <  2:87142/149‑1 (MQ=255)
                                                        tCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTaa                                                                                           <  1:255544/149‑1 (MQ=255)
                                                               tcgtcgAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGGGTGCCATAAAGCGACTCGAACGTTAAGGAAGCt                                                                                    >  2:313752/1‑149 (MQ=255)
                                                                  tcgAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCtgtg                                                                                 <  2:304968/149‑1 (MQ=255)
                                                                    gagaGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTACATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGCAGCTGTGGt                                                                               >  1:75112/1‑149 (MQ=255)
                                                                      gagCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAg                                                                             <  2:317557/149‑1 (MQ=255)
                                                                                   aCCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGa                                                                >  4:254949/1‑149 (MQ=255)
                                                                                      aTTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGTGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACAt                                                             >  2:367311/1‑149 (MQ=255)
                                                                                           ccGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGa                                                        >  1:136956/1‑149 (MQ=255)
                                                                                           ccGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGa                                                        >  4:200154/1‑149 (MQ=255)
                                                                                           ccGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGa                                                        >  1:166056/1‑149 (MQ=255)
                                                                                            cGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTCATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAg                                                       >  3:110014/1‑149 (MQ=255)
                                                                                               aaaaaaGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGc                                                    <  2:166056/149‑1 (MQ=255)
                                                                                                 aaaaGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCtg                                                                                   <  3:415212/116‑1 (MQ=255)
                                                                                                 aaaaGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCtg                                                                                   >  4:415212/1‑116 (MQ=255)
                                                                                                      tgtCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTCATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCaaaa                                             <  4:110014/149‑1 (MQ=255)
                                                                                                        tCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGTGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGt                                           >  2:57666/1‑149 (MQ=255)
                                                                                                         ccTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGAATCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGtt                                          <  4:279219/149‑1 (MQ=255)
                                                                                                           tGGGCAATGATATATGCATCAAAATCAAACTGGAGTTTGAGAGCTAGCCAGCTTTTCAACAAGCCGTGATGTGTGCCGTAATGCGTCTCGATCGTTAAGGAAGCTGTGGTAGTGGCGCTTTCAGACATGCGGAGTGCTTCAAAAGTTaa                                        <  4:444652/149‑1 (MQ=255)
                                                                                                                     cATCTGCATCAAAAGCAGACTGGAGTTTGTGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTa                              <  1:203444/149‑1 (MQ=255)
                                                                                                                       tCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAg                                                                                          <  2:393918/87‑1 (MQ=255)
                                                                                                                       tCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAg                                                                                          >  1:393918/1‑87 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGcc                                                                                                                >  4:466246/1‑59 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGcc                                                                                                                <  3:466246/59‑1 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGtt                                          >  4:227896/1‑129 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGtt                                          <  3:227896/129‑1 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATAtt                      >  3:463516/1‑149 (MQ=255)
                                                                                                                             tCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGGGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATAtt                      >  2:4310/1‑149 (MQ=255)
                                                                                                                               aaaaGCAGACTGGAGTTTGTGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTaa                    <  1:57666/149‑1 (MQ=255)
                                                                                                                               aaaaGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTaa                    <  2:75112/149‑1 (MQ=255)
                                                                                                                                   gCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAGCGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTc                           >  3:440294/1‑138 (MQ=255)
                                                                                                                                   gCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAGCGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTc                           <  4:440294/138‑1 (MQ=255)
                                                                                                                                           gAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGt        >  1:178443/1‑149 (MQ=255)
                                                                                                                                              tttGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGTGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTAt     <  1:367311/149‑1 (MQ=255)
                                                                                                                                                tGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCCTTAAGCGACCCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGCTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGCATGGTTATCa   <  4:309615/149‑1 (MQ=255)
                                                                                                                                                tGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCa   <  1:423594/149‑1 (MQ=255)
                                                                                                                                                 gcgcGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCAt  <  2:136956/149‑1 (MQ=255)
                                                                                                                                                 gcgcGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCAt  <  3:254949/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCAT  >  NC_000913/4518001‑4518294

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