Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,655,091 G→A 62.1% A58T (GCC→ACC)  purN → phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,655,0910GA62.1% 38.5 / 50.9 79A58T (GCC→ACC) purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base A (21/28);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

AATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGT  >  minE/1655025‑1655151
                                                                  |                                                            
aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                            <  1:2173565/69‑1 (MQ=255)
aaTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                            <  1:1244578/69‑1 (MQ=255)
 aTAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                            <  1:1511107/68‑1 (MQ=255)
   aaGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgc                                                         >  1:2520467/1‑69 (MQ=255)
    aGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcg                                                        >  1:2692031/1‑69 (MQ=255)
         gACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGttt                                                   <  1:660343/69‑1 (MQ=255)
                    ccTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTg                                         <  1:1840009/68‑1 (MQ=255)
                           cgcgCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                 <  1:2742034/69‑1 (MQ=255)
                           cgcgCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                 <  1:2357024/69‑1 (MQ=255)
                               cccgccAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                             <  1:264059/69‑1 (MQ=255)
                                ccgccAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                             <  1:120670/68‑1 (MQ=255)
                                ccgccAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                             <  1:2926926/68‑1 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAgc                                                         >  1:1099288/1‑38 (MQ=38)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTc                                            >  1:2206733/1‑51 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:305062/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:3357031/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:1901885/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:31523/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:22821/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:34999/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:442839/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:2552039/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:154456/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:2623724/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:1319882/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:2719902/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:1202095/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:82364/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:847030/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCga                           >  1:1047758/1‑67 (MQ=255)
                                  gccAGGCGGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCga                           >  1:549114/1‑67 (MQ=255)
                                  gccAGGCGGGTATAGCAACACATACGCTCATCAcc                                                            >  1:2586273/1‑35 (MQ=25)
                                  gccAGGCCGGTATTGCAACACATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgag                          >  1:3352838/1‑67 (MQ=255)
                                     aGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                        >  1:2163584/1‑68 (MQ=255)
                                      ggCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                      <  1:3011209/69‑1 (MQ=255)
                                      ggCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                      <  1:1837509/69‑1 (MQ=255)
                                                gCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAttg                 >  1:483755/1‑63 (MQ=255)
                                                gCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaa            >  1:1498652/1‑69 (MQ=255)
                                                gCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaa            >  1:1955184/1‑69 (MQ=255)
                                                gCAACGCATACGCTCATCGCCAGAGCGTTTGACAGTCGTGAAGc                                     >  1:1759740/1‑44 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaaa           >  1:3293650/1‑65 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       <  1:1217433/69‑1 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       >  1:2327154/1‑69 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       <  1:2315334/69‑1 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       >  1:477755/1‑69 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       >  1:2057068/1‑69 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       <  1:2108643/69‑1 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       <  1:532996/69‑1 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       >  1:2244263/1‑69 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       <  1:2947627/69‑1 (MQ=255)
                                                     gCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa       >  1:1432672/1‑69 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:1081023/69‑1 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:313438/69‑1 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:3420023/69‑1 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:2492764/69‑1 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:2397815/69‑1 (MQ=255)
                                                      cATACGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGAc      <  1:2772913/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2060546/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:874088/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1133522/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1179281/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:662666/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1182433/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:541871/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1540870/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1599949/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1682964/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2108653/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:215422/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2194747/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2255520/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:3185242/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2291459/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2353991/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:268881/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:2738675/69‑1 (MQ=255)
                                                          cGCTCATCACCAGCGCGTTTGAAAGTCGTGAGGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1257993/69‑1 (MQ=255)
                                                           gCTCATCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:1569337/68‑1 (MQ=255)
                                                               aTCACCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGAGAGTTGATTCATGAAATCGACATGt  <  1:562167/64‑1 (MQ=255)
                                                                  |                                                            
AATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGT  >  minE/1655025‑1655151

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