Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,488 C→A 17.3% D982E (GAC→GAA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,4880CA17.3% 122.8 / 10.8 46D982E (GAC→GAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (14/24);  new base A (4/4);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCG  >  minE/834422‑834556
                                                                  |                                                                    
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)
         cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATa                                                                  >  1:1111526/1‑62 (MQ=255)
           aCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgcg                                                       <  1:98880/71‑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)
                                                                 aCCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgcg  >  1:3373149/1‑70 (MQ=255)
                                                                  |                                                                    
ACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCG  >  minE/834422‑834556

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: