Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,443 A→T 24.2% V967V (GTA→GTT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,4430AT24.2% 113.2 / 37.7 62V967V (GTA→GTTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (29/18);  new base T (10/5);  total (39/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 = 8.48e-01

TGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATA  >  Exported/834376‑834510
                                                                   |                                                                   
tGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTTAGGTAGcc                                                                  >  1:3078595/1‑71 (MQ=255)
 gCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCg                                                                 >  1:524524/1‑71 (MQ=255)
  cTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCAGAAACCAACGGTCAGGTAGCCGt                                                                >  1:2728700/1‑71 (MQ=255)
   ttgacacCGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGt                                                                >  1:3927219/4‑70 (MQ=255)
     gacacCGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGtca                                                              >  1:2513408/2‑68 (MQ=255)
     aaCACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAg                                                                    >  1:2180490/1‑64 (MQ=255)
     aaCACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTg                                                               >  1:1513117/1‑69 (MQ=255)
     aaCACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGa                                                              >  1:1088208/1‑70 (MQ=255)
     aaCACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGa                                                              >  1:1921058/1‑70 (MQ=255)
        acCGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGc                                                          <  1:5314949/71‑1 (MQ=255)
        acCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCACCGGTCAGGTAGCCGTGAAAGc                                                          <  1:5325265/71‑1 (MQ=255)
        acCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGc                                                          <  1:2065873/71‑1 (MQ=255)
         ccGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGcg                                                         <  1:4829589/71‑2 (MQ=255)
         ccGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGcg                                                         <  1:2907865/71‑2 (MQ=255)
         ccGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGcc                                                         <  1:2940260/71‑1 (MQ=255)
          cGCCATTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGcg                                                         <  1:2075062/70‑2 (MQ=255)
            ccaTTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTgg                                                      >  1:4249939/3‑71 (MQ=255)
            ccaTTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTgg                                                      >  1:1319085/3‑71 (MQ=255)
             caTTGACGCTTCAGAAGTGATTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTggc                                                     <  1:2837042/70‑2 (MQ=255)
                tGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTg                                                       >  1:5306018/1‑66 (MQ=255)
                 gATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTg                                                       >  1:5443950/1‑65 (MQ=255)
                 gATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGc                                                 >  1:3230775/1‑71 (MQ=255)
                 gATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGc                                                 >  1:4143126/1‑71 (MQ=255)
                         agagATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGc                                          >  1:1048035/1‑70 (MQ=255)
                         agagATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCg                                         <  1:5142443/71‑1 (MQ=255)
                             atgatCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGaa                                       <  1:5637538/69‑1 (MQ=255)
                               gatTCTGGCACTGGCCCCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTggca                                                    >  1:4797438/1‑52 (MQ=255)
                               gatTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTggca                                                    >  1:1999719/1‑52 (MQ=255)
                               gatTCTGGCACTGGCACCGGAAACCAACGGTCATGTTGCAGTCAAAGCGTgg                                                      >  1:3563902/1‑52 (MQ=255)
                               gatCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGaa                                       >  1:5447706/1‑67 (MQ=255)
                                         cTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCt                                                        <  1:166028/40‑1 (MQ=255)
                                         cTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCt                                                        <  1:1479846/40‑1 (MQ=255)
                                           gtcACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGa                                            >  1:1986110/3‑50 (MQ=255)
                                           ggCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGa                                            >  1:2867190/1‑50 (MQ=255)
                                              aCCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                           <  1:2993730/64‑1 (MQ=255)
                                                cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTggg                                                     >  1:920425/1‑36 (MQ=255)
                                                cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTggg                                                     >  1:1487802/1‑36 (MQ=255)
                                                cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTggg                                                     >  1:1033794/1‑36 (MQ=255)
                                                 ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACg                  >  1:1766373/1‑70 (MQ=255)
                                                  gAAACCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATAc                   >  1:3798077/1‑68 (MQ=255)
                                                  gAAACCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCa                >  1:1453460/1‑71 (MQ=255)
                                                  gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGc                 >  1:976663/1‑70 (MQ=255)
                                                  gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCa                <  1:1737948/71‑1 (MQ=255)
                                                  gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCa                >  1:139084/1‑71 (MQ=255)
                                                   aaaCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                           <  1:4408041/59‑1 (MQ=255)
                                                   aaaCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                           >  1:3198639/1‑59 (MQ=255)
                                                   aaaCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTg                          >  1:518618/1‑60 (MQ=255)
                                                     aCCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCATCt             >  1:836391/1‑71 (MQ=255)
                                                     aCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGa                         <  1:3147622/59‑1 (MQ=255)
                                                     aCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGa                         <  1:3742897/59‑1 (MQ=255)
                                                      ccAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCa                      <  1:4148986/61‑1 (MQ=255)
                                                      ccAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTg            <  1:4765093/71‑1 (MQ=255)
                                                       cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                           >  1:2733766/1‑55 (MQ=255)
                                                       cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTg            <  1:4758660/70‑1 (MQ=255)
                                                         aCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAAtt                                     >  1:3574503/1‑43 (MQ=255)
                                                         aCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg         <  1:1617288/71‑1 (MQ=255)
                                                          cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg         >  1:1671944/1‑70 (MQ=255)
                                                          cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg         >  1:1951297/1‑70 (MQ=255)
                                                           ggTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCa                      <  1:4297539/56‑1 (MQ=255)
                                                            gTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTg      <  1:1907621/71‑1 (MQ=255)
                                                              cAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCg                            >  1:2598438/1‑47 (MQ=255)
                                                                 gTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATa  >  1:2461451/1‑70 (MQ=255)
                                                                   aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGa                         >  1:5233561/1‑45 (MQ=255)
                                                                   aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTg      >  1:331218/1‑64 (MQ=255)
                                                                   aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATa  >  1:5537256/1‑68 (MQ=255)
                                                                   |                                                                   
TGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATA  >  Exported/834376‑834510

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: