Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,029 C→T 19.2% A829A (GCC→GCT) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0290CT19.2% 109.6 / 19.7 52A829A (GCC→GCT) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base T (4/6);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

TTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTA  >  minE/833962‑834099
                                                                   |                                                                      
ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc                                                                     >  1:1281837/1‑71 (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:2037592/1‑68 (MQ=255)
      cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                                  >  1:1838794/1‑68 (MQ=255)
       gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                                                               >  1:1536389/1‑70 (MQ=255)
           tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                                  <  1:2063502/63‑1 (MQ=255)
           tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                                                  <  1:2575807/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)
                     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:597351/1‑71 (MQ=255)
                                  ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                                   >  1:1586483/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:1517705/71‑1 (MQ=255)
                                        gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                             <  1:3754192/71‑1 (MQ=255)
                                          aaaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaaa                           <  1:3922874/71‑1 (MQ=255)
                                            aGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCaaa                           >  1:2475190/1‑69 (MQ=255)
                                             gCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGa                         <  1:3540844/70‑1 (MQ=255)
                                                aCTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGAAAcc                     <  1:880791/71‑1 (MQ=255)
                                                 cTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                             <  1:2596003/62‑1 (MQ=255)
                                                  tGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGaaa                       >  1:3043177/1‑67 (MQ=255)
                                                    ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGAc                  >  1:1845737/1‑70 (MQ=255)
                                                    ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa                 >  1:2880556/1‑71 (MQ=255)
                                                    ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa                 <  1:716772/71‑1 (MQ=255)
                                                    ggAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACa                 >  1:378323/1‑71 (MQ=255)
                                                     gAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGCGTGGGTCATCTTGGCAAAGAAACCGAcg                 >  1:3090338/1‑69 (MQ=255)
                                                      aaaTCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATc               <  1:3558644/71‑1 (MQ=255)
                                                       aaTCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCg              <  1:2328367/71‑1 (MQ=255)
                                                          cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc                                   >  1:335780/1‑47 (MQ=255)
                                                           tACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                       >  1:2993752/1‑58 (MQ=255)
                                                            aCAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa                             <  1:40990/51‑1 (MQ=255)
                                                            aCAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGaaa                       >  1:1087941/1‑57 (MQ=255)
                                                               taGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt                                  <  1:354189/42‑1 (MQ=255)
                                                               aaGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCt       <  1:1274669/70‑1 (MQ=255)
                                                                   cATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa  >  1:3232936/1‑71 (MQ=255)
                                                                   |                                                                      
TTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTA  >  minE/833962‑834099

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: