Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,243,565 G→A 16.5% I11I (ATC→ATT nfi ← endonuclease V

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,243,5650GA16.5% 215.0 / 23.9 126I11I (ATC→ATTnfiendonuclease V
Reads supporting (aligned to +/- strand):  ref base G (42/63);  new base A (9/12);  total (51/75)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 126 total aligned reads displayed.

TCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTC  >  minE/2243499‑2243632
                                                                  |                                                                   
tCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGAt                                                                   >  1:3364348/1‑69 (MQ=255)
tCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGAt                                                                   >  1:3290127/1‑69 (MQ=255)
tCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGAt                                                                   >  1:2428563/1‑69 (MQ=255)
tCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGAt                                                                   >  1:2459005/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGt                                                      <  1:2902817/68‑1 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:2271069/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:2307592/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:1545546/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:2417916/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:2819574/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:138853/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:2829143/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:3265987/1‑69 (MQ=255)
              ccGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:510591/1‑69 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:324991/1‑64 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:3266198/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:609763/1‑64 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:979724/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:110769/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:1463857/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:1409984/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:1314171/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     <  1:1293190/64‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:159839/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:914553/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:2478063/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:424656/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:3410566/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:2811179/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:972227/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:1942086/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:3130994/69‑1 (MQ=255)
                   ggATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                <  1:3083420/69‑1 (MQ=255)
                    gATCTTTATCGAGTCGATCCTCGCGGATCACAGATGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:855484/1‑68 (MQ=255)
                    gATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:2871000/1‑68 (MQ=255)
                    gATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:1835795/1‑68 (MQ=255)
                    gATCTTTATAGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:3003221/1‑68 (MQ=255)
                    gATCTTTATAGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:53641/1‑68 (MQ=255)
                    gATCTTTATAGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:3023782/1‑68 (MQ=255)
                    gATCTTTATAGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:1735970/1‑68 (MQ=255)
                    gATCTTTATAGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCAATTTGTTGAGCGCGTAATGAc                                                >  1:1753097/1‑68 (MQ=255)
                      tCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGt                                                      >  1:794892/1‑60 (MQ=255)
                      tCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTa                                                     >  1:875688/1‑61 (MQ=255)
                                gTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATa                                    <  1:1186905/68‑1 (MQ=255)
                                 tCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAAt                                  <  1:2081211/69‑1 (MQ=255)
                                 tCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAAt                                  <  1:105982/69‑1 (MQ=255)
                                 tCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAAt                                  <  1:2894545/69‑1 (MQ=255)
                                 tCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAAt                                  <  1:2316584/69‑1 (MQ=255)
                                 tCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAAt                                  <  1:2211812/69‑1 (MQ=255)
                                     tCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCAc                               <  1:2757582/68‑1 (MQ=255)
                                      ccTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTc                             >  1:56901/1‑69 (MQ=255)
                                      ccTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTc                             >  1:89061/1‑69 (MQ=255)
                                      ccTCGCGGATCACAGAAGAAGCCAGTTCGAATTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTc                             >  1:2096840/1‑69 (MQ=255)
                                            ggATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGc                       <  1:2931518/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:3057621/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:2885057/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:23067/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:3001440/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:1115209/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:928983/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGt                    <  1:2191773/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTACTTAGCGGt                    <  1:200006/69‑1 (MQ=255)
                                               tCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCAATAATCACTCCTTAGCGGt                    <  1:43251/69‑1 (MQ=255)
                                                 acaGAAGAAGCTAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1976963/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATGACTCCTTAGCGGTGa                  <  1:1940395/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:3128434/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:321101/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1278856/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:3289899/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1143108/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:3036347/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1691130/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2557030/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1698529/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1592114/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:247652/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2492356/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:3011533/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2719264/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2946636/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2862432/69‑1 (MQ=255)
                                                 acaGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:1337562/69‑1 (MQ=255)
                                                  caGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGa                  <  1:2863449/68‑1 (MQ=255)
                                                   agaagaAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATa                <  1:2068641/69‑1 (MQ=255)
                                                        aagCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGt            <  1:1347964/68‑1 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:290780/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:2236580/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:157429/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:1535767/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:1484187/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:2461219/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:2625556/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:1181301/1‑42 (MQ=255)
                                                         agCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCAt                                     >  1:298276/1‑42 (MQ=255)
                                                                tCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATa                <  1:2744917/56‑1 (MQ=255)
                                                                 cGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTc  <  1:567072/69‑1 (MQ=255)
                                                                 cGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTc  <  1:2136416/69‑1 (MQ=255)
                                                                 cGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTc  <  1:2348264/69‑1 (MQ=255)
                                                                 cGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTc  <  1:2775154/69‑1 (MQ=255)
                                                                 cGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTc  <  1:2823230/69‑1 (MQ=255)
                                                                  |                                                                   
TCCGGCGATCAGATCCGGTGGATCTTTATCGAGTCGATCCTCGCGGATCACAGAAGAAGCCAGTTCGATTTGTTGAGCGCGTAATGACGCGAGATCCATAATCACTCCTTAGCGGTGATACTGTTCAGACAGTC  >  minE/2243499‑2243632

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