Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,482 T→G 17.1% G980G (GGT→GGG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,4820TG17.1% 75.9 / 19.0 47G980G (GGT→GGGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (12/27);  new base G (4/4);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAG  >  minE/834417‑834551
                                                                 |                                                                     
cTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                                                                    <  1:2422971/69‑1 (MQ=255)
     aCCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATa                                                             >  1:779164/1‑71 (MQ=255)
       cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATAc                                                            <  1:3126901/70‑1 (MQ=255)
        ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGt                                                                       <  1:493808/58‑1 (MQ=255)
        ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGc                                                          <  1:2654637/71‑1 (MQ=255)
        ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGc                                                          <  1:2134806/71‑1 (MQ=255)
         gAAACCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCa                                                         >  1:1739912/1‑71 (MQ=255)
           aaCCAACGGTCATGTTGCAGTCAAAGCGTGGCAGGCGCTGGGCGAGATCACCGGGCGCGAACATACCCATc                                                       >  1:2476558/1‑71 (MQ=255)
              cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTg                                                                   <  1:1609248/56‑1 (MQ=255)
              cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATa                                                             >  1:1111526/1‑62 (MQ=255)
                aCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                                  <  1:98880/71‑1 (MQ=255)
                 cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                                                                    <  1:1740392/52‑1 (MQ=255)
                    tCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCt                                                <  1:1918581/69‑1 (MQ=255)
                       ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATa                                           <  1:3896094/71‑1 (MQ=255)
                        gTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATa                                           >  1:1021293/1‑70 (MQ=255)
                          aGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATa                                           >  1:2346364/1‑68 (MQ=255)
                           gCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAaag                                        <  1:2065168/70‑1 (MQ=255)
                            ccGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCa                                                               <  1:2867127/46‑1 (MQ=255)
                              gTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATaa                                          >  1:2907642/1‑65 (MQ=255)
                              gTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAagaaga                                    <  1:2445255/71‑1 (MQ=255)
                               tGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTgg                                                    <  1:2459053/54‑1 (MQ=255)
                               tGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAAt                                            <  1:3687186/62‑1 (MQ=255)
                                 aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATa                                                             >  1:1644837/1‑43 (MQ=255)
                                 aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAAt                                            >  1:834706/1‑60 (MQ=255)
                                 aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACg                                  <  1:1527250/70‑1 (MQ=255)
                                      cTGGGCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATAAAGAAGACgag                                >  1:800475/1‑67 (MQ=255)
                                       tGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCt                                                      <  1:863779/44‑1 (MQ=255)
                                              cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGAt                           >  1:2332549/1‑64 (MQ=255)
                                              cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                     <  1:2870825/70‑1 (MQ=255)
                                              cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                     <  1:1227682/70‑1 (MQ=255)
                                              cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCt                      <  1:3683722/69‑1 (MQ=255)
                                               ccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                     <  1:485273/69‑1 (MQ=255)
                                                 tGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcg                  >  1:243645/1‑70 (MQ=255)
                                                   aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa               >  1:288234/1‑71 (MQ=255)
                                                    gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa               >  1:210116/1‑70 (MQ=255)
                                                     cGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat              <  1:2344688/70‑1 (MQ=255)
                                                       aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAAt                                            <  1:2719780/38‑1 (MQ=255)
                                                          tcaCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAg        <  1:544803/69‑1 (MQ=255)
                                                           caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg       >  1:833000/2‑71 (MQ=255)
                                                           caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg       >  1:1978356/2‑71 (MQ=255)
                                                             ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt           <  1:3300136/65‑1 (MQ=255)
                                                             ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCt                      <  1:1319235/54‑1 (MQ=255)
                                                             ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg     <  1:525403/71‑2 (MQ=255)
                                                             ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg     <  1:1850715/71‑2 (MQ=255)
                                                              cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat            <  1:3829395/63‑1 (MQ=255)
                                                              cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac    <  1:708986/71‑1 (MQ=255)
                                                                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                   <  1:3301163/38‑1 (MQ=255)
                                                                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAg  <  1:3711096/71‑1 (MQ=255)
                                                                 |                                                                     
CTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAG  >  minE/834417‑834551

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: