Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,279 A→C 15.4% V579V (GTA→GTC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2790AC15.4% 72.8 / 12.5 39V579V (GTA→GTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (13/20);  new base C (2/4);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGT  >  minE/833217‑833345
                                                              |                                                                  
aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAg                                                              >  1:1017707/1‑69 (MQ=255)
  cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                                           >  1:101332/1‑70 (MQ=255)
  cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                                           >  1:2285684/1‑70 (MQ=255)
  cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                                           >  1:1672617/1‑70 (MQ=255)
   aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaa                                                         <  1:2775460/71‑1 (MQ=255)
    aCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaaa                                                        <  1:3908201/71‑1 (MQ=255)
    aCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaaa                                                        <  1:2876246/71‑1 (MQ=255)
     cATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                                       >  1:2136151/1‑71 (MQ=255)
      aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACt                                                      >  1:1388559/1‑71 (MQ=255)
       tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                                                     <  1:2421184/71‑1 (MQ=255)
       tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                                                     >  1:2992050/1‑71 (MQ=255)
         cTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCg                                                   <  1:1008393/71‑1 (MQ=255)
           gATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGt                                                  >  1:3474505/1‑70 (MQ=255)
                tCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                          <  1:3760383/57‑1 (MQ=255)
                 cTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCcaca                                             <  1:780955/69‑1 (MQ=255)
                 cTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                                            <  1:2464834/70‑1 (MQ=255)
                   gCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAAcc                                         >  1:2004451/1‑71 (MQ=255)
                      gttgtgTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGc                                      <  1:2975396/69‑1 (MQ=255)
                      gCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                          <  1:1182187/51‑1 (MQ=255)
                            tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGc                                  <  1:3643274/69‑1 (MQ=255)
                            tCGGGCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGGCCGCAAACCGGCTGGCAg                                <  1:3522510/71‑1 (MQ=255)
                                tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGc                            <  1:852018/71‑1 (MQ=255)
                                tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGc                            <  1:2587079/71‑1 (MQ=255)
                                   gagCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTcc                                                >  1:3593555/1‑48 (MQ=255)
                                   gagCGGGGGCGGCTGGGCGAACTATGTAGGTCAGGAAAAACTGCGTCCGCa                                             <  1:1821906/51‑1 (MQ=255)
                                    agCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCAGCTGGc                        <  1:2408/71‑1 (MQ=255)
                                     gCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGcc                       <  1:3840797/71‑2 (MQ=255)
                                     gCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCg                       <  1:3947740/71‑1 (MQ=255)
                                      cGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaa                                                          <  1:111556/35‑1 (MQ=255)
                                      cGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGt                      <  1:3939621/71‑1 (MQ=255)
                                         gggCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTg                   >  1:1326532/1‑71 (MQ=255)
                                               cTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt              >  1:1619670/1‑70 (MQ=255)
                                               cTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt              >  1:1917754/1‑70 (MQ=255)
                                                tGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCg                       <  1:2454234/60‑1 (MQ=255)
                                                 gggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa                                            <  1:1017521/38‑1 (MQ=255)
                                                  ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc           <  1:2329447/70‑1 (MQ=255)
                                                  ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACt          >  1:1744295/1‑71 (MQ=255)
                                                          aTGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg                                   <  1:2728697/38‑1 (MQ=255)
                                                          aTGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGt  >  1:1583429/1‑71 (MQ=255)
                                                              |                                                                  
ATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGT  >  minE/833217‑833345

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: