Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,232 A→G 19.7% Q230Q (CAA→CAG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,2320AG19.7% 138.8 / 29.0 71Q230Q (CAA→CAG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (6/8);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCC  >  Exported/832163‑832298
                                                                     |                                                                  
ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                                                   >  1:5052580/1‑71 (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)
                     cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                              <  1:1722190/71‑1 (MQ=255)
                       tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                             <  1:3382378/70‑1 (MQ=255)
                       tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt                                             <  1:3375919/70‑1 (MQ=255)
                       tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                         >  1:5256285/1‑58 (MQ=255)
                       tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa                                                         >  1:2239762/1‑58 (MQ=255)
                        aaGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                            <  1:1006427/70‑1 (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:2842309/42‑1 (MQ=255)
                                          cTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                      <  1:3609614/42‑1 (MQ=255)
                                            gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                        >  1:3097322/1‑70 (MQ=255)
                                            gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                        >  1:180735/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:945171/71‑1 (MQ=255)
                                             cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa                      <  1:4997379/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:5526511/1‑70 (MQ=255)
                                                        ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTa            >  1:3282601/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)
                                                                 cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc     <  1:4710229/68‑1 (MQ=255)
                                                                 cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc     <  1:3833639/68‑1 (MQ=255)
                                                                 cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc     <  1:1900274/68‑1 (MQ=255)
                                                                 cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGcc  <  1:5186408/71‑1 (MQ=255)
                                                                 cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGcc  <  1:3036599/71‑1 (MQ=255)
                                                                     |                                                                  
CTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCC  >  Exported/832163‑832298

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: