Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,514,054 G→A 34.5% L479L (CTC→CTT purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,514,0540GA34.5% 100.9 / 74.2 130L479L (CTC→CTTpurFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/65);  new base A (28/17);  total (48/82)
Fisher's exact test for biased strand distribution p-value = 2.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

TTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAA  >  minE/1513990‑1514122
                                                                |                                                                    
tttCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                                                  >  1:3170/1‑69 (MQ=255)
tttCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                                                  >  1:3197977/1‑69 (MQ=255)
tttCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                                                  >  1:1395409/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCg                                                             >  1:1830335/1‑68 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCg                                                             >  1:3207894/1‑68 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:3319767/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:617607/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:54423/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:1925691/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2057438/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:3340272/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2184654/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:1028691/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2269804/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2363571/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2987142/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:2907528/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:1189688/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:1462073/1‑69 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:900879/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:2321138/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:1250265/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:540564/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:268146/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:1505657/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:1577571/68‑1 (MQ=255)
      gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:839451/68‑1 (MQ=255)
      gaTTTTCCACTTCATTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                            >  1:3087875/1‑69 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:1982426/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:1519833/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:1188625/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:1013989/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:2090590/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:1336135/68‑1 (MQ=255)
                    ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                               <  1:573117/68‑1 (MQ=255)
                         aCGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                             <  1:1574104/49‑1 (MQ=255)
                         aCGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATc                                         >  1:1066177/1‑69 (MQ=255)
                                aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:916036/68‑1 (MQ=255)
                                aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:2150684/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGTTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1886549/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:2482994/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1098826/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1146916/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1327528/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:2384956/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:2580171/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:293505/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:3161153/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:333359/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             >  1:125876/1‑68 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:2093022/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:3353418/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1942951/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1933605/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1885526/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:850862/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1424343/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:171894/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:1576849/68‑1 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAg                            >  1:494987/1‑69 (MQ=255)
                                      ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAg                            >  1:1439583/1‑69 (MQ=255)
                                       tGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             <  1:2387347/67‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:2859624/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:944671/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:929222/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:1571749/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:724819/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:1832379/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:2254192/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:2479629/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:2904902/68‑1 (MQ=255)
                                          cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:2838237/68‑1 (MQ=255)
                                              atcatTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGccat                     >  1:652257/1‑66 (MQ=255)
                                              atcatTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCg                      >  1:473672/1‑67 (MQ=255)
                                              atcatTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGt                     >  1:2073931/1‑68 (MQ=255)
                                              atcatTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGt                     >  1:835625/1‑68 (MQ=255)
                                              atcatTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGt                     >  1:1402661/1‑68 (MQ=255)
                                                     cgtaaagtaTCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                         <  1:3275568/57‑1 (MQ=255)
                                                       taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg            >  1:2643284/1‑68 (MQ=255)
                                                       taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg            >  1:1948669/1‑68 (MQ=255)
                                                       taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg            >  1:297911/1‑68 (MQ=255)
                                                       taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg            >  1:1664484/1‑68 (MQ=255)
                                                       taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg            >  1:2726290/1‑68 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:2239454/44‑1 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:2820401/44‑1 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:44434/44‑1 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:1830472/44‑1 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAg          <  1:618897/69‑1 (MQ=255)
                                                        aacgtaTCGAGGAAATCGAGGTAGACCTGATCAACATCTTTGGt                                   <  1:1359174/44‑1 (MQ=255)
                                                         acgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                                   <  1:1698700/43‑1 (MQ=255)
                                                               cGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaa   >  1:2387473/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             >  1:2842164/1‑42 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             >  1:2567167/1‑42 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                             >  1:3150878/1‑42 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:369066/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:2089473/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:2808918/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:3242281/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:3087165/1‑69 (MQ=255)
                                                                gAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCaaa  >  1:1126486/1‑69 (MQ=255)
                                                                |                                                                    
TTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAA  >  minE/1513990‑1514122

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: