Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,070 C→G 10.9% Y176* (TAC→TAG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,0700CG10.9% 171.8 / 10.6 64Y176* (TAC→TAG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (3/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTT  >  Exported/832003‑832140
                                                                   |                                                                      
gACGCGGTGGATTTGTTCGTTCATCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTacac                                                                      >  1:651955/1‑70 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                        >  1:5235618/1‑51 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                        >  1:2212690/1‑51 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCaaaaa                                                             >  1:2992601/1‑62 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAAc                                                            >  1:2268023/1‑63 (MQ=255)
                   ttcttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCaaa                                                               <  1:340700/58‑1 (MQ=255)
                     cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                                >  1:349094/1‑71 (MQ=255)
                        ccTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                                <  1:2145428/68‑1 (MQ=255)
                         cTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCg                                              <  1:4029823/69‑1 (MQ=255)
                                 aggTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGt                                    <  1:215467/71‑1 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTAc                                                         >  1:352943/1‑46 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                                >  1:493275/1‑71 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                                >  1:4393584/1‑71 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                                >  1:4941474/1‑71 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                               >  1:2628453/1‑71 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                               >  1:4040059/1‑71 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                                <  1:131151/70‑1 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtc                                 <  1:5481344/69‑1 (MQ=255)
                                        cgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGc                             >  1:3485223/1‑71 (MQ=255)
                                         gaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGcc                            <  1:4976961/71‑1 (MQ=255)
                                            cTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAt                         <  1:4964170/71‑1 (MQ=255)
                                              gATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                                >  1:1732775/1‑46 (MQ=255)
                                               attgccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATccc                      <  1:3152349/68‑1 (MQ=255)
                                               aTCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                           <  1:4243162/50‑1 (MQ=255)
                                               aTCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTcc                      >  1:936587/1‑71 (MQ=255)
                                                tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                            <  1:4168738/48‑1 (MQ=255)
                                                tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAt                         >  1:1866958/1‑67 (MQ=255)
                                                tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                     >  1:2305049/1‑71 (MQ=255)
                                                 cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTg                                         >  1:3150246/1‑50 (MQ=255)
                                                  gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc                   >  1:4053174/1‑71 (MQ=255)
                                                  gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc                   >  1:634764/1‑71 (MQ=255)
                                                  gccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTCGCTGGTTTCTCGCCAATTCCGGc                   <  1:2273519/71‑1 (MQ=255)
                                                    cgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                 >  1:123709/1‑71 (MQ=255)
                                                     gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                            <  1:3474836/43‑1 (MQ=255)
                                                      cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGtt                                   >  1:2335341/1‑51 (MQ=255)
                                                      cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg               <  1:5015116/71‑1 (MQ=255)
                                                        tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTc                       <  1:4360545/61‑1 (MQ=255)
                                                        tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc                   <  1:993397/65‑1 (MQ=255)
                                                        tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc                   <  1:2610061/65‑1 (MQ=255)
                                                        tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc             >  1:4559810/1‑71 (MQ=255)
                                                        gcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc                   <  1:183876/64‑1 (MQ=255)
                                                        gcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTc             <  1:5170591/70‑1 (MQ=255)
                                                          tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTcc                      >  1:5432170/1‑60 (MQ=255)
                                                          tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc                   <  1:2670855/63‑1 (MQ=255)
                                                          tAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGtt                                   <  1:2404767/47‑1 (MQ=255)
                                                          tAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTtct                                >  1:5250444/1‑50 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                           >  1:1158981/1‑38 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                           >  1:1888533/1‑38 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                            >  1:1733049/1‑37 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                            >  1:3044162/1‑37 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCg            >  1:517407/1‑69 (MQ=255)
                                                           aaCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGAt          >  1:317158/1‑71 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                           <  1:2798108/37‑1 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                     <  1:3344827/59‑1 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                 >  1:2961017/1‑63 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc             <  1:3704031/67‑1 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc             <  1:2022767/67‑1 (MQ=255)
                                                             cGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtctc                               >  1:2013047/1‑48 (MQ=255)
                                                              gTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGt       >  1:1435766/1‑71 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGc                                        <  1:1025657/37‑1 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                       <  1:3938579/38‑1 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                       <  1:3931721/38‑1 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                       <  1:5180480/38‑1 (MQ=255)
                                                                ttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt      >  1:3071241/1‑70 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg               >  1:3673520/1‑60 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt     <  1:3578637/70‑1 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt     <  1:2123143/70‑1 (MQ=255)
                                                                   cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCtt  <  1:1033628/71‑1 (MQ=255)
                                                                   cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt   <  1:3327177/70‑1 (MQ=255)
                                                                   |                                                                      
GACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTT  >  Exported/832003‑832140

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: