Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,725,308 G→A 45.4% F974F (TTC→TTT purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,725,3080GA45.4% 11.6 / 51.3 66F974F (TTC→TTTpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base G (22/14);  new base A (10/20);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 2.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

GCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATA  >  minE/1725244‑1725370
                                                                |                                                              
gCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATa                                                            <  1:3336221/69‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCTGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:950524/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1077318/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:2186005/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:2182138/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:800937/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:755600/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:725972/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:2489550/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:2711690/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1702903/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1698302/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1686813/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1530584/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1178656/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:618505/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:352031/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:905081/61‑1 (MQ=255)
                       gTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1240615/61‑1 (MQ=255)
                        tttCTGCCCACGAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:1271349/60‑1 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:2413421/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:2216990/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:2255742/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:1975857/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:458087/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:2664189/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:3392123/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:1434697/1‑68 (MQ=255)
                          tCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:940311/1‑68 (MQ=255)
                          tCTGCCCACCAGACAAGCAACGTGGTGAGGCTTTCGCTAAATACAGTCTGCCCGTTGGCGGTAATcac                                   >  1:2909166/1‑68 (MQ=255)
                           cTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTAAATACAGTCTGCCCGTTgg                                             <  1:791996/57‑1 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGGCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc                 >  1:3391653/1‑69 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGt                           >  1:1321517/1‑59 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGaaa                  >  1:620044/1‑68 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc                 >  1:1954678/1‑69 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc                 >  1:2420927/1‑69 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc                 >  1:2230933/1‑69 (MQ=255)
                                           cAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc                 >  1:1708444/1‑69 (MQ=255)
                                                 ggtgCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTg           <  1:2590610/69‑1 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:1735668/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:577046/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:1284111/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:1312969/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:1707812/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:2763208/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:264213/1‑69 (MQ=255)
                                                    gCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGccc        >  1:2614221/1‑69 (MQ=255)
                                                       gCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTa      >  1:1415498/1‑68 (MQ=255)
                                                       gCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTa      >  1:2446596/1‑68 (MQ=255)
                                                       gCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTa      >  1:1825012/1‑68 (MQ=255)
                                                       gCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTa      >  1:1836752/1‑68 (MQ=255)
                                                       gCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTa      >  1:1217625/1‑68 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:812005/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:1118461/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:583151/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:3168886/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:284872/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:2337114/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:2084180/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:1974691/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:1784307/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:152130/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACAt   <  1:1512388/69‑1 (MQ=255)
                                                         tttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTCCCCTACAt   <  1:954872/69‑1 (MQ=255)
                                                          ttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATa  >  1:2389325/1‑69 (MQ=255)
                                                          ttCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATa  >  1:1372787/1‑69 (MQ=255)
                                                                |                                                              
GCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATA  >  minE/1725244‑1725370

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: