Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,873 G→A 33.5% P197P (CCG→CCA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported835,8730GA33.5% 88.7 / 66.7 81P197P (CCG→CCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base A (15/12);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGG  >  Exported/835803‑835942
                                                                       |                                                                     
acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                                       <  1:1157043/71‑1 (MQ=255)
acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                                       <  1:3126183/71‑1 (MQ=255)
  acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                                     <  1:1122157/71‑1 (MQ=255)
  acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                                     <  1:2239291/71‑1 (MQ=255)
  acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                                     <  1:3395560/71‑1 (MQ=255)
  acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                                     <  1:2812730/71‑1 (MQ=255)
   ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                                     <  1:5009832/70‑1 (MQ=255)
   ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                                                    >  1:1982968/1‑71 (MQ=255)
      tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGt                                                                 >  1:4422061/1‑71 (MQ=255)
      tCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggc                                                                 >  1:1641382/1‑70 (MQ=255)
       cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggc                                                                 <  1:4020327/70‑2 (MQ=255)
       cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg                                                                  <  1:860058/69‑1 (MQ=255)
        atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg                                                                  <  1:73443/68‑1 (MQ=255)
        atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgc                                                               >  1:3207135/1‑71 (MQ=255)
         tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                                              <  1:5125726/71‑1 (MQ=255)
         tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgct                                                              <  1:5492776/71‑2 (MQ=255)
         tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCACAGCTGCGTGGCGACCTGCCCAAGCGgcgc                                                               >  1:3424952/1‑70 (MQ=255)
           atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                                            <  1:5454100/71‑1 (MQ=255)
           atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                                            <  1:4593942/71‑1 (MQ=255)
           atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgcta                                                             <  1:1889141/70‑3 (MQ=255)
                cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc                                                        >  1:2959457/3‑70 (MQ=255)
                cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc                                                        >  1:965612/3‑70 (MQ=255)
                  tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAg                                                     >  1:4295982/1‑71 (MQ=255)
                   gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                                                    >  1:1914269/1‑69 (MQ=255)
                   gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                                                    >  1:5474769/1‑69 (MQ=255)
                    ccGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                                   >  1:4131440/1‑71 (MQ=255)
                    ccGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGCGACCTGCCCGAGCGGTGCGATTTACAAGCGt                                                  >  1:4347390/1‑71 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                                      >  1:4792923/1‑66 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                                      >  1:3909760/1‑66 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                                      >  1:2057523/1‑66 (MQ=255)
                        gcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTga                                                <  1:1424948/70‑1 (MQ=255)
                        gcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaa                                               <  1:4008683/71‑1 (MQ=255)
                          cTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgct                                                              <  1:1982515/54‑2 (MQ=255)
                           tgtgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaag                                              <  1:5114094/69‑1 (MQ=255)
                            gtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                                               >  1:3075944/1‑67 (MQ=255)
                             tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                                        <  1:2743499/57‑1 (MQ=255)
                             tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:5674081/69‑1 (MQ=255)
                             tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCg                                                                   >  1:2189811/1‑46 (MQ=12)
                             tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                                            <  1:5420271/69‑1 (MQ=255)
                             tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                                            <  1:1690358/69‑1 (MQ=255)
                              gCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                                              >  1:2820836/1‑50 (MQ=255)
                              gCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCAGCGGTGCg                                                              >  1:1934623/1‑49 (MQ=255)
                              gCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                                                     >  1:4783294/1‑58 (MQ=255)
                              gCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                                                     >  1:298998/1‑58 (MQ=255)
                               cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                                   <  1:4920128/60‑1 (MQ=255)
                               cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                                 >  1:1109417/1‑62 (MQ=255)
                                gAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:5406577/1‑68 (MQ=255)
                                 aaCACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                                       >  1:2365420/1‑70 (MQ=255)
                                 aaCACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAagaagacg                                        <  1:809713/69‑3 (MQ=255)
                                     cTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:4544790/61‑1 (MQ=255)
                                     cTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg                                  <  1:3433784/71‑1 (MQ=255)
                                          tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt                              >  1:854048/1‑70 (MQ=255)
                                          tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc                                   >  1:1163641/1‑65 (MQ=255)
                                          tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              >  1:5067144/1‑70 (MQ=255)
                                            aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                              >  1:3137770/1‑52 (MQ=255)
                                            aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGGGAAGAAGATGGCATCGTCCTGAt                           <  1:4859008/71‑1 (MQ=255)
                                               cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                         <  1:4446315/71‑1 (MQ=255)
                                                ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:998695/51‑1 (MQ=255)
                                                 cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                           <  1:3814880/51‑1 (MQ=255)
                                                 cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                       <  1:2858466/71‑1 (MQ=255)
                                                  ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                           <  1:2387615/50‑1 (MQ=255)
                                                  ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                      <  1:4683953/71‑1 (MQ=255)
                                                    cATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                           <  1:4593058/48‑1 (MQ=255)
                                                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:2313650/45‑1 (MQ=255)
                                                      tgcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                  >  1:2592635/4‑71 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:1287185/44‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:1390108/44‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:5629327/44‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:543383/44‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGGGCGATTTACAAGCGTgaagaa                                            <  1:1740031/44‑1 (MQ=255)
                                                       gcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                  >  1:4385442/3‑70 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 >  1:2016248/1‑70 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 >  1:534635/1‑70 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 >  1:1775170/1‑70 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 >  1:4396349/1‑70 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 >  1:456002/1‑70 (MQ=255)
                                                            gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa              >  1:5052181/1‑69 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                               >  1:2232410/1‑47 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg       <  1:4317549/71‑1 (MQ=255)
                                                                   gCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                             >  1:3817871/1‑47 (MQ=255)
                                                                    cccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc     <  1:512664/70‑1 (MQ=255)
                                                                    cccAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc     <  1:2476477/70‑1 (MQ=255)
                                                                      cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg  >  1:5396429/1‑71 (MQ=255)
                                                                       |                                                                     
ACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGG  >  Exported/835803‑835942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: