Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,369 T→C 10.7% T609T (ACT→ACC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,3690TC10.7% 148.6 / 9.8 56T609T (ACT→ACCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (21/29);  new base C (4/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACC  >  minE/833299‑833437
                                                                      |                                                                    
cGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAct                                                                      >  1:3167630/1‑71 (MQ=255)
   aaaCCGGCTGTCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                   >  1:3302308/1‑71 (MQ=255)
   aaaCCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCAcctc                                                                    >  1:2172934/1‑70 (MQ=255)
   aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                   <  1:1320267/71‑1 (MQ=255)
     aCCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCAcctcgtt                                                                 >  1:1776245/1‑68 (MQ=255)
       cGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAtt                                                               <  1:152361/71‑1 (MQ=255)
         gCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAttta                                                             >  1:2680760/1‑70 (MQ=255)
             gCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                          <  1:1040045/70‑1 (MQ=255)
              ttgCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACaa                                                        <  1:2131717/69‑1 (MQ=255)
                  cGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATaa                                                        >  1:2773946/1‑67 (MQ=255)
                   gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACt                                                   >  1:775826/1‑71 (MQ=255)
                     tGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTc                                                  <  1:3227087/70‑1 (MQ=255)
                       gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                               >  1:3170020/4‑71 (MQ=255)
                       gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                               >  1:229897/4‑71 (MQ=255)
                        cGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagc                                              <  1:254711/71‑1 (MQ=255)
                          tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                          <  1:613701/57‑1 (MQ=255)
                          tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                          <  1:2309997/57‑1 (MQ=255)
                            tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAt                                                                >  1:597055/1‑49 (MQ=255)
                            tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTcca                                                >  1:3372621/1‑65 (MQ=255)
                            tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                          >  1:1590333/1‑71 (MQ=255)
                              cccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCa                                                     <  1:184800/58‑1 (MQ=255)
                               ccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCta                                                           >  1:188770/1‑51 (MQ=255)
                                attGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTa                                                                 >  1:999774/2‑44 (MQ=255)
                                 ttGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                          <  1:832556/50‑1 (MQ=255)
                                  tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                   >  1:2874469/1‑40 (MQ=255)
                                  tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                                   >  1:2879497/1‑40 (MQ=255)
                                   gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                                          >  1:22360/1‑48 (MQ=255)
                                   gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGtt                                    <  1:1186222/70‑1 (MQ=255)
                                   gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa                                   >  1:3239879/1‑71 (MQ=255)
                                        gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                              >  1:617901/1‑71 (MQ=255)
                                        gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                              >  1:2589941/1‑71 (MQ=255)
                                           gCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                          <  1:1024171/56‑1 (MQ=255)
                                            cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt                          <  1:1574199/71‑1 (MQ=255)
                                            cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt                          <  1:3222984/71‑1 (MQ=255)
                                             gTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTc                         <  1:2144655/71‑1 (MQ=255)
                                                cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                          >  1:2523068/1‑51 (MQ=255)
                                                cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                                          >  1:685868/1‑51 (MQ=255)
                                                  gcgcGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggc                    <  1:3085662/71‑1 (MQ=255)
                                                   cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc                                             <  1:529854/45‑1 (MQ=255)
                                                   cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                           >  1:1769060/1‑60 (MQ=255)
                                                   cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                           >  1:2353032/1‑60 (MQ=255)
                                                    gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg                  <  1:1408708/71‑1 (MQ=255)
                                                       tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa                               <  1:2789547/55‑1 (MQ=255)
                                                       tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg                  >  1:4022660/1‑68 (MQ=255)
                                                       tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg                      <  1:593061/64‑1 (MQ=255)
                                                       tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag               <  1:3654739/71‑1 (MQ=255)
                                                        cacaTGAACAGCACTTCTTATTTCTATAACCACTccagcc                                             <  1:1225315/40‑1 (MQ=255)
                                                         acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                                 >  1:4035666/1‑35 (MQ=255)
                                                         acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                                 >  1:1666908/1‑35 (MQ=255)
                                                         acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag             <  1:3111514/71‑1 (MQ=255)
                                                           aTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAaga              <  1:3941187/68‑1 (MQ=255)
                                                            tGAACAGCACTTCTTATTTCTATAACCACTccagc                                              <  1:2052576/35‑1 (MQ=255)
                                                            tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg                      <  1:4062324/59‑1 (MQ=255)
                                                            tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGaa                <  1:251319/65‑1 (MQ=255)
                                                            tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg          >  1:2453286/1‑71 (MQ=255)
                                                              aaCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa                                   >  1:1258885/1‑44 (MQ=255)
                                                               aCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa                               <  1:2891241/47‑1 (MQ=255)
                                                                 aGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc     <  1:314178/71‑1 (MQ=255)
                                                                  gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc     <  1:481431/70‑1 (MQ=255)
                                                                  gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCa    <  1:2750492/71‑1 (MQ=255)
                                                                   cACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg          >  1:2171872/1‑64 (MQ=255)
                                                                   cACTTCATATTTCTATAACCACTCCAGCCAGTGGCGTTa                                   >  1:536322/1‑39 (MQ=255)
                                                                    aCTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCAcc  <  1:3387854/71‑1 (MQ=255)
                                                                      |                                                                    
CGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACC  >  minE/833299‑833437

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: