Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,510 C→A 10.8% L410M (CTG→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
*Exported836,5100CA10.8% 284.1 / 13.4 92L410M (CTG→ATG) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (35/47);  new base A (4/6);  total (39/53)
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.36e-01

GTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAAGAG  >  Exported/836441‑836578
                                                                     |                                                                    
gTACAGTATCTGGCGAATCTGCTGACCGCCGGAGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt                                                                     <  1:537934/71‑1 (MQ=255)
  aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt                                                                     <  1:1362446/69‑1 (MQ=255)
   caatATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGc                                                                  >  1:982427/4‑71 (MQ=255)
    aatATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGc                                                                  >  1:5160913/3‑70 (MQ=255)
       aTCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGATg                                                              >  1:295549/1‑71 (MQ=255)
                     ctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAAc                                                 <  1:3179718/70‑1 (MQ=255)
                     ctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                                                >  1:563775/1‑71 (MQ=255)
                      tgACCGCCGGTGATACCAAACCGGTACTGCGTGCACTGAAACGTATGCTGGCGATGCGTCAt                                                        >  1:2935420/1‑62 (MQ=255)
                      tgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCAt                                                        >  1:1198279/1‑62 (MQ=255)
                      tgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCAt                                                        >  1:2774006/1‑62 (MQ=255)
                          cgccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                                            >  1:3694536/1‑70 (MQ=255)
                            ccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTg                                                                    <  1:1610253/44‑1 (MQ=255)
                             cgGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGATGCGCCact                                                       <  1:1111407/56‑3 (MQ=255)
                               gTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAcc                                      >  1:5142085/1‑71 (MQ=255)
                                tGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTCGCGATGCGTCATTACAAACGTGCTGAAACCg                                     <  1:4679175/71‑1 (MQ=255)
                                 gATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGt                                    >  1:3463473/1‑71 (MQ=255)
                                 gATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGt                                    >  1:2490242/1‑71 (MQ=255)
                                  aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGGGATGCGTCATTACAAACGTGCTGAAACCGtt                                   <  1:4286664/71‑1 (MQ=255)
                                    aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGa                                                                >  1:1714345/1‑40 (MQ=255)
                                    aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                                 >  1:200883/1‑71 (MQ=255)
                                     ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGAc                                >  1:4891903/1‑71 (MQ=255)
                                     ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGAc                                >  1:5016305/1‑71 (MQ=255)
                                     ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGAc                                >  1:3452521/1‑71 (MQ=255)
                                       aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCTTTACaaa                                                  >  1:1548044/1‑51 (MQ=255)
                                       aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                                 <  1:2536256/68‑1 (MQ=255)
                                       aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                                 >  1:4140248/1‑68 (MQ=255)
                                       aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                                 <  1:20819/68‑1 (MQ=255)
                                        gtccGGTACTGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGTGGAAggc                             >  1:3326425/3‑70 (MQ=255)
                                        gtccGGTACTGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGTGGAAgg                              <  1:2408245/68‑1 (MQ=255)
                                         aCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGt                             >  1:4052406/1‑70 (MQ=255)
                                         aCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTa                            <  1:1353678/71‑1 (MQ=255)
                                          ccGGTACTGCGCTCACTTAAACGTATGCTGGCGATGCGTGATTACa                                                    <  1:1580587/46‑1 (MQ=38)
                                          ccGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTaa                           <  1:4350344/71‑1 (MQ=255)
                                          ccGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTa                            >  1:5233940/1‑70 (MQ=255)
                                              tACTGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTa                                                      <  1:2596377/40‑1 (MQ=12)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTAc                                                     <  1:4777029/41‑1 (MQ=255)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTAc                                                     <  1:3921070/41‑1 (MQ=255)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:534458/71‑1 (MQ=255)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:1352645/71‑1 (MQ=255)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:981980/71‑1 (MQ=255)
                                              tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:2113090/71‑1 (MQ=255)
                                               aCTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:5254307/70‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                              <  1:5259652/46‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                              <  1:2688951/46‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                              <  1:671869/46‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:5091374/71‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:5137903/71‑1 (MQ=255)
                                                cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:3764270/71‑1 (MQ=255)
                                                 tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                    <  1:2457603/55‑1 (MQ=255)
                                                 tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                    <  1:4503339/55‑1 (MQ=255)
                                                 tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                    <  1:4745926/55‑1 (MQ=255)
                                                 tGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                              <  1:1062416/45‑1 (MQ=255)
                                                  gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACaa                                                   >  1:4222580/1‑39 (MQ=255)
                                                  gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                                            >  1:3547417/1‑46 (MQ=255)
                                                  gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                                            >  1:4805507/1‑46 (MQ=255)
                                                  gcgcgcAATGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa                   >  1:5299741/1‑71 (MQ=255)
                                                   cgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa                   >  1:4197050/1‑70 (MQ=255)
                                                   cgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAc                  >  1:2717866/1‑71 (MQ=255)
                                                         cTGAAACGTATGCTGTCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTg                      <  1:3156322/61‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                                   <  1:2870660/48‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                                   <  1:2356256/48‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                                   <  1:2218809/48‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTg                      <  1:2086649/61‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:3262806/62‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:2798985/62‑1 (MQ=255)
                                                         cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:1230990/62‑1 (MQ=255)
                                                          tGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTTAAACCGTTGACGGTAAAGtt                       <  1:5454375/59‑1 (MQ=255)
                                                          tGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:3389093/59‑1 (MQ=255)
                                                          tGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                       <  1:2103654/59‑1 (MQ=255)
                                                          tGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTgcg           >  1:1802740/1‑71 (MQ=255)
                                                           gAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa                   >  1:2508319/1‑62 (MQ=255)
                                                           gAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa                   >  1:4124399/1‑62 (MQ=255)
                                                           gAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTgcgc          >  1:4077020/1‑71 (MQ=255)
                                                           gAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTgcgc          >  1:2254745/1‑71 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTaa                           <  1:4062642/53‑1 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTaa                           <  1:3286745/53‑1 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:1762796/59‑1 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     >  1:5472843/1‑59 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCt         >  1:5257954/1‑71 (MQ=255)
                                                            aaaCGTATGCTGGCGATGCGTAATTACAAACGTGCTGAAACCGTTGACGGTaa                           <  1:3407170/53‑1 (MQ=255)
                                                               cGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGt                        >  1:5248609/1‑53 (MQ=255)
                                                               cGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGt                        >  1:2676723/1‑53 (MQ=255)
                                                                gTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa     <  1:4576868/71‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:1677112/52‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:2508483/52‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                     <  1:3795849/52‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa     <  1:1645623/68‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa     <  1:5386644/68‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa     >  1:3141503/1‑68 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa     <  1:4291632/68‑1 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAagag  >  1:619894/1‑71 (MQ=255)
                                                                   tGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAagag  >  1:2160778/1‑71 (MQ=255)
                                                                     cTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAg                         >  1:2451687/1‑46 (MQ=255)
                                                                     |                                                                    
GTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAAGAG  >  Exported/836441‑836578

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: