Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,577,897 T→C 41.5% I14V (ATA→GTA)  xapA ← purine nucleoside phosphorylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,577,8970TC41.5% 19.8 / 46.6 65I14V (ATA→GTA) xapApurine nucleoside phosphorylase II
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (9/18);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCT  >  minE/1577845‑1577965
                                                     |                                                                    
aTAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGt                                                        <  1:2468356/67‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGt                                                        <  1:675147/61‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:172741/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2433898/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:248808/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2221754/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2216389/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2096085/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2056593/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2594049/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:259967/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:2865268/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:1830426/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:3044591/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:962850/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:54273/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:449275/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:352833/69‑1 (MQ=255)
      gCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:1702406/69‑1 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAAGAGCGGGTTATGAGAAAATTg                                           >  1:1276208/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAAtt                                            >  1:106172/1‑68 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAAtt                                            >  1:3107760/1‑68 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:442403/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:1715939/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:462521/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:3190125/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:799191/1‑69 (MQ=255)
           tCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTg                                           >  1:3226250/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:3297497/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:1717398/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:1520568/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:2718016/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:1878392/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:849915/69‑1 (MQ=255)
                         tCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              <  1:1116810/69‑1 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAAc                                                               >  1:1878307/1‑35 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              >  1:2614035/1‑68 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              >  1:3057044/1‑68 (MQ=255)
                          cAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGacata                              >  1:2527048/1‑68 (MQ=255)
                           aGGTTTATAAGTCTTGATAATATCTACGCAATACAGTGGGTTATGAGaa                                                <  1:1248341/49‑1 (MQ=255)
                                              atatCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAggg          >  1:1074943/1‑68 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATacatatat                           >  1:1904272/1‑50 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATacatatat                           >  1:2263895/1‑50 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:1712399/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:1587245/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:1535885/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:3011842/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:615637/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:1448237/1‑51 (MQ=255)
                                               tatCTACGCAATACAGTGGGTTATGAGAAAACTGAGCCTGATAcatatata                          >  1:2410769/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)
                                                     |                                                                    
ATAAAGGCCACTCGTGGCGTGAAATCAGGTTTATAAGTCTTGATAATATCTATGCAAAACAGCGGGTTATGAGAAAATTGAACCTGAGACATACATATCCTTTTTCTGTAGGGTGGAATCT  >  minE/1577845‑1577965

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: