Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 649,279 G→A 48.4% A653A (GCG→GCA pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE649,2790GA48.4% ‑0.1 / 53.3 62A653A (GCG→GCApepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base A (14/16);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGAT  >  minE/649212‑649342
                                                                   |                                                               
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:631093/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:61049/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3365004/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3184040/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:870526/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2833225/1‑69 (MQ=255)
tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:162298/1‑69 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:717633/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1889601/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2159660/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2202309/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2270281/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2418521/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2684704/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2699246/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2773834/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1039972/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3146874/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3196030/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3233232/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3425462/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:410802/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:686701/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1640649/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1679176/1‑67 (MQ=255)
  atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:167989/1‑67 (MQ=255)
        ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        <  1:1787745/69‑1 (MQ=255)
        ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1256431/1‑69 (MQ=255)
        ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:3089956/1‑69 (MQ=255)
        ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1384253/1‑69 (MQ=255)
        ccGATTGCTATTGCCGAAGTACGCGAAGCACTCAATCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1746206/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCCCTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2375520/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGTCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:777253/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:1439487/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2977839/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2860771/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2843449/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:1914475/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1687972/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1081225/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2371486/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2697056/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1731785/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2700809/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1939681/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:511292/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2929220/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:3314504/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1781457/1‑69 (MQ=255)
                    gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGAGACTGAACTGGCGGATGAGCt                                                         >  1:2751336/1‑56 (MQ=255)
                    gCCGAAGTACGCGAAGCACTAACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1170248/1‑69 (MQ=255)
                       gAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAc                                         <  1:1209284/69‑1 (MQ=255)
                          gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                      <  1:1241995/69‑1 (MQ=255)
                          gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                      <  1:72064/69‑1 (MQ=255)
                             cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                   >  1:2025204/1‑69 (MQ=255)
                             cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                   >  1:2768018/1‑69 (MQ=255)
                                gAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTAc                                <  1:3148309/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:700612/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2572033/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:1000810/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:411763/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:3156675/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:529008/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:3154873/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2879092/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2122453/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:742896/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2822558/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:27993/69‑1 (MQ=255)
                                   gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2684228/69‑1 (MQ=255)
                                    cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2983182/68‑1 (MQ=255)
                                    cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:1350446/68‑1 (MQ=255)
                                    cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2602838/68‑1 (MQ=255)
                                           tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGTTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1477375/69‑1 (MQ=255)
                                           tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1304859/69‑1 (MQ=255)
                                           tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1753305/69‑1 (MQ=255)
                                           tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:2162719/69‑1 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     >  1:1364016/1‑68 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:890897/1‑69 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:954750/1‑69 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:3324481/1‑69 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:2700229/1‑69 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:2444096/1‑69 (MQ=255)
                                            cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:1752492/1‑69 (MQ=255)
                                                              cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2892202/1‑69 (MQ=255)
                                                              cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2880047/1‑69 (MQ=255)
                                                              cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2843020/1‑69 (MQ=255)
                                                              cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGa   <  1:1300499/68‑1 (MQ=255)
                                                                   |                                                               
TCATCGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGAT  >  minE/649212‑649342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: