Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,798 C→T 11.7% F172F (TTC→TTT 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,7980CT11.7% 201.6 / 13.7 77F172F (TTC→TTTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (27/41);  new base T (7/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 3.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCT  >  Exported/835731‑835868
                                                                   |                                                                       
gTGAGTTTGACAAACTGGCGAAATACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa                                                                      <  1:2148505/71‑1 (MQ=255)
 tGAGTTTGACAAACTGGCGAAAGACAATAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa                                                                       <  1:5368627/69‑1 (MQ=255)
  gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                                    >  1:1527763/1‑71 (MQ=255)
  gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                                    >  1:3329932/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:3098997/53‑1 (MQ=255)
                    aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa                                                                    <  1:693211/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: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: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: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:1687524/71‑1 (MQ=255)
                                                    gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                 <  1:5176742/71‑1 (MQ=255)
                                                    gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                 <  1:2563743/71‑1 (MQ=255)
                                                    gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                 <  1:5643594/71‑1 (MQ=255)
                                                        tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                        <  1:698975/45‑1 (MQ=255)
                                                           aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt          >  1:1436486/1‑71 (MQ=255)
                                                               aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc       >  1:2746859/1‑70 (MQ=255)
                                                                gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGa     >  1:4997169/1‑71 (MQ=255)
                                                                gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa               <  1:370947/61‑1 (MQ=255)
                                                                gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa               <  1:2105302/61‑1 (MQ=255)
                                                                gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg      <  1:1012600/70‑1 (MQ=255)
                                                                gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg      <  1:4811165/70‑1 (MQ=255)
                                                                gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg      <  1:4585217/70‑1 (MQ=255)
                                                                 ttCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                        <  1:1620059/36‑1 (MQ=255)
                                                                 ttCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                        <  1:1471726/36‑1 (MQ=255)
                                                                 ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGc                <  1:4370896/59‑1 (MQ=255)
                                                                  tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACt                              >  1:2540103/1‑45 (MQ=255)
                                                                  tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                       <  1:5537493/52‑1 (MQ=255)
                                                                  tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                       <  1:149834/52‑1 (MQ=255)
                                                                  tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                       <  1:3212023/52‑1 (MQ=255)
                                                                  tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg           >  1:2505592/1‑63 (MQ=255)
                                                                   cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg           >  1:5626662/1‑62 (MQ=255)
                                                                   cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgt            <  1:2833587/61‑1 (MQ=255)
                                                                   cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt  >  1:4604916/1‑71 (MQ=255)
                                                                   cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt  >  1:139118/1‑71 (MQ=255)
                                                                   cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt  >  1:83831/1‑71 (MQ=255)
                                                                   |                                                                       
GTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCT  >  Exported/835731‑835868

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: