Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,485 T→C 17.4% R981R (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,4850TC17.4% 99.5 / 15.2 46R981R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (13/25);  new base C (4/4);  total (17/29)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCC  >  minE/834417‑834553
                                                                    |                                                                    
cTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGt                                                                      <  1:2422971/69‑1 (MQ=255)
     aCCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATa                                                               >  1:779164/1‑71 (MQ=255)
       cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATAc                                                              <  1:3126901/70‑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)
                                                                  cGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGcc  >  1:1505855/1‑71 (MQ=255)
                                                                    |                                                                    
CTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCC  >  minE/834417‑834553

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: