Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,792 C→G 12.9% S170R (AGC→AGG) ‡ 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,7920CG12.9% 158.1 / 13.0 62S170R (AGC→AGG) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (6/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  Exported/835722‑835860
                                                                      |                                                                    
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:832129/1‑71 (MQ=255)
                tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCt                                                      >  1:2327926/1‑71 (MQ=255)
                  aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc                                                    <  1:1505982/71‑1 (MQ=255)
                    aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:4135534/71‑1 (MQ=255)
                    aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat                                                  <  1:4764202/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:3098997/53‑1 (MQ=255)
                             aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                           <  1:693211/53‑1 (MQ=255)
                             aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         >  1:44376/1‑71 (MQ=255)
                             aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         <  1:3531289/71‑1 (MQ=255)
                             aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                         <  1:945759/71‑1 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                             <  1:91621/49‑1 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       >  1:2123703/1‑71 (MQ=255)
                               aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                       >  1:4756403/1‑71 (MQ=255)
                                     gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                            >  1:2371877/1‑60 (MQ=255)
                                     gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa                                            >  1:2749590/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:1707990/71‑1 (MQ=255)
                                             acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                         <  1:5528516/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:5149149/71‑1 (MQ=255)
                                                     cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 >  1:5462064/1‑71 (MQ=255)
                                                      aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                 >  1:4348567/1‑70 (MQ=255)
                                                      aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:1663632/71‑1 (MQ=255)
                                                      aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:1781324/71‑1 (MQ=255)
                                                      aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                <  1:3186099/71‑1 (MQ=255)
                                                       aaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCtca               <  1:3666103/70‑3 (MQ=255)
                                                         aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat             >  1:3079316/1‑70 (MQ=255)
                                                         aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat             >  1:3935238/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:2563743/71‑1 (MQ=255)
                                                             gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:1687524/71‑1 (MQ=255)
                                                             gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg         <  1:5176742/71‑1 (MQ=255)
                                                                 tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                               <  1:698975/45‑1 (MQ=255)
                                                                    aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt  >  1:1436486/1‑71 (MQ=255)
                                                                      |                                                                    
ACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGT  >  Exported/835722‑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: