Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,731 T→C 8.3% R1063R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,7310TC 8.3% 196.6 / 7.9 72R1063R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (26/40);  new base C (4/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGC  >  minE/834661‑834801
                                                                      |                                                                      
ctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGt                                                                        >  1:372095/1‑71 (MQ=255)
 tGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                                       >  1:1663510/1‑71 (MQ=255)
      gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                                  <  1:432688/71‑1 (MQ=255)
      gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                                  <  1:3034171/71‑1 (MQ=255)
         aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                                   >  1:1179273/1‑67 (MQ=255)
         aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                                   >  1:3103029/1‑67 (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:309197/71‑2 (MQ=255)
                   tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGta                                                     <  1:3874277/71‑2 (MQ=255)
                   tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                                       >  1:576257/1‑53 (MQ=255)
                   tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                                       >  1:580911/1‑53 (MQ=255)
                   tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                                     <  1:3686967/71‑1 (MQ=255)
                     aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTc                                                    <  1:984172/70‑1 (MQ=255)
                     aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTc                                                    <  1:426008/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)
                         tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccg                                                                     <  1:3398889/49‑1 (MQ=255)
                         tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccg                                                                     <  1:1893049/49‑1 (MQ=255)
                         tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                               <  1:3705689/71‑1 (MQ=255)
                           aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                                   >  1:3233194/1‑49 (MQ=255)
                           aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                                   >  1:2883648/1‑49 (MQ=255)
                           aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                                   >  1:2888313/1‑49 (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)
                                   aTGCGTGATTTCGGTGAAAGCGTGCTGGTTTATCGTccgcc                                                                   <  1:871147/41‑1 (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:852786/50‑1 (MQ=255)
                                        tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                                     <  1:2024037/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:1183517/47‑1 (MQ=255)
                                                 tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                               <  1:2454449/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:3106249/1‑49 (MQ=255)
                                                       ccTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg                                       >  1:633183/1‑49 (MQ=255)
                                                           ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAcc              <  1:2243051/70‑1 (MQ=255)
                                                           ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc               >  1:2115524/1‑69 (MQ=255)
                                                                 tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa        <  1:1290459/70‑1 (MQ=255)
                                                                 tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa        <  1:1709648/70‑1 (MQ=255)
                                                                 tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa        <  1:1714168/70‑1 (MQ=255)
                                                                 tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa        <  1:695058/70‑1 (MQ=255)
                                                                 tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa        <  1:382810/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)
                                                                      |                                                                      
CTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGC  >  minE/834661‑834801

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: