Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,506 T→G 12.3% R408R (CGT→CGG 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,5060TG12.3% 197.1 / 22.6 89R408R (CGT→CGGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base G (5/6);  total (42/47)
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.72e-01

TCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAA  >  Exported/836437‑836575
                                                                     |                                                                     
tCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa                                                                      >  1:65595/1‑71 (MQ=255)
 ccGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGAt                                                                     >  1:1989129/1‑71 (MQ=255)
 ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt                                                                     >  1:1161640/1‑71 (MQ=255)
 ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt                                                                     >  1:4710526/1‑71 (MQ=255)
    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)
                         ctgACCGCCGGTTATACCAAACCGGTACTGCGCGCACTGAAACGTa                                                                      >  1:4896566/1‑46 (MQ=255)
                         ctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa                                                                      >  1:5451852/1‑46 (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:3452521/1‑71 (MQ=255)
                                         ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGAc                             >  1:5016305/1‑71 (MQ=255)
                                         ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGAc                             >  1:4891903/1‑71 (MQ=255)
                                           aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCTTTACaaa                                               >  1:1548044/1‑51 (MQ=255)
                                           aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                              <  1:2536256/68‑1 (MQ=255)
                                           aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                              <  1:20819/68‑1 (MQ=255)
                                           aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                              >  1:4140248/1‑68 (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:981980/71‑1 (MQ=255)
                                                  tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                    <  1:1352645/71‑1 (MQ=255)
                                                  tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                    <  1:2113090/71‑1 (MQ=255)
                                                  tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                    <  1:534458/71‑1 (MQ=255)
                                                   aCTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt                    <  1:5254307/70‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                           <  1:671869/46‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                           <  1:5259652/46‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                           <  1:2688951/46‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:5091374/71‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:3764270/71‑1 (MQ=255)
                                                    cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:5137903/71‑1 (MQ=255)
                                                     tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                 <  1:2457603/55‑1 (MQ=255)
                                                     tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                 <  1:4745926/55‑1 (MQ=255)
                                                     tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                                 <  1:4503339/55‑1 (MQ=255)
                                                     tGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                                           <  1:1062416/45‑1 (MQ=255)
                                                      gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACaa                                                >  1:4222580/1‑39 (MQ=255)
                                                      gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                                         >  1:4805507/1‑46 (MQ=255)
                                                      gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                                         >  1:3547417/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:2218809/48‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                                <  1:2356256/48‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                                <  1:2870660/48‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTg                   <  1:2086649/61‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:3262806/62‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:1230990/62‑1 (MQ=255)
                                                             cTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa                  <  1:2798985/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:4124399/1‑62 (MQ=255)
                                                               gAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa                >  1:2508319/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:2676723/1‑53 (MQ=255)
                                                                   cGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGt                     >  1:5248609/1‑53 (MQ=255)
                                                                    gTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGaa  <  1:4576868/71‑1 (MQ=255)
                                                                     |                                                                     
TCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGTGCGCTGGAA  >  Exported/836437‑836575

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: