Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,003 A→C 41.9% K821Q (AAA→CAA) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,0030AC41.9% 1.8 / 53.3 43K821Q (AAA→CAA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base C (8/10);  total (20/23)
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.88e-01

TGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGT  >  Exported/833935‑834071
                                                                    |                                                                    
tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                                                      <  1:2070046/69‑1 (MQ=255)
 gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa                                                                     >  1:4454508/1‑69 (MQ=255)
   aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                                 <  1:3035158/71‑1 (MQ=255)
         cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                           >  1:5187583/1‑71 (MQ=255)
          ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg                                                          <  1:1846389/71‑1 (MQ=255)
          ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                           >  1:3689816/1‑70 (MQ=255)
             tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa                                                       <  1:1481207/71‑1 (MQ=255)
             tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa                                                       <  1:2915590/71‑1 (MQ=255)
             tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                                                         <  1:3322576/69‑1 (MQ=255)
                gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa                                                               >  1:2439283/1‑60 (MQ=255)
                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt                                                    <  1:3401193/70‑1 (MQ=255)
                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                                   <  1:771430/71‑1 (MQ=255)
                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                                   <  1:4377260/71‑1 (MQ=255)
                   tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                                 <  1:78578/55‑1 (MQ=255)
                     cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa                                               <  1:1469837/71‑1 (MQ=255)
                     cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa                                                 >  1:4287315/1‑69 (MQ=255)
                     cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa                                                 >  1:3522696/1‑69 (MQ=255)
                      gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt                                                      <  1:4112380/63‑1 (MQ=255)
                       tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                                                      <  1:2210751/46‑1 (MQ=255)
                       tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                                 <  1:5515961/51‑1 (MQ=255)
                        ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg                                             >  1:4407565/1‑69 (MQ=255)
                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                           >  1:5665279/1‑55 (MQ=255)
                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                           >  1:5369519/1‑55 (MQ=255)
                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta                                           >  1:4662936/1‑68 (MQ=255)
                         tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa                                           >  1:3288003/1‑71 (MQ=255)
                         tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa                                           >  1:5064259/1‑71 (MQ=255)
                               tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                     <  1:5645485/70‑3 (MQ=255)
                                  gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                                   >  1:4232274/1‑70 (MQ=255)
                                   ctttCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaa                                    >  1:2376203/4‑68 (MQ=255)
                                               gTCGATCCGGCGGGGGAATCACGCAGCGACTGGGAAATCTACaa                                                <  1:2095785/44‑1 (MQ=255)
                                               gTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgt                     <  1:593408/71‑1 (MQ=255)
                                                 cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGaaa                                >  1:3149122/1‑58 (MQ=255)
                                                  gATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGtgtg                  <  1:4110679/71‑1 (MQ=255)
                                                   aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGc                 <  1:1316784/71‑1 (MQ=255)
                                                   aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc                                             <  1:3609045/43‑1 (MQ=255)
                                                   aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt                   <  1:2696619/69‑1 (MQ=255)
                                                   aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                 >  1:217339/1‑71 (MQ=255)
                                                    tCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                                <  1:2132342/39‑1 (MQ=16)
                                                    tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                >  1:1893072/1‑71 (MQ=255)
                                                    tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                >  1:1222129/1‑71 (MQ=255)
                                                     ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt               >  1:3690035/1‑71 (MQ=255)
                                                      cGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt               <  1:5064572/70‑1 (MQ=255)
                                                        gCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc                                            <  1:349508/39‑1 (MQ=255)
                                                                  cGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTggg   >  1:4726129/1‑70 (MQ=255)
                                                                   gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGa                        <  1:5235413/48‑1 (MQ=255)
                                                                   gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGt  >  1:2200177/1‑70 (MQ=255)
                                                                    |                                                                    
TGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGT  >  Exported/833935‑834071

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: