Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 833,972 G→T 48.6% L810L (CTG→CTT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported833,9720GT48.6% 33.3 / 84.7 74L810L (CTG→CTTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (12/26);  new base T (12/24);  total (24/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAG  >  Exported/833903‑834038
                                                                     |                                                                  
gCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                                    <  1:3934286/70‑1 (MQ=255)
     cGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcg                                                              >  1:182574/1‑71 (MQ=255)
       gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                            <  1:4778080/71‑1 (MQ=255)
       gCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                            <  1:2071345/71‑1 (MQ=255)
        cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                                    <  1:1610348/62‑1 (MQ=255)
            cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTc                                                        >  1:3764463/1‑70 (MQ=255)
            cTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg                                                       <  1:495935/71‑1 (MQ=255)
               gTACGAAAAAGACTATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                                      <  1:2753154/69‑1 (MQ=255)
               gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                                      <  1:774894/69‑1 (MQ=255)
               gTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATc                                                    >  1:5188932/1‑71 (MQ=255)
               gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCTGTCGAt                                                     <  1:5336978/70‑1 (MQ=255)
               gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                                     <  1:3132333/70‑1 (MQ=255)
               gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                                     <  1:2241915/70‑1 (MQ=255)
               gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                                     <  1:4590972/70‑1 (MQ=255)
                 aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg                                                  >  1:93619/1‑71 (MQ=255)
                    aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc                                               >  1:159557/1‑71 (MQ=255)
                     aaaaGACGATATGAACATCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                              <  1:5196116/71‑3 (MQ=255)
                     aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                              <  1:1165135/71‑3 (MQ=255)
                     aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc                                               <  1:2971644/70‑1 (MQ=255)
                      aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg                                             <  1:527492/71‑1 (MQ=255)
                          acgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggc                                                            <  1:4214709/52‑1 (MQ=255)
                            gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                        <  1:5605469/70‑1 (MQ=255)
                            gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                        <  1:2757281/70‑1 (MQ=255)
                            gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                        <  1:673802/70‑1 (MQ=255)
                             ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                             <  1:4542792/46‑1 (MQ=255)
                             ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                             <  1:4017036/46‑1 (MQ=255)
                             ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                             <  1:3953829/46‑1 (MQ=255)
                             ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcgg                                                             <  1:3226024/46‑1 (MQ=255)
                             acATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                                      <  1:295201/55‑1 (MQ=255)
                              cATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTg                                                                    >  1:4451928/1‑40 (MQ=255)
                               aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg                                                           >  1:761349/1‑48 (MQ=255)
                               aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcggcg                                                           >  1:1377363/1‑48 (MQ=255)
                               aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTgcagcg                                                           >  1:1760693/1‑48 (MQ=255)
                               aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCg                                                       <  1:52347/52‑1 (MQ=37)
                               aTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGa                                                      <  1:3625764/53‑1 (MQ=37)
                                tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc                                                    <  1:2116227/54‑1 (MQ=255)
                                tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                     <  1:2070046/69‑1 (MQ=255)
                                 gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa                                    >  1:4454508/1‑69 (MQ=255)
                                  aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                        >  1:1726496/1‑64 (MQ=255)
                                   aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                <  1:3035158/71‑1 (MQ=255)
                                      cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                                      <  1:642867/46‑1 (MQ=255)
                                      cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                                      <  1:317352/46‑1 (MQ=255)
                                      cTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGa                                                      <  1:2661546/46‑1 (MQ=255)
                                       ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa                                         <  1:2503919/57‑1 (MQ=255)
                                       ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCAGCGTGGGaa                                         <  1:2231905/57‑1 (MQ=255)
                                        tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATc                                                    >  1:3815887/1‑46 (MQ=255)
                                         cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                          >  1:5187583/1‑71 (MQ=255)
                                          ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg                         <  1:1846389/71‑1 (MQ=255)
                                          ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                          >  1:3689816/1‑70 (MQ=255)
                                             tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa                      <  1:1481207/71‑1 (MQ=255)
                                             tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa                      <  1:2915590/71‑1 (MQ=255)
                                             tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                        <  1:3322576/69‑1 (MQ=255)
                                               tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCgg                                                 <  1:3744128/42‑1 (MQ=38)
                                               tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                           <  1:1551834/48‑1 (MQ=37)
                                               tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                           <  1:4323523/48‑1 (MQ=37)
                                               tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg                                           <  1:2517544/48‑1 (MQ=37)
                                               tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGc                                       <  1:2124322/52‑1 (MQ=255)
                                                gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa                              >  1:2439283/1‑60 (MQ=255)
                                                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt                   <  1:3401193/70‑1 (MQ=255)
                                                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                  <  1:4377260/71‑1 (MQ=255)
                                                 cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                  <  1:771430/71‑1 (MQ=255)
                                                   tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                <  1:78578/55‑1 (MQ=255)
                                                     cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa              <  1:1469837/71‑1 (MQ=255)
                                                     cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa                >  1:4287315/1‑69 (MQ=255)
                                                     cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                              >  1:2225820/1‑39 (MQ=255)
                                                     cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa                >  1:3522696/1‑69 (MQ=255)
                                                      gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt                     <  1:4112380/63‑1 (MQ=255)
                                                       tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                     <  1:2210751/46‑1 (MQ=255)
                                                       tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                <  1:5515961/51‑1 (MQ=255)
                                                        ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg            >  1:4407565/1‑69 (MQ=255)
                                                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                          >  1:5369519/1‑55 (MQ=255)
                                                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                          >  1:5665279/1‑55 (MQ=255)
                                                         tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta          >  1:4662936/1‑68 (MQ=255)
                                                         tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa          >  1:5064259/1‑71 (MQ=255)
                                                         tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa          >  1:3288003/1‑71 (MQ=255)
                                                               tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca    <  1:5645485/70‑3 (MQ=255)
                                                                  gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag  >  1:4232274/1‑70 (MQ=255)
                                                                     |                                                                  
GCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAG  >  Exported/833903‑834038

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: