Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,152,016 C→T 13.2% E538K (GAG→AAG)  yhfK ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,152,0160CT13.2% 181.4 / 16.8 68E538K (GAG→AAG) yhfKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (21/38);  new base T (4/5);  total (25/43)
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.91e-01

ACAGAGTGTTATGTGCCTGATTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAC  >  W3110S.gb/4151950‑4152083
                                                                  |                                                                   
aCAGAGTGTTATGTGCCTGATTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCt                                                                  >  1:2527480/1‑70 (MQ=255)
aCAGAGTGTTATGTGCCTGATTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCt                                                                  >  1:318113/1‑70 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCTAATCGc                                                <  1:3218729/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:496810/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:2339522/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:3196441/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:1211042/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:1956431/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1309073/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:2129403/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1929530/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:3328783/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1704309/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:2987415/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:899531/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:2192582/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:3098647/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1492764/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:77908/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1242398/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:638259/71‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGGTGCGGATCCTCGCTAAGAATCAAGCGAATCGc                                                <  1:1526756/69‑1 (MQ=255)
                   aTTTACCCGCATTCGCTGCCAGCCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1528077/71‑1 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:2769402/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:91089/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:749047/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:3253659/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:1136981/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:3300875/1‑71 (MQ=255)
                    tttACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTc                                             >  1:3304800/1‑71 (MQ=255)
                       aCCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCt                                           <  1:3318198/70‑1 (MQ=255)
                       aCCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCt                                           <  1:1193536/70‑1 (MQ=255)
                        cccGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:885709/66‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTcc                                            <  1:1517923/62‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTcc                                            <  1:946771/62‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTcc                                            <  1:1168561/62‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTcc                                            <  1:3230556/62‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTcc                                            <  1:939772/62‑1 (MQ=255)
                              ttCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCt                                           <  1:35369/63‑1 (MQ=255)
                                cGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTCCTg                                          >  1:305079/1‑62 (MQ=255)
                                cGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTCCTg                                          >  1:3080792/1‑62 (MQ=255)
                                cGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTCCTg                                          >  1:1188908/1‑62 (MQ=255)
                                cGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTTGCTAAGAATCAAGCGAATCGCTTCCTg                                          >  1:1099576/1‑62 (MQ=255)
                                    gccaggccagTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1666311/54‑1 (MQ=255)
                                                gCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGc                  <  1:645775/70‑1 (MQ=255)
                                                gCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGc                  <  1:2908067/70‑1 (MQ=255)
                                                gCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGc                  <  1:3289000/70‑1 (MQ=255)
                                                gCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGc                  <  1:1117095/70‑1 (MQ=255)
                                                  gTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCtt                                              <  1:1757016/40‑1 (MQ=255)
                                                     ggTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGtttt             >  1:354828/1‑70 (MQ=255)
                                                     ggTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGtttt             >  1:2888400/1‑70 (MQ=255)
                                                     ggTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGtttt             >  1:2880532/1‑70 (MQ=255)
                                                     ggTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGtttt             >  1:2198729/1‑70 (MQ=255)
                                                      gTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGa                                         <  1:2922294/41‑1 (MQ=255)
                                                      gTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGa                                         <  1:915567/41‑1 (MQ=255)
                                                      gTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGa                                         <  1:580942/41‑1 (MQ=255)
                                                      gTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGa                                         <  1:3311852/41‑1 (MQ=255)
                                                      gTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGa                                         <  1:2231430/41‑1 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:2391991/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:1109220/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:1299682/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:254644/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:2701765/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:2895980/1‑71 (MQ=255)
                                                               ccTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:3327489/1‑71 (MQ=255)
                                                               acTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  >  1:1629843/2‑71 (MQ=255)
                                                                cTCGCTAAGAATCAAGCTAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  <  1:2217981/70‑1 (MQ=255)
                                                                 tCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAc  <  1:122267/69‑1 (MQ=255)
                                                                  |                                                                   
ACAGAGTGTTATGTGCCTGATTTACCCGCATTCGCTGCCAGGCCAGTGGCGTAGGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAAC  >  W3110S.gb/4151950‑4152083

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