Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 240,219 A→C 13.6% V244V (GTT→GTG araJ ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE240,2190AC13.6% 160.9 / 17.9 81V244V (GTT→GTGaraJpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (41/29);  new base C (5/6);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CGCAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTT  >  minE/240154‑240284
                                                                 |                                                                 
cgcAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGTCCAACt                                                                <  1:440858/69‑1 (MQ=255)
cgcAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACt                                                                <  1:1864938/69‑1 (MQ=255)
cgcAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACt                                                                <  1:1168261/69‑1 (MQ=255)
  cAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACTaa                                                              <  1:3412908/69‑1 (MQ=255)
  cAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACTaa                                                              <  1:2953962/69‑1 (MQ=255)
       gtgAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACTAAcatc                                                          >  1:1615891/1‑68 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:1006341/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:2607924/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:2488473/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:921077/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:2458212/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:1157579/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:698136/1‑69 (MQ=255)
                                   tAGCATATTTCCCAGCACCATCCCTAGCACAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                             >  1:617143/1‑69 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGt                                      >  1:2182530/1‑53 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACGGGa                       >  1:2080691/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCgg                        >  1:2252229/1‑67 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:619077/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:628333/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:2502342/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:2935576/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:1146279/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:2619823/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:2685907/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:1489269/1‑68 (MQ=255)
                                          tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                       >  1:530439/1‑68 (MQ=255)
                                          ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGaaaa                             >  1:651088/4‑62 (MQ=255)
                                          ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                       >  1:1525415/4‑68 (MQ=255)
                                          ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                       >  1:1698240/4‑68 (MQ=255)
                                          ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                       >  1:1724990/4‑68 (MQ=255)
                                          ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                       >  1:510347/4‑68 (MQ=255)
                                            ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:881905/68‑1 (MQ=255)
                                            ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:686590/68‑1 (MQ=255)
                                            ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:1126642/68‑1 (MQ=255)
                                            ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:779293/68‑1 (MQ=255)
                                            ccccAGCACAATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:3065101/68‑1 (MQ=255)
                                              ccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                    <  1:2942522/67‑1 (MQ=255)
                                              ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                   <  1:1307242/68‑1 (MQ=255)
                                              ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                   <  1:2424924/68‑1 (MQ=255)
                                              ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                   <  1:1879808/68‑1 (MQ=255)
                                               cAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                   <  1:1781244/67‑1 (MQ=255)
                                                aGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                 <  1:2417387/68‑1 (MQ=255)
                                                aGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                 <  1:1496303/68‑1 (MQ=255)
                                                aGCACCATCCCTAGCCCAACTAACATCATAAAAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                 <  1:148944/68‑1 (MQ=255)
                                                  cACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAAcat              <  1:2718540/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:1587379/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:700509/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2938120/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2893170/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2764669/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2749849/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:51405/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:1438787/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2261442/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:3230674/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:2273244/69‑1 (MQ=255)
                                                      aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:1310924/69‑1 (MQ=255)
                                                       tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:3335219/68‑1 (MQ=255)
                                                       tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:1937052/68‑1 (MQ=255)
                                                       tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:951868/68‑1 (MQ=255)
                                                       tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg          <  1:1068730/68‑1 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGt         >  1:969767/1‑67 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa        >  1:1223642/1‑68 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa        >  1:2849537/1‑68 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa        >  1:1843854/1‑68 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa        >  1:2634043/1‑68 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa        >  1:1851211/1‑68 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:400652/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:742400/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:1984740/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2270848/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2390502/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2492293/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2494838/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2677270/1‑69 (MQ=255)
                                                         ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt       >  1:2953455/1‑69 (MQ=255)
                                                            aGCCCAACTAACATCATAATAAAGGTCATCGCCGttt                                    >  1:325695/1‑37 (MQ=255)
                                                              cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt  >  1:1707455/1‑69 (MQ=255)
                                                              cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt  >  1:13471/1‑69 (MQ=255)
                                                              cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt  >  1:11908/1‑69 (MQ=255)
                                                              cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACAGGAAATAAACATCATGTATGGCtt  >  1:1366423/1‑69 (MQ=255)
                                                                 |                                                                 
CGCAGTGGTGAATAACGTCCTGAAATCCTGCCACTTAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTT  >  minE/240154‑240284

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: