Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,214 C→A 10.7% L558M (CTG→ATG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2140CA10.7% 286.5 / 12.8 93L558M (CTG→ATG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base A (5/5);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 102 total aligned reads displayed.

GTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGT  >  minE/833146‑833282
                                                                    |                                                                    
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                      <  1:611356/69‑1 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                      <  1:349528/69‑1 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                      <  1:1454271/69‑1 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                      <  1:5667402/69‑1 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCt                                                                     >  1:1401083/1‑70 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                                                                    <  1:1420865/71‑1 (MQ=255)
gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                                                                    <  1:5611212/71‑1 (MQ=255)
 tCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGa                                                                   >  1:4611107/1‑71 (MQ=255)
 tCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGa                                                                   >  1:5689090/1‑71 (MQ=255)
   gTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGAt                                                                  >  1:1571150/1‑70 (MQ=255)
   gTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATc                                                                 <  1:2996689/71‑1 (MQ=255)
   gTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATc                                                                 <  1:2989894/71‑1 (MQ=255)
    ttCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCa                                                                <  1:1583912/71‑1 (MQ=255)
    ttCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCa                                                                <  1:5547149/71‑1 (MQ=255)
      cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                               >  1:1259031/1‑70 (MQ=255)
      cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                               >  1:5430147/1‑70 (MQ=255)
      cGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCaa                                                               >  1:2968102/1‑70 (MQ=255)
       gatgatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCaa                                                               <  1:5499255/69‑1 (MQ=255)
        atgatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAatat                                                            <  1:4954499/71‑1 (MQ=255)
         tgatTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                               <  1:5172200/67‑1 (MQ=255)
           atTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGc                                                          <  1:3981034/70‑1 (MQ=255)
            ttATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCt                                                         >  1:4704277/1‑70 (MQ=255)
               tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATc                                                                 >  1:1731914/1‑59 (MQ=255)
               tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCa                                                                >  1:5692257/1‑60 (MQ=255)
               tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtt                                                     >  1:5324934/1‑71 (MQ=255)
               tcgtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAt                                                      <  1:4684714/70‑1 (MQ=255)
                 gtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                      <  1:3251542/52‑1 (MQ=255)
                       gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                   >  1:4236686/1‑65 (MQ=255)
                       gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                   >  1:2960422/1‑65 (MQ=255)
                       gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCgg                                             >  1:5051978/1‑71 (MQ=255)
                       gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg                                              >  1:4316618/1‑70 (MQ=255)
                       gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg                                              >  1:5672695/1‑70 (MQ=255)
                          gggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                <  1:5569653/65‑1 (MQ=255)
                          gggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                <  1:1363383/65‑1 (MQ=255)
                           ggCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCtgt                                         <  1:1035629/71‑1 (MQ=255)
                            gCTGAACCACTGGTATGACCTCGATATGAACTATCGTGGTc                                                                      <  1:1450204/41‑1 (MQ=255)
                            gCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCtgtg                                        >  1:3725412/1‑71 (MQ=255)
                              ttaaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTgtgt                                       <  1:3337644/68‑1 (MQ=255)
                               gAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCg                                     >  1:3202465/1‑71 (MQ=255)
                                aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACAt                                                            <  1:3254791/47‑1 (MQ=255)
                                aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACAt                                                            <  1:4789883/47‑1 (MQ=255)
                                aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACAt                                                            <  1:5208398/47‑1 (MQ=255)
                                aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  <  1:4140865/57‑1 (MQ=255)
                                aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  <  1:4707401/57‑1 (MQ=255)
                                aaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgg                                    >  1:2572011/1‑71 (MQ=255)
                                 aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                               <  1:2665655/43‑1 (MQ=255)
                                 aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                   <  1:575325/71‑1 (MQ=255)
                                 aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                   <  1:1917723/71‑1 (MQ=255)
                                 aCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCggg                                   >  1:3094736/1‑70 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGc                                                          >  1:2178360/1‑47 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGc                                                          >  1:1368912/1‑47 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                >  1:4676367/1‑57 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                >  1:2864839/1‑57 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                                   >  1:3930466/1‑70 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTc                                  <  1:1400640/71‑1 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGGGCTGTGTCGGt                                   >  1:5052111/1‑69 (MQ=255)
                                  ccACTGGTATCACCTCGATATGAACTATCGGGGTCTGATCAACATGCTGATTTTCTg                                                >  1:5203024/1‑57 (MQ=255)
                                  ccACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgggc                                  <  1:1205448/71‑3 (MQ=255)
                                  ccACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCgggc                                  <  1:4169593/71‑3 (MQ=255)
                                   cACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                 <  1:1693458/71‑1 (MQ=255)
                                    aCTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                <  1:3824254/71‑1 (MQ=255)
                                    aCTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTGCGGTTGTGTCGGGCa                                 >  1:5614519/1‑70 (MQ=255)
                                        gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                              <  1:763299/69‑1 (MQ=255)
                                           tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  <  1:4202694/46‑1 (MQ=255)
                                           tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  <  1:5452341/46‑1 (MQ=255)
                                           tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                             <  1:3514023/67‑1 (MQ=255)
                                           tCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCCGCTGTGTCGGTCAGAGc                             <  1:577276/67‑1 (MQ=255)
                                           tCACCTCGATATGAACTATCGTGATCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                             <  1:3016432/67‑1 (MQ=255)
                                            cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                            <  1:1675534/51‑1 (MQ=255)
                                            cACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGc                                            <  1:5107298/51‑1 (MQ=255)
                                             aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                   >  1:2173763/1‑43 (MQ=255)
                                             aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTg                                                >  1:5588664/1‑46 (MQ=255)
                                             aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                           >  1:1604126/1‑67 (MQ=255)
                                                 cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt                   >  1:832257/1‑71 (MQ=255)
                                                  gATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                                   >  1:3695424/1‑38 (MQ=255)
                                                  gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  >  1:5221555/1‑39 (MQ=255)
                                                  gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                                  >  1:884748/1‑39 (MQ=255)
                                                  gATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTg                  >  1:1200263/1‑71 (MQ=255)
                                                   atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                 <  1:267822/55‑1 (MQ=255)
                                                   atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCgg                           >  1:3879120/1‑61 (MQ=255)
                                                   atatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGAgg                           >  1:3948894/1‑61 (MQ=255)
                                                    tatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTggg                <  1:4653097/71‑1 (MQ=255)
                                                      tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                <  1:2083719/53‑1 (MQ=255)
                                                       gAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                          >  1:47930/1‑58 (MQ=255)
                                                         aCTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCAc           >  1:3275974/1‑71 (MQ=255)
                                                          cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGc                             >  1:3281777/1‑52 (MQ=255)
                                                           tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                <  1:2566668/48‑1 (MQ=255)
                                                           tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                <  1:2477942/48‑1 (MQ=255)
                                                           tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                                <  1:2191825/48‑1 (MQ=255)
                                                            aTCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTAt        >  1:1231917/1‑71 (MQ=255)
                                                             tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTa         <  1:2326530/69‑1 (MQ=255)
                                                             tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg       >  1:407809/1‑71 (MQ=255)
                                                             tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg       >  1:2480508/1‑71 (MQ=255)
                                                             tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGAGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg       <  1:472933/71‑1 (MQ=255)
                                                                  gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt  <  1:1916905/71‑1 (MQ=255)
                                                                  gTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt  <  1:4566378/71‑1 (MQ=255)
                                                                   tCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGt  >  1:806225/1‑70 (MQ=255)
                                                                    cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                 >  1:3725610/1‑38 (MQ=255)
                                                                    cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                 >  1:4474075/1‑38 (MQ=255)
                                                                    cTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCa                                 >  1:1198842/1‑38 (MQ=255)
                                                                    |                                                                    
GTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGT  >  minE/833146‑833282

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: