Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,843 T→C 6.8% D767D (GAT→GAC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,8430TC 6.8% 163.7 / 3.0 59D767D (GAT→GACnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base C (2/2);  total (32/27)
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.61e-01

GGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGA  >  minE/833776‑833913
                                                                   |                                                                      
ggCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtg                                                                     <  1:2231114/71‑1 (MQ=255)
 gCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtgg                                                                    <  1:1579500/71‑1 (MQ=255)
    aCAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                  <  1:3848923/70‑1 (MQ=255)
        ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                             >  1:1042999/1‑71 (MQ=255)
        ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                             >  1:353916/1‑71 (MQ=255)
        ggcggcGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                                                             >  1:93836/1‑71 (MQ=255)
              gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                       >  1:2099251/1‑71 (MQ=255)
              gTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGa                                                       >  1:3922503/1‑71 (MQ=255)
               tGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtgg                                                                    <  1:218400/57‑1 (MQ=255)
               tGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGc                                                           >  1:1420810/1‑66 (MQ=255)
                 aaGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                    >  1:2965366/1‑71 (MQ=255)
                  aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                  >  1:1121658/1‑56 (MQ=255)
                  aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtg                                                                  >  1:1397650/1‑56 (MQ=255)
                  aGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                    <  1:662166/70‑1 (MQ=255)
                      ggAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTc                                               <  1:1820759/71‑1 (MQ=255)
                           aagTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTc                                           <  1:3367957/70‑1 (MQ=255)
                            agTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGAt                                                                        <  1:3753765/40‑1 (MQ=255)
                             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)
                                        ggACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTg                                                         <  1:3636206/43‑1 (MQ=255)
                                        ggACAATGGTCTGGAAGGCAAGCTGGATCTGCTGGTTACGCTg                                                         <  1:1559453/43‑1 (MQ=255)
                                         gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTg                                             <  1:2608539/54‑1 (MQ=255)
                                         gACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTa                            >  1:1941801/1‑71 (MQ=255)
                                          aCAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACtt                                                    <  1:920669/46‑1 (MQ=255)
                                           cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTa                                                              >  1:230671/1‑35 (MQ=255)
                                           cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTa                                                              >  1:1091855/1‑35 (MQ=255)
                                           cAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGtctc                              >  1:3766115/1‑67 (MQ=255)
                                            aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                  >  1:2038522/1‑46 (MQ=255)
                                            aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                  >  1:3195794/1‑46 (MQ=255)
                                            aaTGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTcc                                                  >  1:2844672/1‑46 (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:1254039/1‑47 (MQ=255)
                                                     gAAGGCAAGCTCGACCTGCTGGTGACGCTCGACTTCCGCATGTCCAg                                        >  1:2319323/1‑47 (MQ=255)
                                                      aaGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCa                                      <  1:843662/48‑1 (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)
                                                                   |                                                                      
GGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGA  >  minE/833776‑833913

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: