Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,226 T→G 21.9% S228S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,2260TG21.9% 107.0 / 43.0 73S228S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base G (6/10);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCG  >  Exported/832157‑832296
                                                                     |                                                                      
cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                                                       <  1:1356289/71‑1 (MQ=255)
   tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                                                                     <  1:3323946/70‑1 (MQ=255)
   tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:5193665/1‑69 (MQ=255)
    ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                    >  1:5165428/1‑70 (MQ=255)
      ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:1235824/1‑66 (MQ=255)
      ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:1526610/1‑66 (MQ=255)
      ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                                                 >  1:5052580/1‑71 (MQ=255)
       tcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgatg                                                                  <  1:761100/69‑4 (MQ=255)
       tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                                >  1:3511473/1‑71 (MQ=255)
         gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                               <  1:2597733/70‑1 (MQ=255)
         gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                               <  1:569964/70‑1 (MQ=255)
            gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                                                               <  1:3422109/67‑1 (MQ=255)
              ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                           >  1:5620596/1‑69 (MQ=255)
              ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                         <  1:1200017/71‑1 (MQ=255)
                ggCGGTACATGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                        <  1:1009608/70‑1 (MQ=255)
                    gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCaa                                                   <  1:3585810/71‑1 (MQ=255)
                      aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                    <  1:663507/68‑1 (MQ=255)
                       cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAAc                                                 <  1:4644879/70‑1 (MQ=255)
                       cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                                <  1:163062/71‑1 (MQ=255)
                        ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                                  <  1:5615593/52‑1 (MQ=255)
                         tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:1723349/1‑47 (MQ=255)
                         tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:4719233/1‑47 (MQ=255)
                         tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                                                      >  1:1707352/1‑47 (MQ=255)
                           cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                            <  1:1722190/71‑1 (MQ=255)
                             tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                           <  1:3382378/70‑1 (MQ=255)
                             tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                           <  1:3375919/70‑1 (MQ=255)
                             tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                       >  1:2239762/1‑58 (MQ=255)
                             tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                       >  1:5256285/1‑58 (MQ=255)
                              aaGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                          <  1:1006427/70‑1 (MQ=255)
                              aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                    >  1:3515463/1‑44 (MQ=255)
                              aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGt                                           <  1:4166571/69‑1 (MQ=255)
                                 tttCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGc                                                         <  1:3806603/51‑1 (MQ=255)
                                   tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa                                    >  1:1655704/1‑71 (MQ=255)
                                   tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                     <  1:3303868/70‑1 (MQ=255)
                                    cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                     <  1:999977/69‑1 (MQ=255)
                                    cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                     <  1:2123804/69‑1 (MQ=255)
                                    cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                   >  1:3042141/1‑71 (MQ=255)
                                     tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                   >  1:1807426/1‑70 (MQ=255)
                                     tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc                                  >  1:1283628/1‑71 (MQ=255)
                                     tACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATc                                  >  1:1396551/1‑71 (MQ=255)
                                      aCGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                   >  1:4469270/1‑69 (MQ=255)
                                       cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATc                                  <  1:5684819/69‑1 (MQ=255)
                                       cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                   <  1:90272/68‑1 (MQ=255)
                                       cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATct                                 >  1:4196601/1‑70 (MQ=255)
                                          cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg                             <  1:1065204/71‑1 (MQ=255)
                                           tGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCgg                                      <  1:1009180/61‑1 (MQ=255)
                                            ggTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                           >  1:5140846/1‑71 (MQ=255)
                                                ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGaa                                    <  1:2171945/57‑1 (MQ=255)
                                                ttGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                           >  1:1859496/2‑67 (MQ=255)
                                                cTGCGTCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACt                          <  1:1194552/68‑1 (MQ=255)
                                                cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                       >  1:3899584/1‑39 (MQ=255)
                                                cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                    <  1:3609614/42‑1 (MQ=255)
                                                cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                    <  1:2842309/42‑1 (MQ=255)
                                                  gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                      >  1:180735/1‑70 (MQ=255)
                                                  gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                      >  1:3097322/1‑70 (MQ=255)
                                                  gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGa                                              >  1:3131241/1‑46 (MQ=255)
                                                  gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                                   >  1:5092699/1‑57 (MQ=255)
                                                  gCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAAt                                   >  1:3352922/1‑57 (MQ=255)
                                                  gCGACATGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTa                      >  1:2420630/1‑70 (MQ=255)
                                                   cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTg                         >  1:2286052/1‑66 (MQ=255)
                                                   cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                     <  1:3566441/70‑1 (MQ=255)
                                                   cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                     <  1:1238824/70‑1 (MQ=255)
                                                   cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa                    <  1:4997379/71‑1 (MQ=255)
                                                   cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa                    <  1:945171/71‑1 (MQ=255)
                                                    gACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                   <  1:3450789/71‑1 (MQ=255)
                                                     aCTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAAc                  >  1:3771725/1‑71 (MQ=255)
                                                       ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                >  1:418955/1‑71 (MQ=255)
                                                          cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAg             >  1:2857378/1‑71 (MQ=255)
                                                           ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                          >  1:4872348/1‑41 (MQ=255)
                                                           ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTa                                          >  1:4903174/1‑41 (MQ=255)
                                                              ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa          >  1:3282601/1‑70 (MQ=255)
                                                              ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa          >  1:5526511/1‑70 (MQ=255)
                                                              ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGATa          >  1:2386521/1‑70 (MQ=255)
                                                                 cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat       >  1:5357351/1‑70 (MQ=255)
                                                                     tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGaa                                    >  1:5625362/1‑37 (MQ=255)
                                                                     tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGaa                                    >  1:267698/1‑37 (MQ=255)
                                                                     tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                     <  1:5655540/36‑1 (MQ=255)
                                                                     tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                     <  1:3440159/36‑1 (MQ=255)
                                                                     tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg  >  1:1374314/1‑71 (MQ=255)
                                                                     |                                                                      
CTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCG  >  Exported/832157‑832296

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: