Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,847,184 A→G 38.1% Q313Q (CAA→CAG iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,847,1840AG38.1% 56.9 / 73.3 110Q313Q (CAA→CAGiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base A (53/15);  new base G (26/16);  total (79/31)
Fisher's exact test for biased strand distribution p-value = 8.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 120 total aligned reads displayed.

CAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGC  >  minE/1847116‑1847248
                                                                    |                                                                
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:2707557/69‑1 (MQ=255)
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:620458/69‑1 (MQ=255)
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:558057/69‑1 (MQ=255)
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:1705044/69‑1 (MQ=255)
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:3153840/69‑1 (MQ=255)
cAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:1556529/69‑1 (MQ=255)
 aaCGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCaa                                                                  <  1:1549724/68‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTTCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1379946/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1765682/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:3120967/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1668222/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1881844/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1887768/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1515332/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1471977/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:144535/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1244246/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:2733458/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCAACAGCatat                                                             <  1:1100448/69‑1 (MQ=255)
     cgcAAAAACAGCCGCATTCCCTGCGGGAAATAGCTGGCATCACGTAAGACTGGATAATCAACAGCatat                                                             <  1:293795/69‑1 (MQ=255)
       cAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTg                                                           <  1:2279696/69‑1 (MQ=255)
                 tGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGctct                                                 >  1:323879/1‑69 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:285036/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:1022419/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:2838386/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:2784096/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:128503/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:1719556/1‑68 (MQ=255)
                     ttCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTcc                                              >  1:1096396/1‑68 (MQ=255)
                               gAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAg                                   <  1:2158945/69‑1 (MQ=255)
                                 aaTAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGaa                                 >  1:3419583/1‑69 (MQ=255)
                                 aaTAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGaa                                 >  1:3123321/1‑69 (MQ=255)
                                     gCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAg                                                      >  1:1847969/1‑44 (MQ=255)
                                     gCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAg                                                      >  1:1489594/1‑44 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:1791680/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:2834846/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:3301074/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:2280977/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:1422302/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:431848/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:2187337/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:2883163/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:1944949/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:799681/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:318301/1‑69 (MQ=255)
                                           aTGACGTAAGACTGGATAATCACCAACATATTGATAGGCTCTTCCTGGAAGAATAGAACGTCGCTGcc                       >  1:2773870/1‑68 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:2927616/1‑69 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:505127/1‑69 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:3230943/1‑69 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:956246/1‑69 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:360800/1‑69 (MQ=255)
                                              aCGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTg                    >  1:3358850/1‑69 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2896468/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:3430403/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:3296411/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:3062606/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:3353613/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2166947/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1726470/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:964718/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1103318/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1331315/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1348644/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:142142/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1653844/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:673970/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2654984/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1726814/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:1898800/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2630414/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:357888/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2414164/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2144691/1‑67 (MQ=255)
                                                 tAAGACTGGATAATCAACAGCATATTGATAAGGCACTTCCTGGAAGAATAGAACGTCGCTGCCGTGat                  >  1:2290092/1‑67 (MQ=255)
                                                  aaGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                <  1:1620158/69‑1 (MQ=255)
                                                  aaGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                <  1:325826/69‑1 (MQ=255)
                                                    gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                >  1:3216789/1‑67 (MQ=255)
                                                    gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                >  1:227461/1‑67 (MQ=255)
                                                    gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                >  1:2468671/1‑67 (MQ=255)
                                                    gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                >  1:3319226/1‑67 (MQ=255)
                                                    gACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGt                >  1:279393/1‑67 (MQ=255)
                                                       tGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATg            >  1:1713157/1‑68 (MQ=255)
                                                                aCCAACATATTGATAAGTCTCTTCCTGGAAGAATAGAACGTc                             <  1:3429489/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:2380438/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:822506/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:504568/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:3011086/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:2875969/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:2418697/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:2208833/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:215652/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:1955027/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:1857160/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:1790628/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:1672652/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:1627867/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             >  1:1371312/1‑42 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTc                             <  1:1343442/42‑1 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGc  >  1:254329/1‑69 (MQ=255)
                                                                aCCAACATATTGATAAGGCTCTTCCTGGAGAATAGAACGTc                             <  1:1084185/41‑1 (MQ=255)
                                                                    |                                                                
CAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGC  >  minE/1847116‑1847248

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: