Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,737 G→T 10.3% P1065P (CCG→CCT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,7370GT10.3% 177.1 / 5.6 58P1065P (CCG→CCTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (17/34);  new base T (4/2);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAAC  >  minE/834667‑834806
                                                                      |                                                                     
gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                                       <  1:432688/71‑1 (MQ=255)
gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                                       <  1:3034171/71‑1 (MQ=255)
    gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg                                                                   <  1:1606177/71‑1 (MQ=255)
       aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca                                                                >  1:817567/1‑71 (MQ=255)
       aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                                  <  1:3585401/69‑1 (MQ=255)
        cTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGacac                                                               >  1:838356/1‑71 (MQ=255)
         tGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc                                                                    <  1:3039403/65‑1 (MQ=255)
          gTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAccc                                                             >  1:719501/1‑71 (MQ=255)
            aTCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCg                                                            >  1:1853893/1‑70 (MQ=255)
             tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCg                                                            >  1:1281944/1‑69 (MQ=255)
             tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGta                                                          <  1:3874277/71‑2 (MQ=255)
             tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGta                                                          <  1:309197/71‑2 (MQ=255)
             tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                                          <  1:3686967/71‑1 (MQ=255)
               aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTc                                                         <  1:426008/70‑1 (MQ=255)
               aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTc                                                         <  1:984172/70‑1 (MQ=255)
                ggATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGtag                                                         >  1:2337017/1‑67 (MQ=255)
                ggATCAACCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCg                                                                   >  1:192/1‑59 (MQ=255)
                 gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt                                                      <  1:2205433/71‑1 (MQ=255)
                 gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt                                                      <  1:2155285/71‑1 (MQ=255)
                   tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                                    <  1:3705689/71‑1 (MQ=255)
                      ccAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGac                                                                 <  1:1158622/55‑1 (MQ=255)
                        aGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGacac                                                               >  1:1098755/1‑55 (MQ=255)
                          tGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                                      >  1:3342206/1‑46 (MQ=255)
                          tGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa                                                                      >  1:1193428/1‑46 (MQ=255)
                             aTGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt                                          <  1:3639100/71‑1 (MQ=255)
                               gCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                                       <  1:2224679/40‑1 (MQ=255)
                               gCGTGATTTCGGTGAAAGCCTGCTGGTGTATCGTccgccg                                                                       <  1:3564993/40‑1 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca                                                                <  1:3750205/44‑1 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                                  >  1:2025409/1‑42 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                                  >  1:1284508/1‑42 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                                          <  1:2024037/50‑1 (MQ=255)
                                  tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                                          <  1:852786/50‑1 (MQ=255)
                                    aTTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAg                                   <  1:2938012/71‑1 (MQ=255)
                                       tCGGTGAAATCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaa                                 <  1:1241974/70‑1 (MQ=255)
                                       tCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa                                >  1:772528/1‑71 (MQ=255)
                                        cGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa                                <  1:707921/70‑1 (MQ=255)
                                          gTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTc                                                         <  1:55652/43‑1 (MQ=255)
                                           tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                                    <  1:2454449/47‑1 (MQ=255)
                                           tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                                    <  1:1183517/47‑1 (MQ=255)
                                            gAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGaa                                                   >  1:3474268/1‑47 (MQ=255)
                                            gAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaagaa                                               >  1:2420037/1‑51 (MQ=255)
                                             aaaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCgg                                                       >  1:1152171/1‑42 (MQ=255)
                                             aaaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                             <  1:948919/52‑1 (MQ=255)
                                              aaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACg                         >  1:3408497/1‑71 (MQ=255)
                                               aGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                             <  1:1688354/50‑1 (MQ=255)
                                                gCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt                                             <  1:2708890/49‑1 (MQ=255)
                                                 ccTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                            >  1:633183/1‑49 (MQ=255)
                                                 ccTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                            >  1:3106249/1‑49 (MQ=255)
                                                     ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAcc                   <  1:2243051/70‑1 (MQ=255)
                                                     ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc                    >  1:2115524/1‑69 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa             <  1:695058/70‑1 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa             <  1:1714168/70‑1 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa             <  1:1709648/70‑1 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa             <  1:382810/70‑1 (MQ=255)
                                                           tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa             <  1:1290459/70‑1 (MQ=255)
                                                            aTCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATc                              <  1:2097804/52‑1 (MQ=255)
                                                             tCGTCCGCCCATCGACACCCGTTCGGTGAAAGAAGTGAt                                          >  1:2184057/1‑39 (MQ=255)
                                                               gTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt                                          <  1:1420284/37‑1 (MQ=255)
                                                               gTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACgg                        >  1:2147811/1‑55 (MQ=255)
                                                               gTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgcg        <  1:878384/71‑1 (MQ=255)
                                                               gTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgc         <  1:2004777/70‑1 (MQ=255)
                                                                tCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAgcgc       <  1:2947793/71‑1 (MQ=255)
                                                                  cgccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCa                                    <  1:3922878/40‑1 (MQ=255)
                                                                  cgccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACTGCAACCAGGaaaaa           >  1:1522981/1‑65 (MQ=255)
                                                                  cgccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTc     <  1:1617917/71‑1 (MQ=255)
                                                                  cgccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTc     <  1:3466315/71‑1 (MQ=255)
                                                                   gccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaa                                 <  1:3651736/42‑1 (MQ=255)
                                                                   gccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGa                                  >  1:1391066/1‑41 (MQ=255)
                                                                   gccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGa                                  >  1:3234374/1‑41 (MQ=255)
                                                                   gccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAg          >  1:3314647/1‑65 (MQ=255)
                                                                    ccgATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCaa   >  1:247537/1‑71 (MQ=255)
                                                                     cgATCGACACACGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAAc  >  1:3405719/1‑71 (MQ=255)
                                                                      |                                                                     
GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAAC  >  minE/834667‑834806

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: