Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 833,252 G→A 29.6% Q570Q (CAG→CAA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported833,2520GA29.6% 95.2 / 49.0 71Q570Q (CAG→CAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (31/19);  new base A (15/6);  total (46/25)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

ACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  Exported/833182‑833319
                                                                      |                                                                   
aCTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:3824254/71‑1 (MQ=255)
    gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                                                                   <  1:763299/69‑1 (MQ=255)
       tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                                  <  1:3514023/67‑1 (MQ=255)
       tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCCGCTGTGTCGGTCAGAGc                                                                  <  1:577276/67‑1 (MQ=255)
       tCACCTCGATATGAACTATCGTGATCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                                  <  1:3016432/67‑1 (MQ=255)
         aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                                                                >  1:1604126/1‑67 (MQ=255)
             cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        >  1:832257/1‑71 (MQ=255)
              gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTg                                                       >  1:1200263/1‑71 (MQ=255)
               atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                                                                >  1:3879120/1‑61 (MQ=255)
               atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGAgg                                                                >  1:3948894/1‑61 (MQ=255)
                tatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTggg                                                     <  1:4653097/71‑1 (MQ=255)
                  tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:2083719/53‑1 (MQ=255)
                   gAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                                                               >  1:47930/1‑58 (MQ=255)
                     aCTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCAc                                                >  1:3275974/1‑71 (MQ=255)
                      cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                                  >  1:3281777/1‑52 (MQ=255)
                       tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:2191825/48‑1 (MQ=255)
                       tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:2566668/48‑1 (MQ=255)
                       tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:2477942/48‑1 (MQ=255)
                        aTCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTAt                                             >  1:1231917/1‑71 (MQ=255)
                         tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTa                                              <  1:2326530/69‑1 (MQ=255)
                         tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                            >  1:2480508/1‑71 (MQ=255)
                         tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                            >  1:407809/1‑71 (MQ=255)
                         tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGAGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                            <  1:472933/71‑1 (MQ=255)
                              gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                       <  1:1916905/71‑1 (MQ=255)
                              gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                       <  1:4566378/71‑1 (MQ=255)
                               tCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                       >  1:806225/1‑70 (MQ=255)
                                 tGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:1527186/38‑1 (MQ=255)
                                   aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAg                                    >  1:2821760/1‑69 (MQ=255)
                                   aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGa                                  <  1:2944028/71‑1 (MQ=255)
                                   aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                                                                   <  1:2509602/38‑1 (MQ=255)
                                    tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:3857956/1‑58 (MQ=255)
                                    tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:2691954/1‑58 (MQ=255)
                                    tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:3851377/1‑58 (MQ=255)
                                    tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:49515/1‑58 (MQ=255)
                                    tCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTa                                              >  1:5152176/1‑58 (MQ=255)
                                     cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgc                                                    >  1:5581015/1‑51 (MQ=255)
                                      aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa                               <  1:3016967/68‑1 (MQ=255)
                                      aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaa                               <  1:4184778/68‑1 (MQ=255)
                                      aataTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaa                                >  1:2981445/4‑70 (MQ=255)
                                      aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                <  1:3493599/70‑1 (MQ=255)
                                      aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                <  1:2167987/70‑1 (MQ=255)
                                         aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                         >  1:1262426/1‑58 (MQ=255)
                                         aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                         >  1:4690972/1‑58 (MQ=255)
                                         aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAg                                         >  1:188112/1‑58 (MQ=255)
                                         aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:3406557/1‑69 (MQ=255)
                                         aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:4464413/1‑69 (MQ=255)
                                         aTGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:4304371/1‑69 (MQ=255)
                                          tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                                             >  1:1402031/1‑37 (MQ=255)
                                          tGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                                             >  1:1931444/1‑37 (MQ=255)
                                          tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACt                                               <  1:5102118/51‑1 (MQ=255)
                                          tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:4945892/1‑68 (MQ=255)
                                          tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                           >  1:3034003/1‑71 (MQ=255)
                                               atcttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGaaaaa                              >  1:3438767/4‑63 (MQ=255)
                                               aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:3063288/37‑1 (MQ=255)
                                               aTTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:5604210/37‑1 (MQ=255)
                                                ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:2520814/36‑1 (MQ=255)
                                                ttttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:4352605/36‑1 (MQ=255)
                                                 tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTgg                                                      >  1:2816984/1‑37 (MQ=255)
                                                 tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCg                     >  1:1400618/1‑70 (MQ=255)
                                                 cttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGccca                     >  1:949659/2‑69 (MQ=255)
                                                   tCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCa                   <  1:5222416/70‑1 (MQ=255)
                                                     tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCa                                                 >  1:1644431/1‑38 (MQ=255)
                                                     tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:2936667/1‑58 (MQ=255)
                                                     tGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:4963301/1‑58 (MQ=255)
                                                      gCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                >  1:3734738/1‑54 (MQ=255)
                                                        ggCTGTTTCGGTCAGAGCGGGGGCGGCTGGGCGCTCTATGTa                                          >  1:2750740/1‑42 (MQ=38)
                                                        ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                       >  1:3400123/1‑45 (MQ=255)
                                                        ggCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt                                       >  1:3760021/1‑45 (MQ=255)
                                                            gtgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg         >  1:5529711/1‑71 (MQ=255)
                                                             tgtCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg         <  1:1003885/70‑1 (MQ=255)
                                                             tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTggt        <  1:4121250/71‑2 (MQ=255)
                                                             tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg         <  1:1119725/70‑1 (MQ=255)
                                                               tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:1166919/1‑48 (MQ=255)
                                                               tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                  <  1:2999885/59‑1 (MQ=255)
                                                               tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                  <  1:1294849/59‑1 (MQ=255)
                                                                cGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGc     >  1:5440272/1‑71 (MQ=255)
                                                                 ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:3649610/1‑46 (MQ=255)
                                                                 ggTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                             >  1:3089178/1‑46 (MQ=255)
                                                                   tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGcccc  <  1:3377181/71‑3 (MQ=255)
                                                                      |                                                                   
ACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  Exported/833182‑833319

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