Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,003 A→C 28.6% K821Q (AAA→CAA) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0030AC28.6% 45.7 / 40.9 49K821Q (AAA→CAA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base C (4/10);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTA  >  minE/833943‑834072
                                                            |                                                                     
tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTg                                                             >  1:2156777/1‑71 (MQ=255)
 cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                            >  1:681132/1‑71 (MQ=255)
 cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                            >  1:3257707/1‑71 (MQ=255)
  ggATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTgg                                                            <  1:1921612/70‑1 (MQ=255)
   gATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                                                          <  1:1009121/71‑1 (MQ=255)
    atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaa                                                         <  1:110971/71‑1 (MQ=255)
     tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa                                                        <  1:2340387/71‑1 (MQ=255)
     tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa                                                        <  1:2129344/71‑1 (MQ=255)
             cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                                 <  1:429712/70‑1 (MQ=255)
             cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                                 <  1:2271006/70‑1 (MQ=255)
                ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc                                              >  1:583890/1‑70 (MQ=255)
                 tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta                                            >  1:1205943/1‑68 (MQ=255)
                 tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt                                             <  1:1102742/70‑3 (MQ=255)
                   ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc                                          >  1:1281837/1‑71 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt                                             <  1:1713221/63‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                      <  1:921612/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                      <  1:781429/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                      <  1:2048209/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                      <  1:292052/70‑3 (MQ=255)
                         cGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                      >  1:3470978/1‑67 (MQ=255)
                         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                       >  1:1838794/1‑68 (MQ=255)
                         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                       >  1:2037592/1‑68 (MQ=255)
                          gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                                    >  1:1536389/1‑70 (MQ=255)
                             gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc                                                      <  1:2554677/49‑1 (MQ=255)
                              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                       <  1:2575807/63‑1 (MQ=255)
                              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                       <  1:2063502/63‑1 (MQ=255)
                                     cggTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                        <  1:1211621/71‑1 (MQ=255)
                                       gTCGATACGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg                       >  1:210192/1‑70 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAACGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAt                                <  1:1941501/60‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa                                                        <  1:3122900/36‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCg                          <  1:1616871/66‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtg                     >  1:2782389/1‑71 (MQ=255)
                                         cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                       <  1:2342961/52‑1 (MQ=255)
                                           aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaaa                                    <  1:3020136/53‑2 (MQ=255)
                                            tCCGGCCTTGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                 <  1:945227/71‑1 (MQ=255)
                                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                >  1:3054163/1‑71 (MQ=255)
                                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                >  1:1433009/1‑71 (MQ=255)
                                                  cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct                             >  1:3820388/1‑53 (MQ=255)
                                                   ctggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                  >  1:2525992/3‑63 (MQ=255)
                                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc        >  1:1586483/1‑71 (MQ=255)
                                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc        >  1:597351/1‑71 (MQ=255)
                                                      gAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt       <  1:2591893/71‑1 (MQ=255)
                                                       aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg      <  1:1974148/71‑1 (MQ=255)
                                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                        >  1:1634938/1‑52 (MQ=255)
                                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg     >  1:2398472/1‑71 (MQ=255)
                                                         gCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTggg    <  1:1754002/71‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg     <  1:2087617/68‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa  <  1:3754192/71‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa  <  1:1517705/71‑1 (MQ=255)
                                                            |                                                                     
TCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTA  >  minE/833943‑834072

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: