Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,876 C→T 12.5% S778S (AGC→AGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8760CT12.5% 168.6 / 13.0 64S778S (AGC→AGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (4/4);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833808‑833943
                                                                    |                                                                   
gACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc                                                                   <  1:353209/71‑1 (MQ=255)
    ggCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGt                                                               >  1:2150872/1‑71 (MQ=255)
     gCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTg                                                                <  1:3879027/69‑1 (MQ=255)
     gCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtc                                                              <  1:2423677/71‑1 (MQ=255)
      cAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtct                                                             <  1:841087/71‑1 (MQ=255)
         gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                                          >  1:1941801/1‑71 (MQ=255)
           cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                                            >  1:3766115/1‑67 (MQ=255)
                gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                                    >  1:628468/1‑70 (MQ=255)
                gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                                    >  1:517880/1‑70 (MQ=255)
                  tttgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGa                                                    <  1:1693588/65‑1 (MQ=255)
                  attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACATGCCTGTTCTCCGatat                                                 <  1:78509/68‑1 (MQ=255)
                   tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                                 >  1:1794463/1‑70 (MQ=255)
                   tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                                 >  1:3906359/1‑70 (MQ=255)
                     gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCAc                                                                   >  1:562351/1‑50 (MQ=255)
                     gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:2543467/70‑1 (MQ=255)
                     gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcataat                                              >  1:1805447/1‑71 (MQ=255)
                      aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCa                                                                    <  1:843662/48‑1 (MQ=255)
                      aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc                                                       >  1:3540995/1‑61 (MQ=255)
                       aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                            <  1:2103152/71‑1 (MQ=255)
                       aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                            <  1:2738845/71‑1 (MQ=255)
                       aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtt                                             >  1:2705971/1‑70 (MQ=255)
                        ggCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtttt                                           <  1:3019309/71‑1 (MQ=255)
                         gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                                            <  1:142174/53‑1 (MQ=255)
                         gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                                     <  1:2905222/60‑1 (MQ=255)
                         gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg                                          >  1:465022/1‑71 (MQ=255)
                          cAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                                          >  1:337403/1‑54 (MQ=255)
                            aGCTGGATCTGGTGGTTACCCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                                                          >  1:2885293/1‑52 (MQ=255)
                             gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa                                      >  1:996399/1‑71 (MQ=255)
                              cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc                                     <  1:2875369/71‑1 (MQ=255)
                                   tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt                                                         <  1:476899/46‑1 (MQ=255)
                                   tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGa                                >  1:1051907/1‑71 (MQ=255)
                                    cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                               >  1:1257131/1‑71 (MQ=255)
                                     tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                                                            <  1:2172931/41‑1 (MQ=255)
                                     tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc                                     >  1:2466527/1‑64 (MQ=255)
                                       gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                             <  1:3694872/70‑1 (MQ=255)
                                       gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                            >  1:1699652/1‑71 (MQ=255)
                                        tggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                           <  1:2478699/71‑1 (MQ=255)
                                         ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                           >  1:2571116/1‑70 (MQ=255)
                                         ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                           >  1:2964215/1‑70 (MQ=255)
                                          gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                                     >  1:2511607/1‑43 (MQ=255)
                                          gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                                     >  1:871411/1‑43 (MQ=255)
                                          gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGt                          <  1:1800876/70‑1 (MQ=255)
                                             aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTAt                                                         >  1:260950/1‑36 (MQ=255)
                                             aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                            >  1:1969392/1‑65 (MQ=255)
                                             aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                            >  1:2002814/1‑65 (MQ=255)
                                              cGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaa                     <  1:1824569/71‑1 (MQ=255)
                                               gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAc                                                   <  1:704334/40‑1 (MQ=255)
                                               gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa                    <  1:1861217/71‑1 (MQ=255)
                                               gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                    <  1:2948002/71‑1 (MQ=255)
                                                 tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                    >  1:3614129/3‑69 (MQ=255)
                                                 tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa                     >  1:1667589/3‑68 (MQ=255)
                                                 tGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                               >  1:1817918/1‑58 (MQ=255)
                                                  cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGACACCTGGTACGaaa                    >  1:2271443/2‑68 (MQ=255)
                                                   gACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaat                 <  1:1768595/70‑2 (MQ=255)
                                                   gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga                >  1:1895780/1‑71 (MQ=255)
                                                      ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:3573381/37‑1 (MQ=255)
                                                      ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:3568828/37‑1 (MQ=255)
                                                      ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                             <  1:1625578/55‑1 (MQ=255)
                                                      ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                             <  1:2947024/55‑1 (MQ=255)
                                                             tGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                    <  1:2073247/57‑1 (MQ=255)
                                                               tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa       <  1:3617365/68‑1 (MQ=255)
                                                               tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATAc    >  1:1785461/1‑71 (MQ=255)
                                                                 gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACtt  >  1:1507433/1‑71 (MQ=255)
                                                                  aGCACCTGTCTCTATTCCGACATCATTTTGCCGACgg                                   >  1:2317658/1‑37 (MQ=255)
                                                                    |                                                                   
GACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833808‑833943

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: