Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,388 G→A 33.0% V36M (GTG→ATG)  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,3880GA33.0% 124.5 / 86.1 106V36M (GTG→ATG) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (25/46);  new base A (17/18);  total (42/64)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 108 total aligned reads displayed.

GATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAA  >  Exported/835322‑835454
                                                                  |                                                                  
gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                                 <  1:2055536/70‑1 (MQ=255)
gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                                 >  1:2555624/1‑70 (MQ=255)
gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                                 >  1:5387749/1‑70 (MQ=255)
gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                                 <  1:757757/70‑1 (MQ=255)
gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGCACCAGCCGTGAAGGCGTgg                                                                 <  1:4789397/70‑1 (MQ=255)
  tAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACAAGCCGTGAAGGCGTGGAat                                                              <  1:5430758/71‑1 (MQ=255)
    aGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                             <  1:5479613/70‑1 (MQ=255)
     gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                           <  1:1222241/71‑1 (MQ=255)
      tgtaTCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGa                                                                >  1:1399268/4‑65 (MQ=255)
      tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGc                                                          <  1:2175775/71‑1 (MQ=255)
        taTCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGcat                                                        >  1:4991255/2‑69 (MQ=255)
         aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                       <  1:4955137/71‑1 (MQ=255)
          tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTg                                                       >  1:3219286/1‑70 (MQ=255)
            ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                                                    <  1:1970209/71‑1 (MQ=255)
            ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGt                                                     <  1:2100921/70‑1 (MQ=255)
            ggCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGt                                                     <  1:5213276/70‑1 (MQ=255)
            ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                                                     >  1:2009209/1‑70 (MQ=255)
              cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                  <  1:1467275/71‑1 (MQ=255)
              cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                  <  1:4868144/71‑1 (MQ=255)
                 ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaaca                                               <  1:2643060/71‑1 (MQ=255)
                 ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaaca                                               <  1:814686/71‑1 (MQ=255)
                 ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                <  1:471589/70‑1 (MQ=255)
                  catacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                <  1:4192014/66‑1 (MQ=255)
                  catacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                <  1:3604901/66‑1 (MQ=255)
                  cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                              <  1:3044269/71‑1 (MQ=255)
                  cacaCCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                              <  1:2951129/71‑1 (MQ=255)
                    tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGag                                                               <  1:2068490/51‑2 (MQ=255)
                    tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGag                                                               <  1:170916/51‑2 (MQ=255)
                    cacCTGTTCAGTTACCTGTAAAAACGTCTGGCCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                            <  1:4458463/71‑1 (MQ=255)
                    cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                            >  1:4128877/1‑71 (MQ=255)
                    cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                            <  1:5591519/71‑1 (MQ=255)
                     acCTGTTCGGCGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGttta                                                  <  1:4527608/64‑3 (MQ=255)
                     acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAat                                                              <  1:2963414/52‑1 (MQ=255)
                     acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGt                                           >  1:5084967/1‑71 (MQ=255)
                      ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAgta                                                             >  1:5330937/1‑52 (MQ=255)
                      ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGca                                                         >  1:63110/1‑55 (MQ=255)
                      ccTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGca                                                         >  1:3193917/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)
                         gTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGag                                                               <  1:2268169/47‑2 (MQ=255)
                          ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                      >  1:1289810/1‑71 (MQ=255)
                          ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                      >  1:2403902/1‑71 (MQ=255)
                          ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCg                                                                    >  1:5080133/1‑41 (MQ=255)
                          ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                             >  1:4977295/1‑48 (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:1586002/1‑71 (MQ=255)
                               ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                 >  1:5082360/1‑71 (MQ=255)
                               ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                 >  1:3842573/1‑71 (MQ=255)
                                tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCg                                                                    <  1:871703/35‑1 (MQ=255)
                                gaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaaac                                >  1:3144337/2‑69 (MQ=255)
                                  ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                                   <  1:314561/66‑1 (MQ=255)
                                     gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaacaa                                              >  1:4785181/1‑52 (MQ=255)
                                     gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaacaa                                              >  1:2961760/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:1812442/51‑1 (MQ=255)
                                     gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                               <  1:4768250/51‑1 (MQ=255)
                                     gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                               <  1:1076882/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)
                                      tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                             >  1:5629587/1‑36 (MQ=255)
                                      tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                             >  1:5130005/1‑36 (MQ=255)
                                      tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                        <  1:2172006/57‑1 (MQ=255)
                                      tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                        <  1:2614785/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:3860438/52‑1 (MQ=255)
                                         aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                          <  1:3669423/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:5587888/69‑1 (MQ=255)
                                                  gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                <  1:116334/69‑1 (MQ=255)
                                                       gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac         <  1:2235549/71‑1 (MQ=255)
                                                       gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac         <  1:5170359/71‑1 (MQ=255)
                                                       gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac         <  1:4536268/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)
                                                                 cGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                   <  1:4216180/51‑1 (MQ=255)
                                                                  |                                                                  
GATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAA  >  Exported/835322‑835454

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: