Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 260,129 G→A 41.9% 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
*minE260,1290GA41.9% 35.2 / 67.0 93F118F (TTC→TTTtsxnucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (11/28);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 5.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCC  >  minE/260066‑260183
                                                               |                                                      
aTCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAg                                                     >  1:858429/1‑67 (MQ=255)
aTCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAg                                                     >  1:2550596/1‑67 (MQ=255)
aTCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAg                                                     >  1:3405341/1‑67 (MQ=255)
   tttATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTaa                                                 >  1:2619852/1‑68 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:3321102/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2126499/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2367326/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2398740/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1974729/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2460227/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2741816/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:2890019/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:3283244/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1591338/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1422684/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:470741/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:472150/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:128460/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1219225/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:988900/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1182272/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:723701/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:1051840/68‑1 (MQ=255)
         aCGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGGCAg                                           <  1:2237250/68‑1 (MQ=255)
          cGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:524091/67‑1 (MQ=255)
             cccATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      <  1:2509245/69‑1 (MQ=255)
             cccATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      <  1:2812025/69‑1 (MQ=255)
             cccATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      <  1:1690260/69‑1 (MQ=255)
             ccaaTGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      <  1:418902/66‑1 (MQ=255)
               cATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      >  1:2113820/1‑67 (MQ=255)
               cATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      >  1:1176893/1‑67 (MQ=255)
               cATGTCGTAAATGTAGTTGTTCGAGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      >  1:1872317/1‑67 (MQ=255)
                        aaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:660837/69‑1 (MQ=255)
                        aaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:933997/69‑1 (MQ=255)
                        aaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:3079355/69‑1 (MQ=255)
                        aaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:957283/69‑1 (MQ=255)
                        aaTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:1281689/69‑1 (MQ=255)
                         aTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:1520218/68‑1 (MQ=255)
                         aTGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTc                           <  1:2127926/68‑1 (MQ=255)
                            tAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTg                                      <  1:2015543/54‑1 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTgg                                     >  1:1363374/1‑53 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:729185/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:2723734/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:2623838/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:291124/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:3076695/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:395221/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:2594194/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:737716/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATgg                      >  1:2011946/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATg                       >  1:2015472/1‑67 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCACCTTGTCGATgg                      >  1:1966862/1‑68 (MQ=255)
                              gttgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGATCAGTATTGGTCAGCTTGTCGATgg                      >  1:3307298/1‑68 (MQ=255)
                                tgttCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAg                                           <  1:37175/45‑1 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAACCGTGGTTCGa       >  1:2494092/1‑69 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa       >  1:2799711/1‑69 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa       >  1:2330711/1‑69 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa       >  1:1883875/1‑69 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa       >  1:2956155/1‑69 (MQ=255)
                                            aCCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGa       >  1:2997326/1‑69 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGTTTCGAt      <  1:2070/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:371031/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:3288073/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:430296/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:449847/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:90608/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:203437/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:617650/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2407609/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:168918/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2915837/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:1921938/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2053986/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2742251/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2583310/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:2702921/69‑1 (MQ=255)
                                             ccACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATCGAGAAACGCGGTTCGAt      <  1:3124731/69‑1 (MQ=255)
                                              cACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:1623907/68‑1 (MQ=255)
                                              cACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:388624/68‑1 (MQ=255)
                                              cACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:1980913/68‑1 (MQ=255)
                                              cACTCTTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGAt      <  1:519589/68‑1 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:9905/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:618893/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:609489/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:2838634/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:2576087/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:2324120/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:2262480/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:1783649/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:1354297/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:1212373/1‑68 (MQ=255)
                                                 tctTTGAACGGACCAAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGCGGTTCGATTTc   >  1:1066145/1‑68 (MQ=255)
                                                  ctTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTcc  >  1:897619/1‑68 (MQ=255)
                                                               |                                                      
ATCTTTATTACGACCCATGTCGTAAATGTAGTTGTTCGCGAAGTACCACTCTTTGAACGGACCGAAGCTAAGGTCAGTATTGGTCAGCTTGTCGATGGAGAAACGTGGTTCGATTTCC  >  minE/260066‑260183

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: