Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,218 T→C 23.4% D892D (GAT→GAC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,2180TC23.4% 110.6 / 31.4 64D892D (GAT→GACnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (6/9);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAG  >  minE/834150‑834288
                                                                    |                                                                      
tGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTa                                                                      >  1:2027666/1‑71 (MQ=255)
    gcaaaGGCGAATGCGATCTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTAccct                                                                  >  1:2952907/3‑70 (MQ=255)
     aaaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCgg                                                                 >  1:220715/1‑71 (MQ=255)
       aaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTAt                                                                     <  1:5422594/65‑1 (MQ=255)
         ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGa                                                              >  1:1320910/1‑70 (MQ=255)
          gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt                                                            <  1:2612768/71‑1 (MQ=255)
          gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt                                                            <  1:2308008/71‑1 (MQ=255)
             aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                         <  1:1073281/71‑1 (MQ=255)
             aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                         <  1:4876532/71‑1 (MQ=255)
              gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTa                                                          >  1:5092015/1‑69 (MQ=255)
              gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                         >  1:1805810/1‑70 (MQ=255)
                   atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                     <  1:1171900/67‑1 (MQ=255)
                   atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                     <  1:1531459/67‑1 (MQ=255)
                   atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                     <  1:1759713/67‑1 (MQ=255)
                      tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATcc                                                                   <  1:230778/52‑1 (MQ=255)
                      tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc                                                     <  1:3983531/66‑1 (MQ=255)
                        aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTATAGCGCGATTATCCGGCGACTTACGAACGCtt                                                 <  1:1439717/68‑1 (MQ=255)
                        aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCtt                                                 <  1:1762731/68‑1 (MQ=255)
                        aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:1664951/1‑71 (MQ=255)
                        aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:3137129/1‑71 (MQ=255)
                            cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                          <  1:2961830/71‑1 (MQ=255)
                            cAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTacc                                                                    >  1:3652708/1‑43 (MQ=255)
                             aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                         >  1:4131261/1‑71 (MQ=255)
                             aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGc                                                   <  1:2274748/61‑1 (MQ=255)
                             aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACctcgc                                         <  1:1793489/71‑3 (MQ=255)
                              ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                      >  1:4401786/1‑57 (MQ=255)
                               gTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCttt                                                <  1:3205612/62‑1 (MQ=255)
                                 aaaaCCTCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                     >  1:3373877/1‑71 (MQ=255)
                                 aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                     >  1:4059875/1‑71 (MQ=255)
                                 aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                     >  1:3264816/1‑71 (MQ=255)
                                 aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                     <  1:537943/71‑1 (MQ=255)
                                 aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                     >  1:2028724/1‑71 (MQ=255)
                                       gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:2475642/1‑56 (MQ=255)
                                       gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:2682170/1‑56 (MQ=255)
                                       gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                               <  1:151129/71‑1 (MQ=255)
                                           cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:1043999/1‑52 (MQ=255)
                                             cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctc                                           <  1:5323215/53‑1 (MQ=255)
                                              acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                          <  1:38355/53‑1 (MQ=255)
                                              acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                           >  1:2989700/1‑68 (MQ=255)
                                               caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                         >  1:3775983/1‑53 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                      >  1:4049950/1‑39 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                      >  1:2117296/1‑39 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                      >  1:2688210/1‑39 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGa                                                       >  1:1524131/1‑38 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc                                                     <  1:4409047/40‑1 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:1337985/1‑47 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:4700031/1‑47 (MQ=255)
                                                aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                              >  1:4276509/1‑47 (MQ=255)
                                                 ttATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaa                     <  1:276779/71‑1 (MQ=255)
                                                    tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt                                        >  1:1417611/1‑49 (MQ=255)
                                                      gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                <  1:512139/71‑1 (MQ=255)
                                                      gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                <  1:2863443/71‑1 (MQ=255)
                                                       tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGGGAAAATCGGt               <  1:918393/71‑1 (MQ=255)
                                                       tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt               <  1:4618116/71‑1 (MQ=255)
                                                       tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt               <  1:2213556/71‑1 (MQ=255)
                                                         gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa             >  1:2423121/1‑71 (MQ=255)
                                                           agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg           >  1:2292926/1‑71 (MQ=255)
                                                           agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg           >  1:1718139/1‑71 (MQ=255)
                                                            gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg          <  1:2770130/71‑1 (MQ=255)
                                                            gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg          <  1:2331647/71‑1 (MQ=255)
                                                             agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc         <  1:2732854/71‑1 (MQ=255)
                                                               cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt                   <  1:2180099/59‑1 (MQ=255)
                                                                 cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt      >  1:2191912/1‑70 (MQ=255)
                                                                   aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                <  1:3378964/58‑1 (MQ=255)
                                                                    ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGc                                    >  1:3829941/1‑37 (MQ=255)
                                                                    ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa             >  1:5053776/1‑60 (MQ=255)
                                                                    ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt      >  1:2175296/1‑67 (MQ=255)
                                                                    ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAg  >  1:4055910/1‑71 (MQ=255)
                                                                    |                                                                      
TGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAG  >  minE/834150‑834288

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: