Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,790 A→G 12.5% S170G (AGC→GGC) ‡ 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,7900AG12.5% 158.4 / 12.5 64S170G (AGC→GGC) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base G (6/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  Exported/835721‑835860
                                                                     |                                                                      
gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg                                                                       >  1:4795866/1‑71 (MQ=255)
gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg                                                                       >  1:795753/1‑71 (MQ=255)
 aCCTGGGCGGTGTGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc                                                                      <  1:1143967/71‑1 (MQ=255)
          gTGAGTTTGACAAACTGGCGAAATACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                             <  1:2148505/71‑1 (MQ=255)
           tGAGTTTGACAAACTGGCGAAAGACAATAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa                                                              <  1:5368627/69‑1 (MQ=255)
            gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc                                                                      <  1:2461712/60‑1 (MQ=255)
            gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           >  1:3329932/1‑71 (MQ=255)
            gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           >  1:1527763/1‑71 (MQ=255)
                 tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCt                                                      >  1:2327926/1‑71 (MQ=255)
                 tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCt                                                      >  1:832129/1‑71 (MQ=255)
                   aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc                                                    <  1:1505982/71‑1 (MQ=255)
                     aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:4764202/71‑1 (MQ=255)
                     aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:4135534/71‑1 (MQ=255)
                          ggCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg                                              >  1:4345582/1‑70 (MQ=255)
                           gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg                                              <  1:194374/69‑1 (MQ=255)
                           gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                            <  1:3107848/71‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           <  1:693211/53‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           <  1:3098997/53‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         <  1:945759/71‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         <  1:3531289/71‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         >  1:44376/1‑71 (MQ=255)
                                aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                             <  1:91621/49‑1 (MQ=255)
                                aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       >  1:4756403/1‑71 (MQ=255)
                                aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       >  1:2123703/1‑71 (MQ=255)
                                      gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                            >  1:2749590/1‑60 (MQ=255)
                                      gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                            >  1:2371877/1‑60 (MQ=255)
                                      caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCgcgc                                  >  1:594367/2‑70 (MQ=255)
                                      caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                 >  1:3923764/2‑71 (MQ=255)
                                        aCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt                                           >  1:125370/1‑59 (MQ=255)
                                        aCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt                                           >  1:1884637/1‑59 (MQ=255)
                                         cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                    <  1:4191217/65‑1 (MQ=255)
                                           tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc                                                    >  1:1805003/1‑47 (MQ=255)
                                           tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt                                           >  1:2531262/1‑56 (MQ=255)
                                           tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgtg                             >  1:773307/1‑70 (MQ=255)
                                           tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                            >  1:2789212/1‑71 (MQ=255)
                                            cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                  <  1:3201256/64‑1 (MQ=255)
                                            cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                            >  1:631081/1‑70 (MQ=255)
                                              acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGag                         <  1:686382/71‑2 (MQ=255)
                                              acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                         <  1:5528516/71‑1 (MQ=255)
                                              acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                         <  1:1707990/71‑1 (MQ=255)
                                               caacaTGCAGAAGGCGATGTATAGCCAGTTCGAAAacac                                                        >  1:2496/1‑39 (MQ=255)
                                               caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAacac                                                        >  1:3657043/1‑39 (MQ=255)
                                                aacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca                       <  1:4207407/71‑1 (MQ=255)
                                                      caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                 >  1:2853531/4‑71 (MQ=255)
                                                      caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                 >  1:432562/4‑71 (MQ=255)
                                                      cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                  <  1:1291312/70‑1 (MQ=255)
                                                      cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 >  1:5462064/1‑71 (MQ=255)
                                                      cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 <  1:5149149/71‑1 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 >  1:4348567/1‑70 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:3186099/71‑1 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:1663632/71‑1 (MQ=255)
                                                       aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:1781324/71‑1 (MQ=255)
                                                        aaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCtca               <  1:3666103/70‑3 (MQ=255)
                                                          aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat             >  1:3935238/1‑70 (MQ=255)
                                                          aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat             >  1:3079316/1‑70 (MQ=255)
                                                            gagATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAAtcc           <  1:5372620/69‑1 (MQ=255)
                                                            gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc           <  1:560082/71‑1 (MQ=255)
                                                            gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc            <  1:691822/70‑1 (MQ=255)
                                                             cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 >  1:3961053/1‑64 (MQ=255)
                                                              gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                               >  1:2120516/1‑49 (MQ=255)
                                                              gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:5643594/71‑1 (MQ=255)
                                                              gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:5176742/71‑1 (MQ=255)
                                                              gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:2563743/71‑1 (MQ=255)
                                                              gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:1687524/71‑1 (MQ=255)
                                                                  tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                               <  1:698975/45‑1 (MQ=255)
                                                                     aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt  >  1:1436486/1‑71 (MQ=255)
                                                                     |                                                                      
GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  Exported/835721‑835860

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: