Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,175 C→G 10.5% L545V (CTG→GTG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,1750CG10.5% 159.3 / 8.2 57L545V (CTG→GTG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (3/3);  total (30/27)
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.96e-01

GCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCG  >  minE/833109‑833238
                                                                  |                                                               
gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGaa                                                             >  1:4025038/1‑71 (MQ=255)
gCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGaa                                                             >  1:4020747/1‑71 (MQ=255)
       aaTTTGCCGATACCGCCCATAAAACGCATGGGCGCTCGATGATTATCCTCGGCGCAGGTGTTAACCACTgg                                                      >  1:2400404/1‑71 (MQ=255)
       aaTTTGCCGATAACGCTGATAAAACGCACTGTCGTTCGATGATTATCGTCGGTGCTGGGCTGAACCACTgg                                                      <  1:2529984/71‑1 (MQ=255)
         tttGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGt                                                     >  1:533650/1‑70 (MQ=255)
         tttGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTa                                                    <  1:4128758/71‑1 (MQ=255)
             ccGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCAc                                                >  1:3109312/1‑71 (MQ=255)
                aTAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc                                             <  1:325252/71‑1 (MQ=255)
                      cTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATg                                       <  1:3213866/71‑1 (MQ=255)
                          tAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAAc                                    >  1:30764/1‑70 (MQ=255)
                           aaaaCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTa                                  >  1:4029251/1‑71 (MQ=255)
                            aaaCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGaa                                     >  1:418991/1‑67 (MQ=255)
                            aaaCGCACGGTCGTTCGATGATTATCGTCGGTGCCGGGCTGAACCACTGGTATCACCTCGata                                         >  1:711676/1‑63 (MQ=255)
                             aaCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATc                                <  1:1195329/71‑1 (MQ=255)
                              aCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCa                                                          >  1:4017881/1‑44 (MQ=255)
                              aCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACt                                                        >  1:2229296/1‑46 (MQ=255)
                              aCGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCg                               >  1:2684304/1‑71 (MQ=255)
                                gCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTc                                             <  1:2365217/55‑1 (MQ=255)
                                 cACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTgg                            <  1:1031779/71‑1 (MQ=255)
                                   cGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                          >  1:3584905/1‑71 (MQ=255)
                                    ggTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGa                                           <  1:2747546/53‑1 (MQ=255)
                                     gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGa                                           <  1:3022604/52‑1 (MQ=255)
                                     gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                        <  1:1479669/71‑1 (MQ=255)
                                       cGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGAt                      >  1:3966255/1‑71 (MQ=255)
                                        gctCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATc                     <  1:1547010/69‑1 (MQ=255)
                                        gTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGTTATCACCTCGa                                           <  1:2060640/49‑1 (MQ=255)
                                         ctCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCa                    <  1:3507793/70‑1 (MQ=255)
                                           cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                   >  1:1957518/1‑70 (MQ=255)
                                           cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                   >  1:743291/1‑70 (MQ=255)
                                           cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                   <  1:838513/70‑1 (MQ=255)
                                           cGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTa                                  >  1:1993408/1‑55 (MQ=255)
                                           agatgatTATCGTCGGTGCGGGGCTGAACCACTGGTa                                                    >  1:2614957/2‑37 (MQ=255)
                                            gatgatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCaa                   <  1:3581984/69‑1 (MQ=255)
                                               gatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGc              >  1:957155/1‑71 (MQ=255)
                                                 ttATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGaa                                     >  1:1225912/1‑46 (MQ=255)
                                                 ttATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGa                       >  1:4110667/1‑60 (MQ=255)
                                                  tATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATc                                >  1:3151537/1‑50 (MQ=255)
                                                   aTCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                          >  1:3046272/1‑55 (MQ=255)
                                                   aTCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAt          >  1:2868683/1‑71 (MQ=255)
                                                      gtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTgg                            <  1:2633458/50‑1 (MQ=255)
                                                      gtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                           <  1:3842512/51‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                             <  1:1120413/48‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                             <  1:1881647/48‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                             <  1:2505111/48‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                           <  1:3373450/50‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                           <  1:382336/50‑1 (MQ=255)
                                                       tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc      >  1:3716220/1‑71 (MQ=255)
                                                        cgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                          >  1:598088/1‑50 (MQ=255)
                                                        cgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCt                         >  1:2940499/1‑51 (MQ=255)
                                                         ggTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                           >  1:2095725/1‑48 (MQ=255)
                                                         ggTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCt                         >  1:103783/1‑50 (MQ=255)
                                                            gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                   <  1:2675433/53‑1 (MQ=255)
                                                            gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg  >  1:3018486/1‑70 (MQ=255)
                                                             cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                        <  1:4102981/47‑1 (MQ=255)
                                                             cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                        <  1:1163247/47‑1 (MQ=255)
                                                                 gCTGAACCACTGGTATCACCTCGATATGAACTATCGTggg                           <  1:3460443/40‑2 (MQ=255)
                                                                 gCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                           <  1:539666/40‑1 (MQ=255)
                                                                  |                                                               
GCCCGTGAATTTGCCGATAACGCTGATAAAACGCACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCG  >  minE/833109‑833238

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: