Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,632,970 A→G 55.6% K421R (AAA→AGA)  livM → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,632,9700AG55.6% 8.1 / 31.4 27K421R (AAA→AGA) livMleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base G (6/9);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCCTGCTG  >  minE/2632900‑2633039
                                                                      |                                                                     
atgatCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGaa                                                                       <  1:1609891/71‑1 (MQ=255)
     cTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGagag                                                                     >  1:2020022/1‑68 (MQ=255)
        gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCgagaggagag                                                                >  1:1760883/1‑70 (MQ=255)
        gCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAgg                                                                    <  1:1793727/66‑1 (MQ=255)
           tCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAgg                                                            >  1:221127/1‑71 (MQ=255)
           tCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAgg                                                            <  1:2051450/71‑1 (MQ=255)
           tCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAgg                                                            <  1:1683645/71‑1 (MQ=255)
             cGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCa                                                          >  1:411041/1‑71 (MQ=255)
             cGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCa                                                          <  1:1309367/71‑1 (MQ=255)
                       gctgcCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGc                                                 >  1:2131236/1‑70 (MQ=255)
                            ccATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGACCAGGCATGAGTCAGCCattat                                           <  1:1005140/71‑1 (MQ=255)
                             cATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCattat                                           >  1:1594243/1‑70 (MQ=255)
                                  cgcgCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCa                                                              >  1:1194476/1‑46 (MQ=255)
                                     gcCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGtt                                   <  1:674347/70‑1 (MQ=255)
                                     gcCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGtt                                   <  1:1103875/70‑1 (MQ=255)
                                       ccGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                      <  1:1340427/65‑1 (MQ=255)
                                       ccGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                      <  1:1652878/65‑1 (MQ=255)
                                       ccGCAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                      <  1:1740768/65‑1 (MQ=255)
                                       ccGCAACTGAAGCTGAAAAAGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                      <  1:1350621/64‑1 (MQ=255)
                                        cGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCat                                              <  1:542865/56‑1 (MQ=255)
                                        cGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCat                                              <  1:2148107/56‑1 (MQ=255)
                                         tcAACTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCattatt                                          <  1:2025757/58‑1 (MQ=255)
                                             cTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGtt                                   >  1:1127922/1‑62 (MQ=255)
                                             cTGAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGt                                    >  1:2090310/1‑61 (MQ=255)
                                             cTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCtg                          <  1:1295932/71‑1 (MQ=255)
                                               gAAGCTGAAAAACGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTaa                                 >  1:2175123/1‑62 (MQ=255)
                                               gAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCtgat                        >  1:667225/1‑71 (MQ=255)
                                                        aaaCGGCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCtt               >  1:688408/1‑71 (MQ=255)
                                                            ggCGCAGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCt                                      <  1:2080212/44‑1 (MQ=255)
                                                                 aGCGAGAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACgg                              >  1:1031180/1‑47 (MQ=255)
                                                                     aaaGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCctgctg  >  1:144120/1‑71 (MQ=255)
                                                                     aaaGGAGAGCAGGCATGAGTAAGCCATTATTATCTGtt                                   <  1:1069541/38‑1 (MQ=255)
                                                                      |                                                                     
ATGATCTGGCGTCCGCAGGGCTTGCTGCCCATGACGCGCCCGCAACTGAAGCTGAAAAACGGCGCAGCGAAAGGAGAGCAGGCATGAGTCAGCCATTATTATCTGTTAACGGCCTGATGATGCGCTTCGGCGGCCTGCTG  >  minE/2632900‑2633039

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: