Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,737 G→T 10.2% P1065P (CCG→CCT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,7370GT10.2% 244.1 / 9.5 78P1065P (CCG→CCTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (27/43);  new base T (6/2);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 6.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAG  >  Exported/834667‑834798
                                                                      |                                                             
gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                               <  1:4491433/71‑1 (MQ=255)
   agcagcTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATc                                                            <  1:2868745/71‑1 (MQ=255)
   agcagcTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATc                                                            <  1:136506/71‑1 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                              >  1:3327716/1‑69 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                              >  1:3815079/1‑69 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                              >  1:1724886/1‑69 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGAt                                                             >  1:154109/1‑70 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGAt                                                             >  1:2261944/1‑70 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                            <  1:4205993/71‑1 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                            <  1:247980/71‑1 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                            <  1:5529276/71‑1 (MQ=255)
   aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                            >  1:2032962/1‑71 (MQ=255)
    gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg                                                           <  1:4420269/71‑1 (MQ=255)
     cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg                                                           <  1:837594/70‑1 (MQ=255)
     cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                          <  1:4825885/71‑1 (MQ=255)
     cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                          <  1:4325118/71‑1 (MQ=255)
       aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca                                                        >  1:4505128/1‑71 (MQ=255)
          gTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAccc                                                     >  1:4966644/1‑71 (MQ=255)
          ctATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACAccc                                                     >  1:3472373/2‑71 (MQ=255)
          ctATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACAcc                                                      >  1:3260613/2‑70 (MQ=255)
           tATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                              <  1:2438894/61‑1 (MQ=255)
             tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTAt                                                             >  1:1388478/1‑60 (MQ=255)
              cAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGtag                                                 >  1:557796/1‑69 (MQ=255)
               aGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGtagc                                                >  1:2117312/1‑68 (MQ=255)
               aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCg                                                <  1:334152/71‑1 (MQ=255)
               aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCg                                                <  1:577294/71‑1 (MQ=255)
               aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCAGTTCg                                                <  1:2738612/71‑1 (MQ=255)
                ggATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGac                                                         <  1:12983/61‑1 (MQ=255)
                 gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt                                              <  1:3629170/71‑1 (MQ=255)
                   tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                            <  1:1917789/71‑1 (MQ=255)
                   tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                            <  1:3104137/71‑1 (MQ=255)
                     aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGaa                                           <  1:3637375/70‑1 (MQ=255)
                       cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaaga                                        <  1:120565/71‑1 (MQ=255)
                       cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaaga                                        <  1:3081619/71‑1 (MQ=255)
                       cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaag                                         >  1:3904650/1‑70 (MQ=255)
                          tGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                            <  1:2354310/48‑1 (MQ=255)
                            gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAcc                                                      >  1:4474786/1‑52 (MQ=255)
                            gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                    <  1:4011875/70‑1 (MQ=255)
                            gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                    <  1:929108/70‑1 (MQ=255)
                             aTGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt                                  <  1:2580975/71‑1 (MQ=255)
                             aTGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt                                  <  1:4214503/71‑1 (MQ=255)
                                cGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:2438277/65‑1 (MQ=255)
                                cGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:3909307/65‑1 (MQ=255)
                                 gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca                                                        >  1:1074618/1‑45 (MQ=255)
                                 gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGc                              >  1:5599783/1‑71 (MQ=255)
                                 gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGc                              <  1:4730299/71‑1 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                               <  1:2063339/37‑1 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                               <  1:415935/37‑1 (MQ=255)
                                    cttttGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGta                                                  >  1:868645/2‑47 (MQ=255)
                                     tttCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGa                          <  1:2058021/71‑1 (MQ=255)
                                      ttCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg                                                           >  1:4734083/1‑37 (MQ=255)
                                          gTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATcc                     <  1:3319256/71‑1 (MQ=255)
                                           tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   <  1:4452795/56‑1 (MQ=255)
                                           tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   <  1:4493381/56‑1 (MQ=255)
                                             aaaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaagaa                                       >  1:3234401/1‑50 (MQ=255)
                                              aaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:4145517/51‑1 (MQ=255)
                                                gCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGc               <  1:4618459/71‑1 (MQ=255)
                                                  ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa                   >  1:2953691/1‑65 (MQ=255)
                                                  ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa                   >  1:5376213/1‑65 (MQ=255)
                                                  ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa                   >  1:2972975/1‑65 (MQ=255)
                                                   tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:4062973/46‑1 (MQ=255)
                                                   tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:4365230/46‑1 (MQ=255)
                                                   tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:39515/46‑1 (MQ=255)
                                                    gctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc            >  1:5526011/1‑70 (MQ=255)
                                                    gctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc            >  1:5248953/1‑70 (MQ=255)
                                                     ctgGTTTATCTTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   >  1:3998311/1‑46 (MQ=255)
                                                     ctgGTTTATCTTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   >  1:2626553/1‑46 (MQ=255)
                                                     ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa                        <  1:5093938/57‑1 (MQ=255)
                                                     ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc            <  1:1172881/69‑1 (MQ=255)
                                                      tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                     <  1:5152740/43‑1 (MQ=255)
                                                      tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   <  1:3383669/45‑1 (MQ=255)
                                                      tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   <  1:2605016/45‑1 (MQ=255)
                                                         tttATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                    <  1:1546638/41‑1 (MQ=255)
                                                          ttATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa      >  1:3051587/1‑70 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   >  1:4965810/1‑40 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa                                   >  1:4625639/1‑40 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa     >  1:2027858/1‑70 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa     <  1:5466708/70‑1 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa      >  1:1439398/1‑69 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa      >  1:1148456/1‑69 (MQ=255)
                                                             tCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAg  <  1:668103/71‑1 (MQ=255)
                                                                 ccgccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaaaa   >  1:3365911/1‑66 (MQ=255)
                                                                     cgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa                        >  1:653541/1‑41 (MQ=255)
                                                                     cgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa                        >  1:2422102/1‑41 (MQ=255)
                                                                      |                                                             
GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAG  >  Exported/834667‑834798

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: