Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 833,133 G→C 13.6% D531H (GAT→CAT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported833,1330GC13.6% 188.6 / 17.4 73D531H (GAT→CAT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (8/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGA  >  Exported/833067‑833197
                                                                  |                                                                
gAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaa                                                              <  1:998684/71‑1 (MQ=255)
gAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaa                                                              <  1:812631/71‑1 (MQ=255)
 aGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATaaa                                                             <  1:392896/71‑1 (MQ=255)
   cAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAAc                                                           >  1:4503973/1‑71 (MQ=255)
          ccGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACgg                                                     >  1:2061527/1‑70 (MQ=255)
          ccGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGt                                                    >  1:989685/1‑71 (MQ=255)
            ggCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  >  1:1446035/1‑71 (MQ=255)
                tttCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5458287/71‑1 (MQ=255)
                    tcGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatga                                          >  1:2032117/1‑71 (MQ=255)
                    tcGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatga                                          >  1:4160535/1‑71 (MQ=255)
                     cGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatgat                                         >  1:4495796/1‑71 (MQ=255)
                      gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCgatgat                                         >  1:3294832/1‑70 (MQ=255)
                      gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGAtt                                        <  1:1660017/71‑1 (MQ=255)
                      gCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGAtt                                        <  1:4631708/71‑1 (MQ=255)
                            aGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:4424897/1‑71 (MQ=255)
                              attattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtcg                                 >  1:3876212/1‑70 (MQ=255)
                              attattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCgg                                >  1:5486270/1‑71 (MQ=255)
                                tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:129051/1‑36 (MQ=255)
                                tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:3794782/1‑36 (MQ=255)
                                tattCGCATCGCCCGTGAATTTGCCGATAACGCTGa                                                                 >  1:2120276/1‑36 (MQ=255)
                                tattCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGt                                                    >  1:3776606/1‑49 (MQ=255)
                                 attCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGc                             >  1:2310739/1‑71 (MQ=255)
                                  ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGc                             <  1:4779345/70‑1 (MQ=255)
                                  ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCg                            >  1:448964/1‑71 (MQ=255)
                                  ttCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCg                            >  1:4875783/1‑71 (MQ=255)
                                  ttCGCATCGCACGTGAATTTGCCGATAACGCTGATaa                                                              <  1:388856/37‑1 (MQ=255)
                                     ccATCGCCCGTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCagg                          >  1:3470908/2‑70 (MQ=255)
                                        tCGCCCGTGAATTTGCCGATACCGCCCATAAAACgt                                                         <  1:5321883/36‑2 (MQ=21)
                                        tCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCt                       >  1:5553360/1‑70 (MQ=255)
                                          gCCCGTGAATTTGCCGATAGCGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:4249936/1‑57 (MQ=255)
                                          gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:1539867/1‑57 (MQ=255)
                                          gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGa                     <  1:4305386/70‑1 (MQ=255)
                                           cccGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAAc                   <  1:1926014/71‑1 (MQ=255)
                                              gTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCgatg                                           >  1:4957753/1‑44 (MQ=16)
                                              gTGAATTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCAc                <  1:3346124/71‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:2804959/36‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:776077/36‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCg                                                  <  1:3733505/36‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5602389/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4947925/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4379060/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:5063282/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4147449/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:4130145/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:1478583/40‑1 (MQ=255)
                                               tGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCg                                              <  1:2529688/40‑1 (MQ=255)
                                                gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:725909/63‑1 (MQ=255)
                                                gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:2100598/63‑1 (MQ=255)
                                                gAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTg              >  1:3616516/1‑71 (MQ=255)
                                                 aaTTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtcctc                                  >  1:1112654/1‑50 (MQ=255)
                                                 aaTTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtcctc                                  >  1:2958871/1‑50 (MQ=255)
                                                 aaTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTg                      <  1:3599317/62‑1 (MQ=255)
                                                 aaTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTgg             >  1:3454848/1‑71 (MQ=255)
                                                  aTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTAtcgtc                                  >  1:963441/1‑49 (MQ=255)
                                                  aTTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGt            <  1:5615322/71‑1 (MQ=255)
                                                    ttGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCgg                                >  1:2504301/1‑49 (MQ=255)
                                                     tggcGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTAt          >  1:5308718/4‑70 (MQ=255)
                                                     tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTAtc                                     >  1:1594714/1‑43 (MQ=16)
                                                     tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCaggtg                        >  1:964578/1‑54 (MQ=255)
                                                     tGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCa                 >  1:3023558/1‑63 (MQ=255)
                                                     tGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTAt          >  1:2204781/1‑70 (MQ=255)
                                                          aTAACGCTGATAAAACGCACGGTCGTTCGATGATTAt                                      >  1:5328840/1‑37 (MQ=255)
                                                          aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:511904/1‑71 (MQ=255)
                                                          aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:5594759/1‑71 (MQ=255)
                                                          aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:3836548/1‑71 (MQ=255)
                                                          aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    >  1:2843909/1‑71 (MQ=255)
                                                            aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:4694658/69‑1 (MQ=255)
                                                            aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:3817943/69‑1 (MQ=255)
                                                            aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc    <  1:3811269/69‑1 (MQ=255)
                                                            aaCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGa  >  1:3457660/1‑71 (MQ=255)
                                                             aCGCTGATAAAACGCACGGTCGTTCGATGATTATCGTAGGTGGGGGGCTGAACCACTGGTATCACCTc    <  1:5588011/68‑1 (MQ=255)
                                                              cGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCg   >  1:3274504/1‑68 (MQ=255)
                                                               gCTGATAAAACGCACGGTCGTTCGATGATTAtcgt                                   <  1:2252131/35‑1 (MQ=255)
                                                                  |                                                                
GAGCAGATTACCGGCGTTTCTCGCAGCCAGATTATTCGCATCGCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGA  >  Exported/833067‑833197

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: