Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,724 T→C 10.4% V1061A (GTT→GCT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,7240TC10.4% 175.0 / 10.9 67V1061A (GTT→GCT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base C (5/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACC  >  minE/834655‑834789
                                                                     |                                                                 
cGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGtt                                                                  <  1:1766230/71‑1 (MQ=255)
cGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGt                                                                   <  1:1069052/70‑1 (MQ=255)
 gCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGctt                                                                 <  1:1594662/71‑1 (MQ=255)
 gCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGtt                                                                  >  1:3119623/1‑70 (MQ=255)
  ctctctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTa                                                                <  1:211357/71‑1 (MQ=255)
      ctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGt                                                            >  1:372095/1‑71 (MQ=255)
       tGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                           >  1:1663510/1‑71 (MQ=255)
            gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                      <  1:3034171/71‑1 (MQ=255)
            gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg                                                      <  1:432688/71‑1 (MQ=255)
               aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                       >  1:1179273/1‑67 (MQ=255)
               aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                       >  1:3103029/1‑67 (MQ=255)
                gCAACTGTATCAGGATCACCAGTGGATTCGTGATTTCGGTGAAAGCCTGCTGGTTTATCg                                                             <  1:2242873/60‑1 (MQ=255)
                gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg                                                  <  1:1606177/71‑1 (MQ=255)
                 cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCGTGCTGGTTTa                                                                >  1:3890062/1‑56 (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)
                         tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                           >  1:576257/1‑53 (MQ=255)
                         tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                                                           >  1:580911/1‑53 (MQ=255)
                         tCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGtt                                         <  1:3686967/71‑1 (MQ=255)
                           aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATc                                                              >  1:44742/1‑48 (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)
                             gATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCgcc                                                           >  1:3401231/1‑47 (MQ=255)
                             gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt                                     <  1:2205433/71‑1 (MQ=255)
                             gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt                                     <  1:2155285/71‑1 (MQ=255)
                               tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccg                                                         <  1:1893049/49‑1 (MQ=255)
                               tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccg                                                         <  1:3398889/49‑1 (MQ=255)
                               tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa                                   <  1:3705689/71‑1 (MQ=255)
                                cACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGttt                                                                 <  1:2154747/40‑1 (MQ=255)
                                 aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                       >  1:2883648/1‑49 (MQ=255)
                                 aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                       >  1:2888313/1‑49 (MQ=255)
                                 aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                                       >  1:3233194/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:1284508/1‑42 (MQ=255)
                                              tGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa                                                 >  1:2025409/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: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)
                                                                     |                                                                 
CGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACC  >  minE/834655‑834789

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: