Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,475,814 C→T 19.8% R293R (CGC→CGT 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,8140CT19.8% 158.2 / 26.3 105R293R (CGC→CGTnrdBribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base C (26/58);  new base T (7/14);  total (33/72)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 105 total aligned reads displayed.

CTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATAT  >  minE/1475754‑1475879
                                                            |                                                                 
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:1843541/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                          <  1:567001/68‑1 (MQ=255)
                  ggCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATa                                          <  1:412162/68‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                        <  1:3136420/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                        <  1:1320838/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                        <  1:755952/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATaaa                                        <  1:1028597/69‑1 (MQ=255)
                   gCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTATGAATa                                          <  1:3372966/67‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:236776/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:2632634/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:1185060/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:2739810/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:2810064/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:1208362/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:1903535/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:389885/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:632681/69‑1 (MQ=255)
                     aGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGa                                      <  1:785087/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:3238037/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:3399947/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2788111/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:2694006/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2311848/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2509102/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2006339/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:1338844/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:1611132/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2361539/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:137718/69‑1 (MQ=255)
                       cTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACa                                    <  1:2475949/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:1303355/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2066153/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2017231/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:3042603/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2492153/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2860241/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2409953/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:1530382/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:2332742/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGcc                           >  1:1461382/1‑68 (MQ=255)
                                 gaAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCCGAATAAAGACATTCTCTGcc                           >  1:1186594/1‑68 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                <  1:3222564/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                <  1:3428180/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                <  1:2174150/68‑1 (MQ=255)
                                            cGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGaa                <  1:2670709/68‑1 (MQ=255)
                                             tgATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa              <  1:343183/68‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa              <  1:3187083/69‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa              <  1:3135101/69‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa              <  1:2633919/69‑1 (MQ=255)
                                             ggATTATCTGTTCCGCGACGGTTCGATGATTGGGCTGAATAAAGACATTCTCTGCCAGTACGTTGAATa              <  1:2234381/69‑1 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACg                     >  1:273982/1‑59 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt           >  1:1456852/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt           >  1:84037/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt           >  1:1049193/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt           >  1:3133081/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACAt           >  1:2886230/1‑69 (MQ=255)
                                                ttATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTATCTGCCAGTACGTTGAATACAt           >  1:1340160/1‑69 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:740642/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:3254579/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:3270247/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:807894/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:2043673/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:1474872/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:3004880/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:74737/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:1575544/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:2988126/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:2744610/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:1793799/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:179388/69‑1 (MQ=255)
                                                         ccGTGATGGTTCGATGATTGGTCTGAATAAAGACATCCTCTGCCAGTACGTTGAATACATCACCAatat  <  1:1842255/69‑1 (MQ=255)
                                                            |                                                                 
CTATGACCTGTTTGTTCAGGCAGCTCAACAGGAGAAAGACTGGGCGGATTATCTGTTCCGCGACGGTTCGATGATTGGTCTGAATAAAGACATTCTCTGCCAGTACGTTGAATACATCACCAATAT  >  minE/1475754‑1475879

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: