Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,882 T→C 14.3% C780C (TGT→TGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8820TC14.3% 136.5 / 13.6 56C780C (TGT→TGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (4/4);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833812‑833943
                                                                      |                                                             
ggCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGt                                                               >  1:2150872/1‑71 (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:3906359/1‑70 (MQ=255)
               tGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcat                                                 >  1:1794463/1‑70 (MQ=255)
                 gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:2543467/70‑1 (MQ=255)
                 gAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGAcataat                                              >  1:1805447/1‑71 (MQ=255)
                  aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTc                                                       >  1:3540995/1‑61 (MQ=255)
                   aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                            <  1:2738845/71‑1 (MQ=255)
                   aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                            <  1:2103152/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: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:3573381/37‑1 (MQ=255)
                                                  ttCCGTCTGTCGAGCACCTGTCTCTATTCCGAcatca                                               <  1:3568828/37‑1 (MQ=255)
                                                  ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                             <  1:2947024/55‑1 (MQ=255)
                                                  ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                             <  1:1625578/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)
                                                                    tGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATg         <  1:2630184/57‑1 (MQ=255)
                                                                     gTCTCTATTCCGACATCATTTTGCCGACGGCGAcc                              <  1:1438584/35‑1 (MQ=255)
                                                                      |                                                             
GGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTT  >  minE/833812‑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: