Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,217,545 A→G 12.6% N57N (AAT→AAC yjaB ← predicted acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,217,5450AG12.6% 211.1 / 19.5 79N57N (AAT→AACyjaBpredicted acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (47/22);  new base G (8/2);  total (55/24)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGA  >  W3110S.gb/4217497‑4217611
                                                |                                                                  
cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGgc                                               <  1:1528757/70‑1 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGa                                                          >  1:807998/1‑55 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:93796/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:2888568/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:447751/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:2273817/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:2344/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:1437890/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTc                                          >  1:327937/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACgg                                                >  1:2158906/1‑65 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTc                                          >  1:1720882/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTc                                          >  1:1561834/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTc                                          >  1:181357/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTc                                          >  1:1494267/1‑71 (MQ=255)
    tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTc                                          >  1:844411/1‑71 (MQ=255)
      ctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCg                                        >  1:1381287/1‑71 (MQ=255)
               tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaag                                >  1:1310387/1‑70 (MQ=255)
               tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaag                                >  1:2637259/1‑70 (MQ=255)
               tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaag                                >  1:684702/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:455489/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:1988922/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:172536/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:1961808/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:695204/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:3288434/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:3236387/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagg                               >  1:3171323/1‑70 (MQ=255)
                aGCAACATAAATCCAACCGGCTGATCCCGCTCAATaa                                                                >  1:1706183/1‑37 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:2234336/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:2614875/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:273111/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:99493/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1912220/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1221764/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1222776/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1279727/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:879929/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1539459/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:83372/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGa                                                          >  1:1875821/1‑41 (MQ=255)
                  cAACATAAATCCAACCGGCTGATCCCGCTCATTAAACGCGa                                                          >  1:376173/1‑41 (MQ=255)
                    aCATAAATCCAACCGGCTGCTCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacg                           >  1:3210507/1‑70 (MQ=255)
                    aCATAAATCCAACCGGCTGATCCCGCTCGTTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacga                          <  1:1442549/71‑1 (MQ=255)
                    aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGCCCCACAACGGCGCTTCCGGCAGGAAGgaacg                           >  1:2013956/1‑70 (MQ=255)
                    aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacg                           >  1:619786/1‑70 (MQ=255)
                    aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacg                           >  1:2785460/1‑70 (MQ=255)
                      aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:1542998/71‑1 (MQ=255)
                      aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:3212089/71‑1 (MQ=255)
                      aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:1265739/71‑1 (MQ=255)
                      aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:3207871/71‑1 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCAGAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:2979496/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:499543/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:1093974/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:3198437/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:3221117/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:797948/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:3049754/70‑1 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        >  1:2556228/1‑70 (MQ=255)
                       tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                        <  1:2030461/70‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:192768/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:1305437/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:1613050/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:738579/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:1616868/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:173119/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:1831/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:3187232/71‑1 (MQ=255)
                         aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:2244893/71‑1 (MQ=255)
                          aTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAgg                    <  1:2016843/71‑1 (MQ=255)
                                          ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc    <  1:1994939/71‑1 (MQ=255)
                                          ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc    <  1:2555017/71‑1 (MQ=255)
                                          ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc    <  1:1170849/71‑1 (MQ=255)
                                          ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc    <  1:1003682/71‑1 (MQ=255)
                                            gCTCATTAACCTCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTcc    <  1:1541013/69‑1 (MQ=255)
                                             cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGa  >  1:2359895/1‑70 (MQ=255)
                                             cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGa  >  1:2498773/1‑70 (MQ=255)
                                             cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGa  >  1:1593740/1‑70 (MQ=255)
                                             cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGa  >  1:1396267/1‑70 (MQ=255)
                                               cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                     <  1:2893506/49‑1 (MQ=255)
                                                |                                                                  
CATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGA  >  W3110S.gb/4217497‑4217611

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