Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,249 T→G 16.2% G569G (GGT→GGG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2490TG16.2% 73.4 / 13.7 43G569G (GGT→GGGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base G (2/5);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  minE/833179‑833319
                                                                      |                                                                      
aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        <  1:532157/71‑1 (MQ=255)
aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        <  1:2039939/71‑1 (MQ=255)
  cACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                                      <  1:2452171/71‑1 (MQ=255)
   aCTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                                                     <  1:229776/71‑1 (MQ=255)
         aTCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                                                               >  1:3192615/1‑71 (MQ=255)
          tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgggg                                                              >  1:3019126/1‑71 (MQ=255)
                cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        >  1:707151/1‑71 (MQ=255)
                 gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                                                        <  1:593856/70‑1 (MQ=255)
                 gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTg                                                       <  1:2358624/71‑1 (MQ=255)
                     tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                                                      >  1:2618477/1‑52 (MQ=255)
                     tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgc                                                    <  1:1425439/70‑1 (MQ=255)
                      gAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggggg                                                             <  1:2117250/60‑1 (MQ=255)
                       aaCTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                                                                  <  1:3379657/54‑1 (MQ=255)
                            tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                            >  1:3229928/1‑71 (MQ=255)
                             cGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg                                            >  1:937595/1‑70 (MQ=255)
                              gTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                                                        >  1:2489336/1‑41 (MQ=255)
                              gTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGt                                           <  1:3169175/70‑1 (MQ=255)
                                 gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgggg                                                              >  1:2546758/1‑48 (MQ=255)
                                      aTCAATATGCTGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAg                                    >  1:1017707/1‑69 (MQ=255)
                                      aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCg                                                                 <  1:235335/40‑1 (MQ=255)
                                      aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggcgg                                                          >  1:2196467/1‑47 (MQ=255)
                                      aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggcgg                                                          >  1:514187/1‑47 (MQ=255)
                                      aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggcgg                                                          >  1:2610281/1‑47 (MQ=255)
                                      aTCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggcgg                                                          >  1:2649595/1‑47 (MQ=255)
                                        cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                 >  1:2285684/1‑70 (MQ=255)
                                        cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                 >  1:101332/1‑70 (MQ=255)
                                        cAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaa                                 >  1:1672617/1‑70 (MQ=255)
                                         aaCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaa                               <  1:2775460/71‑1 (MQ=255)
                                          aCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaaa                              <  1:2876246/71‑1 (MQ=255)
                                          aCATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGaaaaa                              <  1:3908201/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)
                                                tGATCTTCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGgc                                                    <  1:320456/43‑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:2587079/71‑1 (MQ=255)
                                                                      tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGc  <  1:852018/71‑1 (MQ=255)
                                                                      |                                                                      
ACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGC  >  minE/833179‑833319

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: