Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,966 C→T 34.2% H808H (CAC→CAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,9660CT34.2% 38.4 / 30.9 38H808H (CAC→CATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (5/8);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATC  >  minE/833897‑834032
                                                                     |                                                                  
catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAc                                                                    >  1:1547707/1‑70 (MQ=255)
          aCGGCGACCTGGTACGAAAAAGACGACATGAATACTTAGGATATGCATCCGTTTATTCACCCGCTGTCTgc                                                         >  1:184044/1‑71 (MQ=255)
              cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggcg                                                     <  1:1316571/71‑1 (MQ=255)
                aCCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGt                                                   <  1:3340782/71‑1 (MQ=255)
                 ccTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTc                                                  >  1:1082343/1‑71 (MQ=255)
                 ccTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                  >  1:523411/1‑71 (MQ=255)
                 ccTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                  >  1:645282/1‑71 (MQ=255)
                  cTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAc                                                                    <  1:2702202/52‑1 (MQ=255)
                  cTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCg                                                 <  1:1147388/71‑1 (MQ=255)
                   tGGTACGAAAAACACGACATGAATACTTCGGATATGCATCCGTTTATTCAcc                                                                   <  1:382861/52‑1 (MQ=255)
                    ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt                                               >  1:539908/1‑71 (MQ=255)
                     gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                               <  1:2894652/70‑1 (MQ=255)
                       aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg                                            >  1:1174350/1‑71 (MQ=255)
                          aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCg                                                                 <  1:3899757/47‑1 (MQ=255)
                           aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                        <  1:2898435/71‑3 (MQ=255)
                           aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                        >  1:3946166/1‑71 (MQ=255)
                           aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                        <  1:2567963/71‑1 (MQ=255)
                            aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg                                       <  1:3714343/71‑1 (MQ=255)
                            aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                  >  1:1896657/1‑60 (MQ=255)
                            aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg                                       <  1:942542/71‑1 (MQ=255)
                            aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg                                       >  1:2116689/1‑71 (MQ=255)
                                 cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa                                   <  1:3126815/70‑1 (MQ=255)
                                 cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                  <  1:404217/71‑1 (MQ=255)
                                   ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa                                   <  1:2987246/66‑1 (MQ=255)
                                         aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTgg                                      <  1:1288601/59‑1 (MQ=255)
                                           aCTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGa                                    <  1:3113469/59‑1 (MQ=255)
                                              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:2271006/70‑1 (MQ=255)
                                                           cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa         <  1:429712/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)
                                                               tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa                                   >  1:2600892/1‑40 (MQ=255)
                                                                 ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc  >  1:1281837/1‑71 (MQ=255)
                                                                     |                                                                  
CATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATC  >  minE/833897‑834032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: