Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,385 C→T 7.4% V281V (GTC→GTT) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,3850CT 7.4% 236.7 / 6.4 81V281V (GTC→GTT) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base T (4/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGG  >  Exported/832321‑832455
                                                                |                                                                      
cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                                                                  <  1:568557/71‑1 (MQ=255)
 gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                                                    >  1:4665990/1‑68 (MQ=255)
   gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGAc                                                               <  1:962101/71‑1 (MQ=255)
      cgccGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTa                                                             >  1:1199696/3‑70 (MQ=255)
         cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                         >  1:2756726/1‑71 (MQ=255)
         cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                         >  1:3479817/1‑71 (MQ=255)
         cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActc                                                          >  1:259288/1‑68 (MQ=255)
          ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGc                                                          >  1:1864301/1‑69 (MQ=255)
          ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                                                        <  1:5512562/71‑1 (MQ=255)
          ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                                                        <  1:2691524/71‑1 (MQ=255)
           gATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTg                                                        >  1:689125/1‑70 (MQ=255)
            aTGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                                                      <  1:4310725/71‑1 (MQ=255)
             tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                     <  1:2307398/71‑1 (MQ=255)
             tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                                                      >  1:3174580/1‑70 (MQ=255)
               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                       >  1:5423924/1‑51 (MQ=255)
               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                       >  1:5139528/1‑51 (MQ=255)
               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                                                       >  1:3486655/1‑51 (MQ=255)
                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                                                  >  1:3303168/1‑71 (MQ=255)
                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                                                             <  1:3966909/59‑1 (MQ=255)
                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                                                             <  1:5307474/59‑1 (MQ=255)
                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGa                                                                >  1:5292437/1‑55 (MQ=255)
                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                                 >  1:2624974/1‑70 (MQ=255)
                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                                 >  1:2640918/1‑70 (MQ=255)
                     tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                             >  1:2717180/1‑71 (MQ=255)
                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                             <  1:2724908/70‑1 (MQ=255)
                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                            >  1:3395995/1‑71 (MQ=255)
                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                            >  1:1454583/1‑71 (MQ=255)
                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                            >  1:270543/1‑71 (MQ=255)
                      cttTACTGAAGTGCGTTACAAAGGGACAAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                                                   >  1:3524381/1‑64 (MQ=255)
                         tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtg                                          <  1:2003354/70‑1 (MQ=255)
                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa                                              <  1:4275270/65‑1 (MQ=255)
                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg                                        >  1:2901675/1‑71 (MQ=255)
                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg                                        >  1:5422094/1‑71 (MQ=255)
                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                         >  1:88103/1‑70 (MQ=255)
                           cTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGc                                       <  1:5286533/71‑1 (MQ=255)
                             gAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                                     >  1:5670754/1‑71 (MQ=255)
                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                            >  1:5057104/1‑63 (MQ=255)
                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                                      <  1:1851614/69‑1 (MQ=255)
                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                                      <  1:1822553/69‑1 (MQ=255)
                               aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                                   <  1:1794494/71‑1 (MQ=255)
                               aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                                   >  1:1294859/1‑71 (MQ=255)
                                 tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                                     <  1:517313/67‑1 (MQ=255)
                                 tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                 >  1:2701369/1‑71 (MQ=255)
                                   cGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                     >  1:988673/1‑49 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGtaa                                                     <  1:2958244/47‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                     <  1:95249/47‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                     <  1:2893078/47‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                         <  1:1713914/59‑1 (MQ=255)
                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                         <  1:617241/59‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                >  1:4834499/1‑50 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                >  1:4828489/1‑50 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                                                   >  1:1235976/1‑47 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGcc                                                >  1:2869033/1‑50 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                             >  1:1711066/1‑53 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                               >  1:5662090/1‑51 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                               <  1:5438222/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                               <  1:3993389/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                               <  1:3783935/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                               <  1:638277/51‑1 (MQ=255)
                                       aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggct                            <  1:2914552/70‑1 (MQ=255)
                                         aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                          >  1:2249693/1‑70 (MQ=255)
                                         aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                          >  1:4894099/1‑70 (MQ=255)
                                          aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                          <  1:4086787/69‑1 (MQ=255)
                                           aGGGACCGAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAc                       <  1:572442/71‑1 (MQ=255)
                                             ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                      >  1:4125620/1‑70 (MQ=255)
                                             ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                      >  1:1196310/1‑70 (MQ=255)
                                              gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                    <  1:403498/71‑1 (MQ=255)
                                              gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                    <  1:2448142/71‑1 (MQ=255)
                                               aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa                   <  1:4646452/71‑1 (MQ=255)
                                               aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa                   <  1:2999605/71‑1 (MQ=255)
                                                ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                             <  1:63704/44‑1 (MQ=255)
                                                ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                  <  1:4162029/71‑1 (MQ=255)
                                                  aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                  >  1:1259197/1‑69 (MQ=255)
                                                   aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctggc                         <  1:2482576/61‑1 (MQ=255)
                                                   aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACa                >  1:5118343/1‑70 (MQ=255)
                                                    aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                            >  1:2172988/1‑41 (MQ=255)
                                                        gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc          >  1:4448240/1‑71 (MQ=255)
                                                         ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg              <  1:3690683/66‑1 (MQ=255)
                                                         ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg              <  1:669984/66‑1 (MQ=255)
                                                         ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc         <  1:1910169/71‑1 (MQ=255)
                                                            ccGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt      <  1:3196267/71‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                <  1:3756818/58‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                <  1:894360/58‑1 (MQ=255)
                                                               tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGa       <  1:3866823/67‑1 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                 >  1:2355399/1‑40 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                 >  1:2958692/1‑40 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc            >  1:2585835/1‑61 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg  >  1:4275793/1‑71 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg  >  1:2149322/1‑71 (MQ=255)
                                                                cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg  >  1:2155976/1‑71 (MQ=255)
                                                                |                                                                      
CGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGG  >  Exported/832321‑832455

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: