Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,277,555 T→C 26.2% E398E (GAA→GAG btuB ← vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,277,5550TC26.2% 60.8 / 19.9 57E398E (GAA→GAGbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (14/28);  new base C (4/11);  total (18/39)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTC  >  minE/2277491‑2277611
                                                                |                                                        
gCCCCAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTccc                                                      <  1:1478500/69‑1 (MQ=255)
  cccAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCa                                                     <  1:3296656/68‑1 (MQ=255)
  cccAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCa                                                     <  1:2651555/68‑1 (MQ=255)
  cccAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCa                                                     <  1:2339538/68‑1 (MQ=255)
  cccAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCa                                                     <  1:858752/68‑1 (MQ=255)
    aaGATTTGGTGCCTTATAAGATGTTCCGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAAcc                                                  >  1:298350/2‑69 (MQ=255)
    aaGATTTGGTGCCTTATAAGATGTTCCGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAAcc                                                  >  1:1842369/2‑69 (MQ=255)
    aaGATTTGGTGCCTTATAAGATGTTCCGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAAcc                                                  >  1:2737387/2‑69 (MQ=255)
    aaGATTTGGTGCCTTATAAGATGTTCCGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAAc                                                   >  1:370574/2‑68 (MQ=255)
         ttGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGgcgc                                             >  1:2459959/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:1599941/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:2550440/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:60552/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:558927/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:2251240/1‑69 (MQ=255)
                        aTGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              >  1:2179826/1‑69 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:2551495/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:902010/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:2958202/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:2607605/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:307051/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:2993824/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:1829107/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:1797997/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:637438/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCGGCGCTGGTTTGc                                     <  1:1599714/57‑1 (MQ=255)
                          gttccGTAGGAAGCAATAAAGCGATAACCTTCGATGAACTCCCAACCAGCGCTGGTTTGc                                     <  1:15644/57‑1 (MQ=255)
                           tCCCGTAGGAAGCAATGAAGCGATAACCTTCGATTAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCAt                           <  1:657398/69‑1 (MQ=255)
                                     aGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              <  1:1337766/56‑1 (MQ=255)
                                     aGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              <  1:218826/56‑1 (MQ=255)
                                     aGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              <  1:502965/56‑1 (MQ=255)
                                     aGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTc                              <  1:1135569/56‑1 (MQ=255)
                                              gCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt        <  1:1156755/69‑1 (MQ=255)
                                              gCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt        <  1:2910390/69‑1 (MQ=255)
                                                 aTAACCTTCGATGAATTCCCAACCGGCTCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc     >  1:1156907/1‑69 (MQ=255)
                                                 aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt        <  1:1718846/66‑1 (MQ=255)
                                                 aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc     <  1:2015337/69‑1 (MQ=255)
                                                 aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc     >  1:2518156/1‑69 (MQ=255)
                                                   taCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:175132/68‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:176609/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:96220/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   >  1:1208041/1‑69 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:747021/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:163970/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:1692704/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   >  1:636213/1‑69 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:515162/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:49696/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:376700/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:1822995/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:3136614/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:189982/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:2272650/69‑1 (MQ=255)
                                                   aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt   <  1:2322560/69‑1 (MQ=255)
                                                    aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc  >  1:1841791/1‑69 (MQ=255)
                                                    aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc  >  1:807695/1‑69 (MQ=255)
                                                    aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc  >  1:2017643/1‑69 (MQ=255)
                                                                |                                                        
GCCCCAGATTTGGTGCCTTATAAGATGTCCCGTAGGAAGCAATGAAGCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTC  >  minE/2277491‑2277611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: