Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,405 C→T 22.0% F41F (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,4050CT22.0% 110.6 / 26.0 59F41F (TTC→TTTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base T (10/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 3.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATA  >  Exported/835336‑835467
                                                                     |                                                              
cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                               <  1:4868144/71‑1 (MQ=255)
cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                               <  1:1467275/71‑1 (MQ=255)
   ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaaca                                                            <  1:814686/71‑1 (MQ=255)
   ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaaca                                                            <  1:2643060/71‑1 (MQ=255)
   ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                             <  1:471589/70‑1 (MQ=255)
    catacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                             <  1:3604901/66‑1 (MQ=255)
    catacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                             <  1:4192014/66‑1 (MQ=255)
    cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                           <  1:2951129/71‑1 (MQ=255)
    cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                           <  1:3044269/71‑1 (MQ=255)
    cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                           <  1:2106645/71‑1 (MQ=255)
      cacCTGTTCAGTTACCTGTAAAAACGTCTGGCCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                                         <  1:4458463/71‑1 (MQ=255)
      cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                                         <  1:5591519/71‑1 (MQ=255)
      cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                                         >  1:4128877/1‑71 (MQ=255)
       acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGt                                                        >  1:5084967/1‑71 (MQ=255)
        ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                       >  1:3684700/1‑70 (MQ=255)
        ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                       >  1:57590/1‑70 (MQ=255)
         cTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACTCATGGTTTAACAACGtc                                                       >  1:1973877/1‑69 (MQ=255)
            ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                   >  1:2403902/1‑71 (MQ=255)
            ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                   >  1:1289810/1‑71 (MQ=255)
            ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                     <  1:1053184/69‑1 (MQ=255)
              cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                                >  1:4228091/1‑56 (MQ=255)
              cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                                >  1:2077286/1‑56 (MQ=255)
                gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                               <  1:3125305/68‑1 (MQ=255)
                gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                               <  1:2572197/68‑1 (MQ=255)
                 ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                              >  1:1586002/1‑71 (MQ=255)
                 ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                              >  1:3842573/1‑71 (MQ=255)
                 ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                              >  1:5082360/1‑71 (MQ=255)
                  gaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaaac                                             >  1:3144337/2‑69 (MQ=255)
                   aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                             >  1:1615938/1‑54 (MQ=255)
                    ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                                                <  1:314561/66‑1 (MQ=255)
                      tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                         >  1:5021954/1‑71 (MQ=255)
                       gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaacaa                                                           >  1:2961760/1‑52 (MQ=255)
                       gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaacaa                                                           >  1:4785181/1‑52 (MQ=255)
                       gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                             <  1:128933/50‑1 (MQ=255)
                       gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                             <  1:3115841/50‑1 (MQ=255)
                       gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCggg                                         >  1:3579405/1‑70 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                            <  1:4768250/51‑1 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                            <  1:1076882/51‑1 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                            <  1:1812442/51‑1 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                     >  1:2549367/1‑58 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                     >  1:1070581/1‑58 (MQ=255)
                       gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                         >  1:65983/1‑70 (MQ=255)
                        tAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCggg                                         <  1:5453935/69‑1 (MQ=255)
                        tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                     <  1:2614785/57‑1 (MQ=255)
                        tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                     <  1:2172006/57‑1 (MQ=255)
                        tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                        <  1:5391119/70‑1 (MQ=255)
                         aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac                                                          >  1:3419345/1‑51 (MQ=255)
                         aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                        >  1:763680/1‑69 (MQ=255)
                          aaaaCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGt                                                        >  1:2524200/1‑52 (MQ=255)
                          aaaaCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCGGGcatt                                     >  1:3831433/1‑68 (MQ=255)
                          aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                             <  1:1976516/47‑1 (MQ=255)
                          aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                               >  1:2311803/1‑61 (MQ=255)
                          aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                               >  1:4597067/1‑61 (MQ=255)
                          aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAg                                     <  1:5291826/71‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                       <  1:3669423/52‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                       <  1:3860438/52‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                       <  1:3225637/52‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                       <  1:2925466/52‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCGGGGCCa                                      <  1:7934/69‑1 (MQ=255)
                           aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                                      <  1:1014431/69‑1 (MQ=255)
                             aCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                   >  1:1029114/1‑54 (MQ=255)
                               gTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                >  1:1468629/1‑71 (MQ=255)
                                  tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                             >  1:2553376/1‑71 (MQ=255)
                                    gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                             <  1:116334/69‑1 (MQ=255)
                                    gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                             <  1:5587888/69‑1 (MQ=255)
                                         gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                      <  1:4536268/71‑1 (MQ=255)
                                         gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                      <  1:5170359/71‑1 (MQ=255)
                                         gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                      <  1:2235549/71‑1 (MQ=255)
                                          ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                     <  1:2013259/71‑1 (MQ=255)
                                            gTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTgg                   <  1:1987430/71‑1 (MQ=255)
                                               aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGa                 <  1:591039/70‑1 (MQ=255)
                                                aGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa               <  1:4420474/71‑1 (MQ=255)
                                                   cGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                <  1:4216180/51‑1 (MQ=255)
                                                     tGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                                             <  1:3736977/36‑1 (MQ=255)
                                                     tGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                                       <  1:2714899/42‑1 (MQ=255)
                                                      ggAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAg          >  1:719327/1‑70 (MQ=255)
                                                       gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                             >  1:2942501/1‑50 (MQ=255)
                                                       gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTcc                              >  1:3956612/1‑49 (MQ=255)
                                                            tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaaa    <  1:165773/70‑1 (MQ=255)
                                                             gCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa  >  1:2757157/1‑71 (MQ=255)
                                                                     |                                                              
CTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATA  >  Exported/835336‑835467

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