Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 836,494 C→G 21.3% R404R (CGC→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,4940CG21.3% 158.8 / 38.4 81R404R (CGC→CGGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base G (7/10);  total (35/46)
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.29e-01

AGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAC  >  Exported/836426‑836562
                                                                    |                                                                    
aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                                                                     >  1:1105263/1‑70 (MQ=255)
   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:4710526/1‑71 (MQ=255)
            ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt                                                        >  1:1161640/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:603611/1‑52 (MQ=255)
                    gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                                   >  1:2672883/1‑52 (MQ=255)
                    gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                                                   >  1:4812093/1‑52 (MQ=255)
                      aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               >  1:236465/1‑54 (MQ=255)
                      aTCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGATg                                              >  1:295549/1‑71 (MQ=255)
                         tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:1838486/51‑1 (MQ=255)
                         tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGa                                                               <  1:1812489/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: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:20819/68‑1 (MQ=255)
                                                      aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                 <  1:2536256/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:534458/71‑1 (MQ=255)
                                                             tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt       <  1:981980/71‑1 (MQ=255)
                                                             tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt       <  1:2113090/71‑1 (MQ=255)
                                                             tACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt       <  1:1352645/71‑1 (MQ=255)
                                                              aCTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt       <  1:5254307/70‑1 (MQ=255)
                                                               cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                              <  1:2688951/46‑1 (MQ=255)
                                                               cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                              <  1:5259652/46‑1 (MQ=255)
                                                               cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg                              <  1:671869/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:4503339/55‑1 (MQ=255)
                                                                tGCGGGCGCTGAAACGGATGATGGCGATGCGCCACTATATGCGTTCACAAACCGt                    <  1:2457603/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: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)
                                                                    |                                                                    
AGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAC  >  Exported/836426‑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: