Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 833,912 G→C 28.2% A790A (GCG→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported833,9120GC28.2% 129.3 / 59.6 85A790A (GCG→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (36/25);  new base C (16/8);  total (52/33)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

GTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGC  >  Exported/833847‑833980
                                                                 |                                                                    
gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGActt                                                                  <  1:3872112/70‑3 (MQ=255)
      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg                                                            <  1:2530617/70‑1 (MQ=255)
      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg                                                            >  1:2414707/1‑70 (MQ=255)
      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg                                                            <  1:4562829/70‑1 (MQ=255)
      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACg                                                            <  1:1473821/70‑1 (MQ=255)
       cGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa                                                          <  1:774887/71‑1 (MQ=255)
        gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa                                                         <  1:5376024/71‑1 (MQ=255)
        gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                                                         <  1:3798898/71‑1 (MQ=255)
          tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaa                                                        >  1:1865983/3‑70 (MQ=255)
           cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa                                                          >  1:2054885/2‑67 (MQ=255)
            gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa                                                       <  1:2869982/69‑1 (MQ=255)
            gACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa                                                       <  1:227333/69‑1 (MQ=255)
            gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaaaa                                                       <  1:2516655/69‑1 (MQ=255)
            gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga                                                     >  1:3490634/1‑71 (MQ=255)
            gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga                                                     >  1:3497426/1‑71 (MQ=255)
               ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAgacga                                                  >  1:2829301/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                                                                    >  1:4541119/1‑51 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:5451873/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:3799207/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:3333857/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:3290762/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:252478/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:1729023/1‑71 (MQ=255)
                 ccGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGCCGACCTGGTACGAAAAAGACGACa                                                >  1:3451766/1‑71 (MQ=255)
                  cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa                                                         >  1:4675448/1‑61 (MQ=255)
                  cGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:590478/1‑70 (MQ=255)
                   gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaag                                                         >  1:5404182/1‑59 (MQ=255)
                   gTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaa                                                        >  1:1978027/1‑61 (MQ=255)
                        tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                                                                    >  1:2661012/1‑44 (MQ=255)
                          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:460371/1‑62 (MQ=255)
                          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:3676956/1‑62 (MQ=255)
                          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:4517095/1‑62 (MQ=255)
                          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACa                                                >  1:2271004/1‑62 (MQ=255)
                           aGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTc                                      >  1:5319269/1‑71 (MQ=255)
                            gCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCg                                     >  1:2561185/1‑71 (MQ=255)
                                 tGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa                                                       <  1:5502763/48‑1 (MQ=255)
                                  gTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATg                               >  1:5295264/1‑71 (MQ=255)
                                   tctcTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa                                                         <  1:4890858/44‑1 (MQ=255)
                                   tctcTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATAc                                         >  1:4163088/1‑60 (MQ=255)
                                   tctcTATTCCGACATCATTTTGCCGACCGCGACCTg                                                                 >  1:2024091/1‑36 (MQ=255)
                                       tATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCa                             <  1:2424132/68‑1 (MQ=255)
                                       tATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATcc                          >  1:2999031/1‑71 (MQ=255)
                                         ttCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCg                         <  1:3539171/70‑1 (MQ=255)
                                         ctCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAg                                                      <  1:541957/40‑1 (MQ=12)
                                         ctCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGaa                                            >  1:1923086/2‑51 (MQ=255)
                                         ctCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGt                        >  1:141150/2‑71 (MQ=255)
                                          tCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGc                              >  1:983598/1‑64 (MQ=255)
                                          tCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGc                              >  1:3065344/1‑64 (MQ=255)
                                           ccGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAgac                                                    >  1:3400195/1‑41 (MQ=16)
                                           ccGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa                                            >  1:2381677/1‑49 (MQ=255)
                                           ccGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa                                            >  1:1487850/1‑49 (MQ=255)
                                           ccGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGtt                       <  1:5072479/70‑1 (MQ=255)
                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAgac                                                    <  1:384137/40‑1 (MQ=16)
                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAgac                                                    <  1:416510/40‑1 (MQ=16)
                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATg                                              <  1:4124788/46‑1 (MQ=255)
                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTa                     >  1:4247396/1‑71 (MQ=255)
                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTa                     >  1:2054840/1‑71 (MQ=255)
                                            cGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTa                     >  1:2743258/1‑71 (MQ=255)
                                             gatatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTa                     <  1:1043388/67‑1 (MQ=255)
                                             gatatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTAt                    >  1:892632/4‑71 (MQ=255)
                                             gatatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTAt                    >  1:3897985/4‑71 (MQ=255)
                                             gatatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTAt                    >  1:892630/4‑71 (MQ=255)
                                             gatatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTAt                    >  1:3568848/4‑71 (MQ=255)
                                              aCATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATa                                          >  1:1793678/1‑48 (MQ=255)
                                              aCATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTa                     >  1:426848/1‑69 (MQ=255)
                                               catcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTc                  <  1:5664120/71‑1 (MQ=255)
                                                atcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGata                                 >  1:4274569/1‑55 (MQ=255)
                                                atcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGata                                 >  1:3466423/1‑55 (MQ=255)
                                                atcgtTCTGCCCACCGCCACCTGCTACGAAAAAGACGATATGaa                                            >  1:4673438/1‑44 (MQ=12)
                                                 tcatTTTGCCGACGGCGACCTGGTACGAAAAAgac                                                    <  1:1829124/35‑1 (MQ=255)
                                                  catTTTGCCGACGGCGACCTGGTACGAAAAAGAAGACATGAATACTTCGGATATGCATCCGTTTATTCAcc               >  1:2473935/1‑71 (MQ=255)
                                                  caatttGCCGACGGCGACCTGGTACGAAAAAGACTACATGAATACTTCGGATATGCATCCGTTTATTCAcc               >  1:1915154/4‑71 (MQ=255)
                                                   atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCg                                     <  1:1370847/48‑1 (MQ=255)
                                                   atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGa                                   <  1:4351762/50‑1 (MQ=255)
                                                   atTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAccc              >  1:1279436/1‑71 (MQ=255)
                                                     tttGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa                                            <  1:4593526/39‑1 (MQ=255)
                                                        gCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg          <  1:3934286/70‑1 (MQ=255)
                                                         ccGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGata                                 <  1:269858/46‑1 (MQ=255)
                                                            aCGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGata                                 <  1:2219525/43‑1 (MQ=255)
                                                             cGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcg    >  1:182574/1‑71 (MQ=255)
                                                              tgCGACGTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGtt                       <  1:5074128/50‑1 (MQ=255)
                                                              ggCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGtt                       <  1:3505476/51‑1 (MQ=255)
                                                               gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc  <  1:4778080/71‑1 (MQ=255)
                                                               gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc  <  1:2071345/71‑1 (MQ=255)
                                                                cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg          <  1:1610348/62‑1 (MQ=255)
                                                                 gACCTGGTACGAAAAAGACGACATGAATACTTCGGata                                 <  1:3961289/38‑1 (MQ=255)
                                                                 gACCTGGTACGAAAAAGACGACATGAATACTTCGGata                                 <  1:1009580/38‑1 (MQ=255)
                                                                 |                                                                    
GTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGC  >  Exported/833847‑833980

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: