Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 430,884 G→A 7.7% F118F (TTC→TTT tsx ← nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb430,8840GA 7.7% 229.9 / 6.4 78F118F (TTC→TTTtsxnucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base A (3/3);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGAA  >  W3110S.gb/430815‑430953
                                                                     |                                                                     
gCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCg                                                                       >  1:3597762/1‑70 (MQ=255)
gCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCg                                                                       >  1:1644421/1‑70 (MQ=255)
gCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGa                                                                      <  1:2877558/71‑1 (MQ=255)
    ccATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCt                                                                  <  1:1135120/71‑1 (MQ=255)
    ccATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCt                                                                  <  1:3415242/71‑1 (MQ=255)
     cATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCt                                                                  >  1:2596360/1‑70 (MQ=255)
     cATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTa                                                                 <  1:2707744/71‑1 (MQ=255)
     cATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTa                                                                 <  1:17698/71‑1 (MQ=255)
       tCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAg                                                               >  1:2860673/1‑71 (MQ=255)
       tCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAg                                                               >  1:3651079/1‑71 (MQ=255)
          ttattaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             >  1:2360529/1‑70 (MQ=255)
          ttattaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             >  1:2963849/1‑70 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTa                                                                 <  1:2264647/63‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             <  1:1086547/67‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             <  1:2676407/67‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             <  1:3254698/67‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGt                                                             <  1:3796520/67‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          <  1:776118/70‑1 (MQ=255)
             ttaCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          <  1:3776494/70‑1 (MQ=255)
               aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCa                                                           >  1:2368685/1‑67 (MQ=255)
               aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCa                                                           >  1:2822113/1‑67 (MQ=255)
                 gACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTAtt                                                      >  1:1794811/1‑70 (MQ=255)
                 gACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                                     <  1:2560377/71‑1 (MQ=255)
                   cccATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGt                                                   <  1:2110693/71‑1 (MQ=255)
                        gTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTc                                                            <  1:1664783/57‑1 (MQ=255)
                         tCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTag                                                                >  1:1214339/1‑51 (MQ=255)
                          cGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          >  1:3297281/1‑57 (MQ=255)
                          cGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          >  1:1491453/1‑57 (MQ=255)
                           gTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTg                                            <  1:3683807/70‑1 (MQ=255)
                           gTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGt                                           >  1:634894/1‑71 (MQ=255)
                           gTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGt                                           >  1:1168258/1‑71 (MQ=255)
                            tAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAg                                                <  1:311610/65‑1 (MQ=255)
                             aaaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                                          >  1:2857941/1‑70 (MQ=255)
                             aaaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                                          >  1:699138/1‑70 (MQ=255)
                             aaaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                                          >  1:923218/1‑70 (MQ=255)
                                tGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        >  1:154884/1‑69 (MQ=255)
                                tGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                                      <  1:512619/71‑1 (MQ=255)
                                  tAGTTGTTCGCGAAGTACCACTCTTTGAACTGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                                      >  1:2923445/1‑69 (MQ=255)
                                  tAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                                      >  1:819063/1‑69 (MQ=255)
                                  tAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                                      >  1:1800512/1‑69 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGgaga                                  <  1:2403322/71‑1 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGgaga                                  <  1:3241660/71‑1 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGgaga                                  <  1:3038494/71‑1 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        <  1:2565535/65‑1 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        <  1:922301/65‑1 (MQ=255)
                                    gttgttCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        <  1:1428912/65‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          <  1:2378082/46‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                                          <  1:2937207/46‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                                     <  1:3532931/51‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                                     <  1:1344032/51‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCt                                              <  1:2722068/58‑1 (MQ=255)
                                     ttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGaa                                 <  1:540600/71‑1 (MQ=255)
                                       gttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAAc                               >  1:3466214/1‑71 (MQ=255)
                                       gttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAAc                               >  1:1480538/1‑71 (MQ=255)
                                        ttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGgaga                                  <  1:2960891/67‑1 (MQ=255)
                                          cgcgAAGTACCACTCTTTGAACGGACCGAAGCTAAAGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGt                             <  1:1128284/70‑1 (MQ=255)
                                          cgcgAAGTAACACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTg                            >  1:2937628/1‑71 (MQ=255)
                                           gcgAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTgg                           <  1:1843837/71‑1 (MQ=255)
                                           gcgAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTg                            <  1:2910464/70‑1 (MQ=255)
                                              aaGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTc                        >  1:1126748/1‑71 (MQ=255)
                                              aaGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTc                        >  1:2508578/1‑71 (MQ=255)
                                              aaGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTAGTCGATGGAGAAACg                              >  1:3531939/1‑65 (MQ=255)
                                                  aCCACTCTTTGAACGGACCGAAGCTAAGCTCAGTATTGGTCAGCt                                              >  1:3525898/1‑45 (MQ=255)
                                                   ccACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGgaga                                  >  1:457620/1‑56 (MQ=255)
                                                   ccACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGAttt                   >  1:1626479/1‑71 (MQ=255)
                                                   ccACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGAttt                   >  1:2447689/1‑71 (MQ=255)
                                                   ccACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGAttt                   >  1:2345441/1‑71 (MQ=255)
                                                   ccACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGAtt                    <  1:1042420/70‑1 (MQ=255)
                                                      ctctTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTc                                          <  1:2105029/45‑1 (MQ=255)
                                                      ctctTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCg                                         <  1:389040/46‑1 (MQ=255)
                                                       tctTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAg                                                <  1:2773910/38‑1 (MQ=255)
                                                        ctTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        >  1:1493811/1‑45 (MQ=255)
                                                         tttGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa                      >  1:2164514/1‑62 (MQ=255)
                                                         tttGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATa              >  1:770225/1‑70 (MQ=255)
                                                         tttGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATa              >  1:3853829/1‑70 (MQ=255)
                                                           tGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATa              <  1:884286/68‑1 (MQ=255)
                                                              aCGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                                      <  1:3411882/41‑1 (MQ=255)
                                                               cGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGa                                        >  1:3310960/1‑38 (MQ=255)
                                                                  aCCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGc        >  1:218836/1‑67 (MQ=255)
                                                                    cGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGaa  >  1:40689/1‑71 (MQ=255)
                                                                     |                                                                     
GCGACCATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCCATAAACAGCGGAGAA  >  W3110S.gb/430815‑430953

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: