Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,873 G→C 15.1% S777S (TCG→TCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8730GC15.1% 159.1 / 18.1 66S777S (TCG→TCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (6/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833805‑833943
                                                                    |                                                                      
gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                                      <  1:1867046/71‑1 (MQ=255)
gTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAg                                                                      >  1:170765/1‑71 (MQ=255)
   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:517880/1‑70 (MQ=255)
                   gTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGa                                                    >  1:628468/1‑70 (MQ=255)
                    tCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCg                                                                        <  1:1354309/49‑1 (MQ=255)
                     tttgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGa                                                    <  1:1693588/65‑1 (MQ=255)
                     attgAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAGTACATGCCTGTTCTCCGatat                                                 <  1:78509/68‑1 (MQ=255)
                      tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                                 >  1:3906359/1‑70 (MQ=255)
                      tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                                 >  1:1794463/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)
                        gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                                                      >  1:2319323/1‑47 (MQ=255)
                        gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                                                      >  1:1254039/1‑47 (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:2964215/1‑70 (MQ=255)
                                            ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                           >  1:2571116/1‑70 (MQ=255)
                                             gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                                     >  1:871411/1‑43 (MQ=255)
                                             gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCg                                                     >  1:2511607/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:3568828/37‑1 (MQ=255)
                                                         ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:3573381/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)
                                                                    |                                                                      
GTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833805‑833943

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: