Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,475,817 C→T 19.5% D294D (GAC→GAT nrdB → ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,475,8170CT19.5% 163.1 / 26.1 107D294D (GAC→GATnrdBribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base C (26/60);  new base T (7/14);  total (33/74)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 107 total aligned reads displayed.

CTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTA  >  minE/1475754‑1475884
                                                               |                                                                   
cTATGACGTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTc                                                                >  1:142934/1‑69 (MQ=255)
cTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTc                                                                >  1:2381828/1‑69 (MQ=255)
cTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTc                                                                >  1:2357149/1‑69 (MQ=255)
  aTGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCg                                                               >  1:574443/1‑68 (MQ=255)
  aTGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCg                                                               >  1:2971407/1‑68 (MQ=255)
  aTGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCg                                                               >  1:1600461/1‑68 (MQ=255)
            tgttCAGTCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACgg                                                                   <  1:2553268/54‑1 (MQ=255)
            tgttCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACgg                                                                   <  1:2505855/54‑1 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTg                                                   >  1:2617094/2‑67 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTg                                                   >  1:1469198/2‑67 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTGa                                                  >  1:87514/2‑68 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTGa                                                  >  1:2693165/2‑68 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTGa                                                  >  1:2272398/2‑68 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTGa                                                  >  1:1379426/2‑68 (MQ=255)
               ccAGGCGGCACAGCAAGAGAAAGACTGGGCGGATTATCTGTTCCGTGATGGTTCGATGATTGGTCTGa                                                  >  1:1126693/2‑68 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGTATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:1511124/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGTTCTGAATa                                               <  1:1859432/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:82896/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:783232/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:2617681/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:567001/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                               <  1:1843541/68‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                             <  1:1028597/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                             <  1:1320838/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                             <  1:3136420/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                             <  1:755952/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTATGAATa                                               <  1:3372966/67‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:1208362/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:389885/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:1185060/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:632681/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:2632634/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:2810064/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:785087/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:2739810/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGACTGGTCTGAATAAAGa                                           <  1:1079293/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCCACGGTTCGATGATTGGTCTGAATAAAGa                                           <  1:3049423/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATTAAGACa                                         <  1:3793/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2694006/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2509102/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2788111/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:3238037/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:3399947/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:356543/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:49979/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:681948/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2475949/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:137718/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2361539/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2006339/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2338374/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2311848/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:1338844/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:1611132/69‑1 (MQ=255)
                        tCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                         <  1:2297892/68‑1 (MQ=255)
                         cAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACAtt                                       <  1:1904444/69‑1 (MQ=255)
                                agaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCgatga                                                           <  1:2807167/42‑1 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:3042603/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2066153/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2017231/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2860241/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:3426843/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2332742/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:1530382/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:1461382/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:1303355/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2409953/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                                >  1:2492153/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCCGAATAAAGACATTCTCTGcc                                >  1:1186594/1‑68 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:1536313/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:3222564/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:2781323/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:2670709/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:2174150/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                     <  1:3428180/68‑1 (MQ=255)
                                             tgATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa                   <  1:343183/68‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa                   <  1:2633919/69‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa                   <  1:3187083/69‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGGCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa                   <  1:2234381/69‑1 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt                >  1:2886230/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt                >  1:1049193/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt                >  1:3133081/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt                >  1:84037/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTATCTGCCAGTACGTTGAATACAt                >  1:1340160/1‑69 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:179388/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:1474872/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:74737/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:740642/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:807894/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:1575544/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:2744610/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:1793799/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:1842255/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:3270247/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:3254579/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:2988126/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:2043673/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat       <  1:3004880/69‑1 (MQ=255)
                                                              aCGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTa  <  1:112820/69‑1 (MQ=255)
                                                              aCGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTa  <  1:30951/69‑1 (MQ=255)
                                                               |                                                                   
CTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATATCCGTA  >  minE/1475754‑1475884

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: