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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,4010AG48.3% ‑3.1 / 50.9 58K224K (AAA→AAGargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (21/9);  new base G (13/15);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACA  >  minE/2144347‑2144457
                                                      |                                                        
tCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGagag                                            >  1:18553/1‑69 (MQ=255)
tCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGagag                                            >  1:1857933/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                   >  1:1398090/1‑69 (MQ=255)
         ccTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                   >  1:598575/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)
         ccTCAACTCCAGCGTCAAAATCGTCAACACGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGAt                                   >  1:894434/1‑69 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:3154003/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:3284775/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:421001/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:2271442/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:2267539/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:573555/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:1859772/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:3173620/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:3177734/1‑67 (MQ=255)
                 ccAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGc                             >  1:595814/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:597003/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:3032854/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:1372882/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:1796152/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:2819347/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:374876/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:1817699/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:2680243/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:377449/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:2351617/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:700208/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCagaa                         <  1:1711002/68‑1 (MQ=255)
                    gCGTCAAAATCGTCAACCCAATTATGGGCGTGAAGTTTTGGGAGGAGAGCGTGAAAATCCCGGCagaa                         <  1:620084/68‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:2281204/69‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:944363/69‑1 (MQ=255)
                     cGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:3191516/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:3353460/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:524333/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:796176/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:2745004/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:3287628/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:1222809/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:1373602/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:2609414/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:2280217/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:1541570/3‑66 (MQ=255)
                                      caaTTATGGGCGTGAAGTTTTGGGATGAGAGCGTGAAAATCCCGGCAGAAGAAGTCACAGTACGCt         >  1:1761053/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:1759323/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       >  1:1604603/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       >  1:1564080/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       >  1:2087630/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:533230/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:1511384/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       >  1:2653744/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:1341315/38‑1 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       >  1:2746420/1‑38 (MQ=255)
                                                    aaaTTCTGGGATGAGAGCGTGAAGATCCCGGCagaaga                       <  1:2765929/38‑1 (MQ=255)
                                                      |                                                        
TCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGATGAGAGCGTGAAGATCCCGGCAGAAGAAGTCACAGTACGCTTTGAACA  >  minE/2144347‑2144457

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: