Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,514,048 C→A 40.1% T481T (ACG→ACT purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,514,0480CA40.1% 81.0 / 97.1 112T481T (ACG→ACTpurFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/55);  new base A (29/16);  total (41/71)
Fisher's exact test for biased strand distribution p-value = 8.61e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 139 total aligned reads displayed.

ATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAG  >  minE/1513985‑1514114
                                                               |                                                                  
aTGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaacgta                                                                 >  1:2560389/1‑67 (MQ=255)
aTGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaacgta                                                                 >  1:909185/1‑67 (MQ=255)
aTGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaacgta                                                                 >  1:1719758/1‑67 (MQ=255)
aTGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaacgta                                                                 >  1:3402256/1‑67 (MQ=255)
aTGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaacgta                                                                 >  1:1027181/1‑67 (MQ=255)
     tttCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                                          >  1:3170/1‑69 (MQ=255)
     tttCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                                          >  1:1395409/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCGTTACGTAAAGTATCAAGGAAATCGa                                                    >  1:3094633/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATc                                                               >  1:1082568/1‑58 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCg                                                     >  1:3207894/1‑68 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:2987142/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:2907528/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:265758/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:3319767/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:2269804/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:54423/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:2203455/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:617607/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1925691/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1847456/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1819918/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:3233560/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1028691/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1413974/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1372049/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:132679/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1462073/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCGTCATTACGTAAAGTATCAAGGAAATCGa                                                    >  1:1088646/1‑69 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:900879/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:1505657/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:839451/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:268146/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:1250265/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:1577571/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:2321138/68‑1 (MQ=255)
           gaTTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:540564/68‑1 (MQ=255)
                         ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                       <  1:3323330/68‑1 (MQ=255)
                         ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                       <  1:1188625/68‑1 (MQ=255)
                         ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                       <  1:1013989/68‑1 (MQ=255)
                         ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                       <  1:156667/68‑1 (MQ=255)
                         ttCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATc                                       <  1:1982426/68‑1 (MQ=255)
                              aCGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                                     <  1:1574104/49‑1 (MQ=255)
                                     aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:2089950/68‑1 (MQ=255)
                                     aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:916036/68‑1 (MQ=255)
                                     aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:573024/68‑1 (MQ=255)
                                     aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:2150684/68‑1 (MQ=255)
                                     aCTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:20409/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGTTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1886549/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1327528/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2093022/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1576849/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1098826/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:3161153/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1146916/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:850862/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:146491/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:293505/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:291716/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:333359/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2820825/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:3353418/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:274679/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1885526/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     >  1:125876/1‑68 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2580171/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1424343/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2482994/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:240526/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2384956/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:2364429/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1933605/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTa                     <  1:1942951/68‑1 (MQ=255)
                                           ttGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAg                    >  1:494987/1‑69 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:724819/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:3348001/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:827700/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:929222/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:944671/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:2867311/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:123968/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:1571749/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:1832379/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:214803/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:2309023/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:2479629/68‑1 (MQ=255)
                                               cgtcgtCATTACGTAAAGTATCAAGGAAATCGAGGTAGCCCTGATCAACGTCTTTGGTGACGTAGACg                 <  1:2859624/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)
                                                            taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg    >  1:2643284/1‑68 (MQ=255)
                                                            taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg    >  1:1948669/1‑68 (MQ=255)
                                                            taacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCg    >  1:1664484/1‑68 (MQ=255)
                                                             aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:2239454/44‑1 (MQ=255)
                                                             aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:2820401/44‑1 (MQ=255)
                                                             aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:1830472/44‑1 (MQ=255)
                                                             aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:1697525/44‑1 (MQ=255)
                                                             aacgtaTCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt                           <  1:3034194/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)
                                                               |                                                                  
ATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAG  >  minE/1513985‑1514114

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: