Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,277,216 G→A 26.8% T511T (ACC→ACT 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,2160GA26.8% 147.9 / 49.4 126T511T (ACC→ACTbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (42/50);  new base A (27/7);  total (69/57)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

Only 100 of 126 total aligned reads displayed.

ATACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCGG  >  minE/2277153‑2277274
                                                               |                                                          
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:910881/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:2014359/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:252356/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:1534245/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:1533379/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:1518527/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:3407860/1‑69 (MQ=255)
aTACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAAt                                                       >  1:894200/1‑69 (MQ=255)
  aCAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCg                                                     <  1:3405583/69‑1 (MQ=255)
  aCAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCg                                                     <  1:188205/69‑1 (MQ=255)
  aCAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCg                                                     <  1:2751965/69‑1 (MQ=255)
  aCAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCg                                                     <  1:2562562/69‑1 (MQ=255)
  aCAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCg                                                     <  1:2125548/69‑1 (MQ=255)
   cAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGc                                                    >  1:1820381/1‑69 (MQ=255)
   cAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGc                                                    >  1:2870250/1‑69 (MQ=255)
       tGCCAGTCGTGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:499796/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTGACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:3056550/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:697624/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:1886991/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:2104406/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:2109679/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:2206486/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:292776/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:1396328/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:1352377/69‑1 (MQ=255)
       tGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTg                                                <  1:3109671/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:830028/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:2023619/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:3249442/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:939625/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:1493436/69‑1 (MQ=255)
           aGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcg                                            <  1:1139697/69‑1 (MQ=255)
            gTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcgc                                           <  1:1350429/69‑1 (MQ=255)
            gTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTgcgcgc                                           <  1:214806/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCa                                          <  1:1522972/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCa                                          <  1:101444/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCa                                          <  1:315992/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCa                                          <  1:329370/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCa                                          <  1:2321313/69‑1 (MQ=255)
             tCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCa                                          <  1:2291887/69‑1 (MQ=255)
              cGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCa                                          >  1:2800587/1‑68 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCAt                                         >  1:1758437/1‑57 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             <  1:2052307/69‑1 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             <  1:2028874/69‑1 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:1296681/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:2344246/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:525660/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:2384256/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:3264115/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:404290/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:83647/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:2699858/1‑69 (MQ=255)
                          tCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCAt                             >  1:2588763/1‑69 (MQ=255)
                                      tAGGACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTa                  <  1:353530/68‑1 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTa                  <  1:289814/68‑1 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTa                  <  1:2898079/68‑1 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTAt                 >  1:2899249/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCAt                                                  >  1:1898663/1‑36 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATa                            >  1:3156100/1‑58 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTa                  >  1:1075474/1‑68 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1958834/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:3109617/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:3119990/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:41867/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:567161/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1178652/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:764050/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:845673/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1694076/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:2823901/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:2802330/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1668006/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:2335862/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:2400802/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1598713/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGTAATCGCATTGCGAGCATCAACATAATCATAACTCACAGTAt                 >  1:358307/1‑69 (MQ=255)
                                      tAGCACGGCGTAACAACGGCGTGTCAGCAATCGCATTGCGCGCATCAACATAATCATAACTCACAGTAt                 >  1:1486536/1‑69 (MQ=255)
                                        gCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATg                <  1:3312348/68‑1 (MQ=255)
                                        gCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATg                <  1:1088349/68‑1 (MQ=255)
                                        gCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATg                <  1:3276677/68‑1 (MQ=255)
                                        gCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATg                <  1:2369958/68‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:120738/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:1921685/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:923558/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:1699369/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:3102069/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:2764155/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:2382041/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:2713169/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:2532475/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:1218428/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:1252893/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAg                    >  1:3308872/1‑51 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:1276787/69‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:1500561/69‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:3214370/69‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:390335/69‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:1818969/69‑1 (MQ=255)
                                                     acGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:3319575/69‑1 (MQ=255)
                                                         cGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCgg  <  1:2014310/65‑1 (MQ=255)
                                                               |                                                          
ATACAACTGCCAGTCGAGCTGGTATTTCACCTGCTGTTTAGCACGGCGTAACAACGGCGTGTCGGTAATCGCATTGCGCGCATCGACATAATCATAACTCACAGTATGCGTCAGTGGTCCGG  >  minE/2277153‑2277274

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: