Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,789 T→C 12.5% Y169Y (TAT→TAC 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,7890TC12.5% 157.3 / 12.7 64Y169Y (TAT→TACnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (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.31e-01

GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGG  >  Exported/835721‑835853
                                                                    |                                                                
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:1527763/1‑71 (MQ=255)
            gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                    >  1:3329932/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:3531289/71‑1 (MQ=255)
                              aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                  <  1:945759/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: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: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)
                                                                    |                                                                
GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGG  >  Exported/835721‑835853

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