Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,497 A→G 20.2% A405A (GCA→GCG 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,4970AG20.2% 163.6 / 36.5 84A405A (GCA→GCGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (7/10);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAC  >  Exported/836429‑836562
                                                                    |                                                                 
cTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCt                                                                 <  1:2254244/71‑1 (MQ=255)
cTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACt                                                                 <  1:2687358/71‑1 (MQ=255)
   cGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGa                                                               <  1:4432431/70‑1 (MQ=255)
   cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGaa                                                              <  1:4110428/71‑1 (MQ=255)
   cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               >  1:3630906/1‑70 (MQ=255)
    gTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGaa                                                              >  1:1837881/1‑70 (MQ=255)
     tATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCt                                                                 <  1:1188187/66‑1 (MQ=255)
       ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg                                                           >  1:401318/1‑70 (MQ=255)
       ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg                                                           <  1:1596027/70‑1 (MQ=255)
        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)
                 gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                                   >  1:4812093/1‑52 (MQ=255)
                 gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                                   >  1:603611/1‑52 (MQ=255)
                 gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                                   >  1:2672883/1‑52 (MQ=255)
                   aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               >  1:236465/1‑54 (MQ=255)
                   aTCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGATg                                              >  1:295549/1‑71 (MQ=255)
                      tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:1812489/51‑1 (MQ=255)
                      tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:1838486/51‑1 (MQ=255)
                      tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:727777/51‑1 (MQ=255)
                           aaTCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGaaa                                                             >  1:5265714/1‑48 (MQ=255)
                           aaTCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:1351984/46‑1 (MQ=255)
                               tgctgACCGCCGGTGATACCAAACCGGTACTGCTCGCACTGa                                                               <  1:4878885/42‑1 (MQ=255)
                               tgctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:4524877/42‑1 (MQ=255)
                               tgctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:2862653/42‑1 (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:2774006/1‑62 (MQ=255)
                                  tgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCAt                                        >  1:1198279/1‑62 (MQ=255)
                                      cgccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCt                            >  1:3694536/1‑70 (MQ=255)
                                        ccgGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                                           <  1:1461203/37‑1 (MQ=255)
                                        ccgGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                                           <  1:4073761/37‑1 (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:1352645/71‑1 (MQ=255)
                                                          tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt       <  1:534458/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:2688951/46‑1 (MQ=255)
                                                            cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                              <  1:671869/46‑1 (MQ=255)
                                                            cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                              <  1:5259652/46‑1 (MQ=255)
                                                            cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa     <  1:5137903/71‑1 (MQ=255)
                                                            cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa     <  1:3764270/71‑1 (MQ=255)
                                                            cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa     <  1:5091374/71‑1 (MQ=255)
                                                             tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                    <  1:4745926/55‑1 (MQ=255)
                                                             tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                    <  1:2457603/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: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)
                                                                    |                                                                 
CTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAC  >  Exported/836429‑836562

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: