Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,144,422 G→A 32.1% K231K (AAG→AAA argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,4220GA32.1% 77.2 / 42.3 87K231K (AAG→AAAargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (44/15);  new base A (13/15);  total (57/30)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTAA  >  minE/2144356‑2144484
                                                                  |                                                              
ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:1098975/1‑69 (MQ=255)
ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:2956815/1‑69 (MQ=255)
ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:598575/1‑69 (MQ=255)
ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:3047498/1‑69 (MQ=255)
ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:1398090/1‑69 (MQ=255)
ccTCAACTCCAGCGTCAAAATCGTCAACACGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                              >  1:894434/1‑69 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:3177734/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:573555/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:3173620/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:595814/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:3154003/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:2267539/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:421001/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:2271442/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:3284775/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:1859772/1‑67 (MQ=255)
        ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:806151/1‑67 (MQ=255)
        ccAGCGCCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:3301281/1‑67 (MQ=255)
        caaGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                                        >  1:1857319/3‑67 (MQ=255)
           gCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaag                                                   <  1:1364873/69‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAACTCCCGGCagaa                                                    <  1:1432165/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:1711002/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:2351617/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:3032854/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:2819347/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:1817699/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:2680243/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:1796152/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:374876/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:700208/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:377449/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:597003/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:1372882/68‑1 (MQ=255)
           gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGAGGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:620084/68‑1 (MQ=255)
            cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:944363/69‑1 (MQ=255)
            cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:3191516/69‑1 (MQ=255)
            cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:2281204/69‑1 (MQ=255)
                           cccGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCt                                    <  1:1343529/68‑1 (MQ=255)
                           cccGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCt                                    <  1:1432976/68‑1 (MQ=255)
                           cccAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                                    <  1:589081/52‑1 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTcac                                            >  1:964357/3‑58 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:524333/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:2609414/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:796176/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:3287628/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:3353460/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:1373602/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:1761053/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:1541570/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:2280217/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:1222809/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                                    >  1:2745004/3‑66 (MQ=255)
                             caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCAAAGTACGCt                                    >  1:2839413/3‑66 (MQ=255)
                                ttATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAAc                              <  1:1180909/69‑1 (MQ=255)
                                ttATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAAc                              <  1:152795/69‑1 (MQ=255)
                                  aTGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACa                             <  1:901397/68‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:1341315/38‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  >  1:2746420/1‑38 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:1511384/38‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:533230/38‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  >  1:1564080/1‑38 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  >  1:1604603/1‑38 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:1759323/38‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  <  1:2765929/38‑1 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  >  1:2653744/1‑38 (MQ=255)
                                           aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                                  >  1:2087630/1‑38 (MQ=255)
                                                cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc              >  1:343636/1‑69 (MQ=255)
                                                cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc              >  1:1545729/1‑69 (MQ=255)
                                                cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc              >  1:1185618/1‑69 (MQ=255)
                                                cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc              >  1:1056292/1‑69 (MQ=255)
                                                cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc              >  1:2247629/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1057499/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1851119/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1157067/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:2363321/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1465404/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:447007/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1870923/1‑69 (MQ=255)
                                                     aTGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGa         >  1:1739995/1‑69 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAGGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:630423/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:2181941/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:688007/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:2160635/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:1543398/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:3199453/1‑68 (MQ=255)
                                                        agagCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAAc       >  1:928600/1‑68 (MQ=255)
                                                            cGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTaa  <  1:3186727/69‑1 (MQ=255)
                                                                  |                                                              
CCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGCGCTGAACGGTAA  >  minE/2144356‑2144484

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: