Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,577,902 A→T 39.1% F12Y (TTT→TAT)  xapA ← purine nucleoside phosphorylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,577,9020AT39.1% 41.3 / 62.8 69F12Y (TTT→TAT) xapApurine nucleoside phosphorylase II
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base T (9/18);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

ATAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTAA  >  minE/1577845‑1577967
                                                          |                                                                 
aTAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGt                                                          <  1:2468356/67‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGt                                                          <  1:675147/61‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:352833/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2594049/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:1830426/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:259967/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2865268/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:3044591/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:172741/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:248808/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:1702406/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2433898/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2056593/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:449275/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2096085/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:54273/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2216389/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:2221754/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:962850/69‑1 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAAGAGCGGGTTATGAGAAAATTg                                             >  1:1276208/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAAtt                                              >  1:3107760/1‑68 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAAtt                                              >  1:106172/1‑68 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:1715939/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:3190125/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:3226250/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:442403/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:462521/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                             >  1:799191/1‑69 (MQ=255)
            cGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATTCAAAACAGCGGGTTATGAGAAAATTGa                                            <  1:736637/69‑1 (MQ=255)
               ggCGTGAAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGaa                                           <  1:2316559/68‑1 (MQ=255)
               ggCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCaaa                                                                   >  1:966154/1‑43 (MQ=255)
               ggCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAAc                                          >  1:1626320/1‑68 (MQ=255)
               ggCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAAc                                          >  1:1802059/1‑68 (MQ=255)
                gCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCt                                        <  1:1764521/69‑1 (MQ=255)
                gCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCt                                        <  1:1768510/69‑1 (MQ=255)
                  gTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTga                                      <  1:861769/69‑1 (MQ=255)
                    gAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTgaga                                    <  1:3064244/69‑1 (MQ=255)
                      aaaTCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTgaga                                    <  1:559471/68‑1 (MQ=255)
                        aTCAGGTTTATAAGTCTTGATAATATATATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGaca                                  <  1:2412777/68‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:1878392/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:3297497/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:2718016/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:1116810/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:849915/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:1520568/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                <  1:1717398/69‑1 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAAc                                                                 >  1:1878307/1‑35 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:2527048/1‑68 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:2614035/1‑68 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:3057044/1‑68 (MQ=255)
                           aGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                  <  1:1248341/49‑1 (MQ=255)
                                              atatCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAggg            >  1:1074943/1‑68 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATacatatat                             >  1:2263895/1‑50 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATacatatat                             >  1:1904272/1‑50 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:1712399/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:3011842/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:1587245/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:615637/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:1535885/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:2410769/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                            >  1:1448237/1‑51 (MQ=255)
                                                 tCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:3097171/1‑45 (MQ=255)
                                                 tCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:2942010/1‑45 (MQ=255)
                                                 tCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                                >  1:2560810/1‑45 (MQ=255)
                                                 tCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGa        >  1:1333117/1‑69 (MQ=255)
                                                     tGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCt    >  1:2408695/1‑69 (MQ=255)
                                                       cAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTaa  >  1:1575435/1‑69 (MQ=255)
                                                       cAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTaa  >  1:3105313/1‑69 (MQ=255)
                                                       cAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTaa  >  1:1135216/1‑69 (MQ=255)
                                                       cAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTaa  >  1:1935356/1‑69 (MQ=255)
                                                          |                                                                 
ATAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCTAA  >  minE/1577845‑1577967

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: