Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,738 C→T 11.9% L71L (CTG→TTG) ‡ narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,7380CT11.9% 112.1 / 5.4 42L71L (CTG→TTG) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base T (2/3);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  minE/837668‑837807
                                                                      |                                                                     
gATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGc                                                                       <  1:2548172/71‑1 (MQ=255)
 aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      >  1:498578/1‑71 (MQ=255)
 aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                                      >  1:174865/1‑71 (MQ=255)
      ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc                                                                 >  1:2300660/1‑71 (MQ=255)
       cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc                                                                 <  1:3764346/70‑1 (MQ=255)
         ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                               >  1:1505350/1‑70 (MQ=255)
          gCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                               <  1:2446114/69‑1 (MQ=255)
             tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                          <  1:2020350/71‑1 (MQ=255)
               gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg                                                         <  1:999754/70‑1 (MQ=255)
               gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTGCGGCATGTTAACGCCGCACTg                                                          <  1:1437546/69‑1 (MQ=255)
               gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGa                                                        >  1:1273085/1‑71 (MQ=255)
                 aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATg                                                      >  1:1463209/1‑71 (MQ=255)
                  tcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                                      >  1:2560582/2‑70 (MQ=255)
                    tGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg                                                         <  1:1092000/65‑1 (MQ=255)
                      ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa                                                 <  1:1064298/71‑1 (MQ=255)
                       tCCATTTCGGGATTCTGGGGATTTTTGTCGGTCAGTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa                                                 <  1:2305298/70‑1 (MQ=255)
                       tCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATg                                                  >  1:3013147/1‑69 (MQ=255)
                            atCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                                    <  1:2855573/46‑1 (MQ=255)
                                      tGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGTa                                                    >  1:723555/1‑52 (MQ=255)
                                        gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                                        >  1:210283/1‑62 (MQ=255)
                                        gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                                        >  1:1730880/1‑62 (MQ=255)
                                        gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                               <  1:766058/71‑1 (MQ=255)
                                         ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc                                              >  1:941868/1‑55 (MQ=255)
                                           aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            <  1:1054019/71‑1 (MQ=255)
                                           aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            <  1:994830/71‑1 (MQ=255)
                                            tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            >  1:3052113/1‑70 (MQ=255)
                                            tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            >  1:3873386/1‑70 (MQ=255)
                                               ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                        >  1:2970484/1‑71 (MQ=255)
                                               ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaa                         <  1:2830759/70‑1 (MQ=255)
                                                tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                            >  1:3079998/1‑50 (MQ=255)
                                                tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                            >  1:3248328/1‑50 (MQ=255)
                                                  tCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                        >  1:2326049/1‑68 (MQ=255)
                                                      tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                                <  1:403390/40‑1 (MQ=255)
                                                      tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                                <  1:304793/40‑1 (MQ=255)
                                                        aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATg               >  1:608373/1‑71 (MQ=255)
                                                          ttcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATgg              >  1:3934213/1‑70 (MQ=255)
                                                            cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa           >  1:3931720/1‑71 (MQ=255)
                                                                ggTATGCTGACGCCGCACCGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                   <  1:257214/59‑1 (MQ=255)
                                                                    tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGa                          <  1:1704555/48‑1 (MQ=255)
                                                                    tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                   <  1:275004/55‑1 (MQ=255)
                                                                    tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                   <  1:270269/55‑1 (MQ=255)
                                                                    tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGc    <  1:640231/70‑1 (MQ=255)
                                                                    tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGc    <  1:1918380/70‑1 (MQ=255)
                                                                      cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg  <  1:2459600/70‑1 (MQ=255)
                                                                      |                                                                     
GATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  minE/837668‑837807

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: