Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,270 G→A 14.6% A576A (GCG→GCA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2700GA14.6% 99.9 / 8.8 41A576A (GCG→GCAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base A (2/4);  total (17/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.28e-01

TCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACT  >  minE/833207‑833337
                                                               |                                                                   
tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                                              >  1:3229928/1‑71 (MQ=255)
 cGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                                              >  1:937595/1‑70 (MQ=255)
  gTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGt                                                             <  1:3169175/70‑1 (MQ=255)
          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)
               cATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGt                                                             >  1:14721/1‑57 (MQ=255)
               cATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAAc                                               >  1:2136151/1‑71 (MQ=255)
                aTGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACt                                              >  1:1388559/1‑71 (MQ=255)
                 tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                                             >  1:2992050/1‑71 (MQ=255)
                 tGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTg                                             <  1:2421184/71‑1 (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:1917754/1‑70 (MQ=255)
                                                         cTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt      >  1:1619670/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)
                                                               |                                                                   
TCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACT  >  minE/833207‑833337

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: