Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,144,404 C→T 41.1% F225F (TTC→TTT argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,4040CT41.1% 26.9 / 46.6 68F225F (TTC→TTTargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/14);  new base T (13/15);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGC  >  minE/2144347‑2144472
                                                         |                                                                    
tCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGagag                                                           >  1:1857933/1‑69 (MQ=255)
tCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGagag                                                           >  1:18553/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:1098975/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:2956815/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:3047498/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:598575/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:1398090/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACACGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                                  >  1:894434/1‑69 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:573555/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:1859772/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:806151/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:3284775/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:2267539/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:2271442/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:421001/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:3177734/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:595814/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:3173620/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                                            >  1:3154003/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:1372882/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:2351617/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:2680243/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:1796152/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:374876/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:377449/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:1711002/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:597003/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:700208/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGAGGAGAGCGTGAAAATCCCGGCagaa                                        <  1:620084/68‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:3191516/69‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:2281204/69‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:944363/69‑1 (MQ=255)
                                    cccGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCt                        <  1:1432976/68‑1 (MQ=255)
                                    cccGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCt                        <  1:1343529/68‑1 (MQ=255)
                                    cccAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                                        <  1:589081/52‑1 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTcac                                >  1:964357/3‑58 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:796176/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:1541570/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:524333/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:1373602/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:1761053/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:3353460/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:3287628/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:1222809/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:2280217/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt                        >  1:2609414/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:2765929/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:1341315/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:1511384/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      >  1:1564080/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:533230/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      >  1:1604603/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      <  1:1759323/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      >  1:2087630/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      >  1:2653744/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                                      >  1:2746420/1‑38 (MQ=255)
                                                         cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc  >  1:1545729/1‑69 (MQ=255)
                                                         cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc  >  1:343636/1‑69 (MQ=255)
                                                         cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc  >  1:2247629/1‑69 (MQ=255)
                                                         cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc  >  1:1185618/1‑69 (MQ=255)
                                                         cTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGgc  >  1:1056292/1‑69 (MQ=255)
                                                         |                                                                    
TCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACAAGGTCATCCGGTGGC  >  minE/2144347‑2144472

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: