Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,494 G→C 28.6% T984T (ACG→ACC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,4940GC28.6% 65.2 / 33.6 49T984T (ACG→ACCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (23/12);  new base C (11/3);  total (34/15)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAA  >  Exported/834425‑834561
                                                                     |                                                                   
ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACg                                                                     >  1:1766373/1‑70 (MQ=255)
 gAAACCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCa                                                                   >  1:1453460/1‑71 (MQ=255)
 gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGc                                                                    >  1:976663/1‑70 (MQ=255)
 gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCa                                                                   >  1:139084/1‑71 (MQ=255)
 gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCa                                                                   <  1:1737948/71‑1 (MQ=255)
    aCCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCATCt                                                                >  1:836391/1‑71 (MQ=255)
     ccAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTg                                                               <  1:4765093/71‑1 (MQ=255)
      cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTg                                                               <  1:4758660/70‑1 (MQ=255)
        aCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                                            <  1:1617288/71‑1 (MQ=255)
         cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                                            >  1:1671944/1‑70 (MQ=255)
         cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                                            >  1:1951297/1‑70 (MQ=255)
           gTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTg                                                         <  1:1907621/71‑1 (MQ=255)
                gTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATa                                                     >  1:2461451/1‑70 (MQ=255)
                  aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTg                                                         >  1:331218/1‑64 (MQ=255)
                  aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATa                                                     >  1:5537256/1‑68 (MQ=255)
                   gCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAaag                                                  >  1:3476574/1‑70 (MQ=255)
                         aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATaa                                                    >  1:3609774/1‑62 (MQ=255)
                         aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAagaaga                                              >  1:5674681/1‑68 (MQ=255)
                         aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                             >  1:1786923/1‑69 (MQ=255)
                         aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                             >  1:2851789/1‑69 (MQ=255)
                          aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgag                                          <  1:5475302/71‑1 (MQ=255)
                                   cTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCAGGCGCTGAATAAAGAAGACGAGAAGATCCGc                                 >  1:4356105/1‑71 (MQ=255)
                                        cTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcg                            >  1:2001859/1‑71 (MQ=255)
                                          gAGCGAATTTACCGGTCGTGACCATACGCATCTGGGGCTGAATAAAGAAGACGAGAAGATCCGCTTCcgc                           >  1:3195853/1‑70 (MQ=255)
                                           aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa                         <  1:1801173/71‑1 (MQ=255)
                                           aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa                         >  1:1934889/1‑71 (MQ=255)
                                           aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa                         <  1:2626879/71‑1 (MQ=255)
                                               aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat                      >  1:3871637/1‑70 (MQ=255)
                                               aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt                     <  1:3332867/71‑1 (MQ=255)
                                                aTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt                     >  1:3796856/1‑70 (MQ=255)
                                                   taacGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCa                   >  1:4472993/4‑69 (MQ=255)
                                                   tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgaga                                         >  1:2373374/1‑47 (MQ=255)
                                                   tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCa                   >  1:1068834/1‑69 (MQ=255)
                                                      cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGata                       <  1:266052/62‑1 (MQ=255)
                                                      cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg               >  1:5275372/1‑69 (MQ=255)
                                                      cGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg               >  1:571306/1‑69 (MQ=255)
                                                       ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcaca             <  1:1678312/71‑1 (MQ=255)
                                                          cGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGata                       >  1:2690224/1‑58 (MQ=255)
                                                              aacATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcg       >  1:5178234/3‑70 (MQ=255)
                                                               acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGc           >  1:2943682/2‑65 (MQ=255)
                                                               acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGc           >  1:650790/2‑65 (MQ=255)
                                                               acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt      >  1:1830055/2‑69 (MQ=255)
                                                               acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt      >  1:4818750/2‑69 (MQ=255)
                                                               acATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcgt      >  1:363581/2‑69 (MQ=255)
                                                                cATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaa    <  1:5646884/71‑1 (MQ=255)
                                                                cATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaa    <  1:4841285/71‑1 (MQ=255)
                                                                cATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaa    <  1:4857332/71‑1 (MQ=255)
                                                                 aTACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCg                                  <  1:2931192/40‑1 (MQ=255)
                                                                  tACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTaaaa  >  1:4450702/1‑71 (MQ=255)
                                                                    cGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgc        >  1:1926493/1‑63 (MQ=255)
                                                                     |                                                                   
GGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAA  >  Exported/834425‑834561

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: