Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,396 T→C 31.2% Y38Y (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,3960TC31.2% 112.8 / 71.8 93Y38Y (TAT→TACnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (23/41);  new base C (15/14);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAA  >  Exported/835327‑835465
                                                                     |                                                                     
gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                                      <  1:1222241/71‑1 (MQ=255)
 tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGc                                                                     <  1:2175775/71‑1 (MQ=255)
   taTCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGcat                                                                   >  1:4991255/2‑69 (MQ=255)
    aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                                  <  1:4955137/71‑1 (MQ=255)
    aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                                  <  1:3262564/71‑1 (MQ=255)
     tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg                                                                 <  1:285264/71‑1 (MQ=255)
     tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                                  >  1:3219286/1‑70 (MQ=255)
       ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                                                               <  1:1970209/71‑1 (MQ=255)
       ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGt                                                                <  1:5213276/70‑1 (MQ=255)
       ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGt                                                                <  1:2100921/70‑1 (MQ=255)
       ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                                                                >  1:2009209/1‑70 (MQ=255)
         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:3044269/71‑1 (MQ=255)
             cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                         <  1:2951129/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)
                acCTGTTCGGCGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGttta                                                             <  1:4527608/64‑3 (MQ=255)
                acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGt                                                      >  1:5084967/1‑71 (MQ=255)
                 ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGca                                                                    >  1:3193917/1‑55 (MQ=255)
                 ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGca                                                                    >  1:63110/1‑55 (MQ=255)
                 ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                     >  1:57590/1‑70 (MQ=255)
                 ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                     >  1:3684700/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:2077286/1‑56 (MQ=255)
                       cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                              >  1:4228091/1‑56 (MQ=255)
                         gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                             <  1:2572197/68‑1 (MQ=255)
                         gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                             <  1:3125305/68‑1 (MQ=255)
                          ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                            >  1:5082360/1‑71 (MQ=255)
                          ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                            >  1:1586002/1‑71 (MQ=255)
                          ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                            >  1:3842573/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:3115841/50‑1 (MQ=255)
                                gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                           <  1:128933/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)
                                 tAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATg                                                                  <  1:353592/42‑1 (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:4597067/1‑61 (MQ=255)
                                   aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                             >  1:2311803/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:2235549/71‑1 (MQ=255)
                                                  gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                    <  1:5170359/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)
                                                                     |                                                                     
GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAA  >  Exported/835327‑835465

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: