Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,740 G→A 11.6% L71L (CTG→CTA) ‡ 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,7400GA11.6% 113.4 / 6.2 43L71L (CTG→CTA) ‡narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (2/3);  total (21/22)
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.46e-01

CCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  minE/837674‑837807
                                                                  |                                                                   
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:3873386/1‑70 (MQ=255)
                                      tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                            >  1:3052113/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:270269/55‑1 (MQ=255)
                                                              tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                   <  1:275004/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)
                                                                 tGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa                 >  1:356391/1‑54 (MQ=255)
                                                                 tGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa                 >  1:318422/1‑54 (MQ=255)
                                                                  |                                                                   
CCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTG  >  minE/837674‑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: