Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,069 A→G 10.9% Y176C (TAC→TGC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,0690AG10.9% 164.8 / 9.4 64Y176C (TAC→TGC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (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 = 1.00e+00

GACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTT  >  Exported/832003‑832137
                                                                  |                                                                    
gACGCGGTGGATTTGTTCGTTCATCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTacac                                                                   >  1:651955/1‑70 (MQ=255)
           tttGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                      <  1:3144920/56‑1 (MQ=255)
           tttGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                      <  1:741651/56‑1 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                     >  1:2212690/1‑51 (MQ=255)
                 cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                     >  1:5235618/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)
                     cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                      >  1:4580473/1‑46 (MQ=255)
                     cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                             >  1:349094/1‑71 (MQ=255)
                        ccTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                             <  1:2145428/68‑1 (MQ=255)
                         cTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCg                                           <  1:4029823/69‑1 (MQ=255)
                            gCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                      >  1:1536165/1‑39 (MQ=255)
                                 aggTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGt                                 <  1:215467/71‑1 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTAc                                                      >  1:352943/1‑46 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                             >  1:4393584/1‑71 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                             >  1:4941474/1‑71 (MQ=255)
                                     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                             >  1:493275/1‑71 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                            >  1:4040059/1‑71 (MQ=255)
                                      aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                            >  1:2628453/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:2022767/67‑1 (MQ=255)
                                                            aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc          <  1:3704031/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:5180480/38‑1 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                    <  1:3938579/38‑1 (MQ=255)
                                                               tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                    <  1:3931721/38‑1 (MQ=255)
                                                                ttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt   >  1:3071241/1‑70 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg            >  1:3673520/1‑60 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt  <  1:2123143/70‑1 (MQ=255)
                                                                 tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt  <  1:3578637/70‑1 (MQ=255)
                                                                  |                                                                    
GACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTT  >  Exported/832003‑832137

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: