Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,451,454 T→A 8.2% L269L (CTT→CTA hdfR → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,451,4540TA 8.2% 200.4 / 8.8 73L269L (CTT→CTAhdfRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (3/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CACTTTCACGGCCGCTTTATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGATGCCGGTAGAA  >  Exported/2451387‑2451517
                                                                   |                                                               
cacTTTCACGGCCGCTTTATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCtttt                                                               >  1:2287162/1‑70 (MQ=255)
 acTTTCACGGCCGCTTTATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTg                                                              >  1:5264047/1‑70 (MQ=255)
       aCGGCCGCTTTATGCCATATGGATGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAAtt                                                        >  1:5115325/1‑70 (MQ=255)
            cGCTTTATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTg                                                  <  1:2457298/71‑1 (MQ=255)
                 tATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTg                                               <  1:4915468/69‑1 (MQ=255)
                      cATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAg                                        >  1:863952/1‑71 (MQ=255)
                       atatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGt                                                   <  1:917573/59‑1 (MQ=255)
                       atatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGt                                                   <  1:656175/59‑1 (MQ=255)
                       atatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGt                                                   <  1:3092972/59‑1 (MQ=255)
                       atatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTg                                               <  1:4011528/63‑1 (MQ=255)
                        tatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGAt                                            <  1:1210011/65‑1 (MQ=255)
                        tatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGAt                                            <  1:3358641/65‑1 (MQ=255)
                        tatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtg                                      <  1:5113653/71‑1 (MQ=255)
                        tatGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtg                                      <  1:3430920/71‑1 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAAtt                                                        >  1:193513/1‑52 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAAtt                                                        >  1:180268/1‑52 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAAtt                                                        >  1:1747555/1‑52 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGc                                                 >  1:5313079/1‑59 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGc                                                 >  1:5648206/1‑59 (MQ=255)
                         atGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtg                                      >  1:2190642/1‑70 (MQ=255)
                          tGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACg                                                    >  1:4076769/1‑55 (MQ=255)
                          tGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTg                                                  <  1:1954079/57‑1 (MQ=255)
                          tGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGc                                                 >  1:3344185/1‑58 (MQ=255)
                          tGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTa                                    >  1:2869907/1‑71 (MQ=255)
                          tGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTa                                    >  1:2075814/1‑71 (MQ=255)
                              tGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGTATGaa                                         <  1:4331349/62‑1 (MQ=255)
                                cAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAAc                                                     >  1:22741/1‑48 (MQ=255)
                                cAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtgt                                     >  1:4266389/1‑64 (MQ=255)
                                cAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTa                                    >  1:2787918/1‑65 (MQ=255)
                                cAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTa                                    >  1:5107518/1‑65 (MQ=255)
                                 aaaaTAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtgt                                     <  1:2134698/63‑1 (MQ=255)
                                 aaaaTAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAatat                             >  1:5060283/1‑71 (MQ=255)
                                 aaaaTAGCGATAAAAATGCGTTGAATCGCGATCTTTTGAAAATTAACGTGCTg                                               <  1:975269/53‑1 (MQ=255)
                                  aaaTAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtgt                                     <  1:3054852/62‑1 (MQ=255)
                                   aaTAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGa                                             >  1:665794/1‑53 (MQ=255)
                                   aaTAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGa                                             >  1:2987512/1‑53 (MQ=255)
                                     tAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTa                                                       <  1:921163/41‑1 (MQ=255)
                                     tAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTa                                                       <  1:4411237/41‑1 (MQ=255)
                                         gATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTaa                                                      <  1:2086788/38‑1 (MQ=255)
                                         gATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTaa                                                      <  1:1923739/38‑1 (MQ=255)
                                         gATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTaa                                                      <  1:5403180/38‑1 (MQ=255)
                                         gATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTaa                                                      <  1:4214594/38‑1 (MQ=255)
                                         gATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTaa                                                      <  1:758246/38‑1 (MQ=255)
                                           tAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAgtg                                      <  1:1376925/52‑1 (MQ=255)
                                            aaaaaTGCGTTGATTCGCGATCTTTTGAAAATTat                                                      <  1:5498115/35‑2 (MQ=255)
                                             aaaaTGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAAtggctggc                  >  1:1484428/1‑70 (MQ=255)
                                             aaaaTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctggc                  <  1:5594853/70‑1 (MQ=255)
                                             aaaaTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctggc                  <  1:1891037/70‑1 (MQ=255)
                                             aaaaTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctggc                  <  1:4973114/70‑1 (MQ=255)
                                               aaTGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAg               <  1:1732527/71‑1 (MQ=255)
                                               aaTGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAg               <  1:1164958/71‑1 (MQ=255)
                                               aaTGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATAAATATGAATGGCTGGCAAg               <  1:1728646/71‑1 (MQ=255)
                                                aTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctg                    >  1:3787375/1‑65 (MQ=255)
                                                aTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctg                    >  1:3977220/1‑65 (MQ=255)
                                                aTGCGTTGATTCGCGATCTATTGAAAATTAACGTGCTGGATGAAGTGTATTAATCTGAAtggctg                    >  1:961039/1‑65 (MQ=255)
                                                  gCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAAtgg                       <  1:1304831/60‑1 (MQ=255)
                                                     ttGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCaa                <  1:5255762/64‑1 (MQ=255)
                                                     ttGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCaa                <  1:1879326/64‑1 (MQ=255)
                                                      tGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAat                               >  1:4523646/1‑48 (MQ=255)
                                                      tGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAat                               >  1:3015435/1‑48 (MQ=255)
                                                      tGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAat                               >  1:2859228/1‑48 (MQ=255)
                                                      tGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGATGCCg        >  1:3999197/1‑71 (MQ=255)
                                                       gATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAg                                        >  1:5026540/1‑38 (MQ=255)
                                                       gATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAata                              >  1:5624068/1‑48 (MQ=255)
                                                       gATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGaa                          >  1:5134135/1‑52 (MQ=255)
                                                       gATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGaa                          >  1:3279930/1‑52 (MQ=255)
                                                        aTTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGa             >  1:938753/1‑64 (MQ=255)
                                                        aTTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGa             >  1:3212467/1‑64 (MQ=255)
                                                         ttCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTaa                                <  1:3083856/44‑1 (MQ=255)
                                                         ttCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAata                              <  1:454889/46‑1 (MQ=255)
                                                           cgcgATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGATGCCGGTAg    >  1:2802642/1‑70 (MQ=255)
                                                            gcgATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGATGCCGGTAGaa  <  1:3850691/71‑1 (MQ=255)
                                                               aTCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTa                                 <  1:5037463/37‑1 (MQ=255)
                                                               aTCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTa                                 <  1:4616324/37‑1 (MQ=255)
                                                                   |                                                               
CACTTTCACGGCCGCTTTATGCCATATGGCTGCAAAATAGCGATAAAAATGCGTTGATTCGCGATCTTTTGAAAATTAACGTGCTGGATGAAGTGTATTAATATGAATGGCTGGCAAGGATGCCGGTAGAA  >  Exported/2451387‑2451517

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