Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,277,651 T→C 17.4% Q366Q (CAA→CAG 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,6510TC17.4% 114.4 / 14.0 74Q366Q (CAA→CAGbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (36/25);  new base C (7/6);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

ACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTG  >  minE/2277588‑2277712
                                                               |                                                             
acgacCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAg                                                          >  1:2220327/1‑69 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:1898742/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:2652026/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:1140641/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:442789/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:358266/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTCGCTGCAg                                                          >  1:387305/1‑68 (MQ=255)
 cgacCAAACTGAGAATTATCGTCGCAGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                                          >  1:861381/1‑68 (MQ=255)
        aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  <  1:323795/69‑1 (MQ=255)
        aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  <  1:416026/69‑1 (MQ=255)
        aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  <  1:1543739/69‑1 (MQ=255)
          tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                                                <  1:1835946/69‑1 (MQ=255)
          tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                                                <  1:3300946/69‑1 (MQ=255)
          tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                                                <  1:1669890/69‑1 (MQ=255)
              aGTTATCGTCACTGCTTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata                                              <  1:2651043/67‑1 (MQ=255)
              aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                                            >  1:2376560/1‑69 (MQ=255)
              aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                                            >  1:1108553/1‑69 (MQ=255)
              aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata                                              <  1:2777580/67‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:2314032/68‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:1743977/68‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:1826605/68‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:448117/68‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:1990297/68‑1 (MQ=255)
               attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                                           <  1:1994254/68‑1 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTTAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:1224645/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:2820658/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:2971345/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:331614/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:3329377/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:1497401/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                                  >  1:1307019/1‑59 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATGGATGcc                                        >  1:1805626/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:2566580/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:2783055/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:385943/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:2764923/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:447381/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:2651802/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:963159/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:1087534/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:1287081/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:2409963/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:1364270/1‑69 (MQ=255)
                  aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:1765728/1‑69 (MQ=255)
                  aTCGTCACTGCGTACTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                                        >  1:28792/1‑69 (MQ=255)
                     gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCgg                                      <  1:1088309/68‑1 (MQ=255)
                     gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                                     >  1:1008764/1‑69 (MQ=255)
                     gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                                     >  1:59919/1‑69 (MQ=255)
                     gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                                     >  1:2892193/1‑69 (MQ=255)
                     gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                                     >  1:2039455/1‑69 (MQ=255)
                            cGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACg                               >  1:671523/1‑68 (MQ=255)
                            cGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACg                               >  1:501596/1‑68 (MQ=255)
                              tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCTGTCagataga                                            >  1:1042312/1‑53 (MQ=255)
                              tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                                            >  1:2663400/1‑53 (MQ=255)
                              tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagatag                                             >  1:2695397/1‑52 (MQ=255)
                              tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                                                >  1:1720942/1‑49 (MQ=255)
                                   cgcCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTAt                                   <  1:2522180/57‑1 (MQ=255)
                                      cTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc                    <  1:2550816/69‑1 (MQ=255)
                                      cTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc                    <  1:598140/69‑1 (MQ=255)
                                       ttCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc                    <  1:2905848/68‑1 (MQ=255)
                                        tCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc                    <  1:2676456/67‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:433451/69‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:54765/69‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:1226990/69‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:2253620/69‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:223648/69‑1 (MQ=255)
                                               tAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:306603/69‑1 (MQ=255)
                                                aaaaTCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           >  1:2421316/1‑68 (MQ=255)
                                                aaaaTCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           <  1:491933/68‑1 (MQ=255)
                                                aaaaTCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCaa           >  1:1824702/1‑68 (MQ=255)
                                                      gCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAAcc    <  1:2498274/69‑1 (MQ=255)
                                                      gCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAAcc    <  1:3020572/69‑1 (MQ=255)
                                                       ccGACTTGTTGCAGCCCGGTCAGATAGATGCCCGTATTACGTTGATCATATCCATCCTCAACATAAcc    <  1:682958/68‑1 (MQ=255)
                                                        cGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  <  1:3435482/69‑1 (MQ=255)
                                                              gttgCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  >  1:1315447/1‑63 (MQ=255)
                                                              gttgCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  >  1:2170319/1‑63 (MQ=255)
                                                              gttgCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  >  1:1163961/1‑63 (MQ=255)
                                                              gttgCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  >  1:1103125/1‑63 (MQ=255)
                                                              gttgCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTg  >  1:3247982/1‑63 (MQ=255)
                                                               |                                                             
ACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATCCATCCTCAACATAACCTG  >  minE/2277588‑2277712

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: