Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,906,404 T→A 39.3% A387A (GCA→GCT treB ← fused trehalose(maltose)‑specific enzyme IIBC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,906,4040TA39.3% 72.8 / 111.2 117A387A (GCA→GCTtreBfused trehalose(maltose)‑specific enzyme IIBC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base A (23/23);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 117 total aligned reads displayed.

ACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTG  >  minE/2906344‑2906469
                                                            |                                                                 
accagaACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTc                                                           <  1:471320/69‑1 (MQ=255)
accagaACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTc                                                           <  1:3187436/69‑1 (MQ=255)
accagaACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTc                                                           <  1:1232862/69‑1 (MQ=255)
 ccagaACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTc                                                           <  1:2973255/68‑1 (MQ=255)
             atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                >  1:3228806/1‑67 (MQ=255)
             atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                >  1:1838635/1‑67 (MQ=255)
             atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                >  1:1956625/1‑67 (MQ=255)
             atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                >  1:666157/1‑67 (MQ=255)
             atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                >  1:1022409/1‑67 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGcc                                             >  1:1303382/1‑68 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:674280/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:2258712/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:2317671/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:328803/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:1568588/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:457436/1‑69 (MQ=255)
               cgcgCACAGCATCGGGAAGAGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                            >  1:1019005/1‑69 (MQ=255)
                   cacaGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCg                                        <  1:1531533/69‑1 (MQ=255)
                   cacaGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCg                                        <  1:2362312/69‑1 (MQ=255)
                    acaGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCg                                        <  1:8956/68‑1 (MQ=255)
                     caGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      <  1:3215992/69‑1 (MQ=255)
                     caGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      <  1:1287043/69‑1 (MQ=255)
                      aGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      <  1:874449/68‑1 (MQ=255)
                       gCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      >  1:2778165/1‑67 (MQ=255)
                       gCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      >  1:3297178/1‑67 (MQ=255)
                        cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGgaga                                    >  1:603789/1‑68 (MQ=255)
                        cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                   >  1:2621181/1‑69 (MQ=255)
                        cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                   >  1:1772158/1‑69 (MQ=255)
                        cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                   >  1:1437456/1‑69 (MQ=255)
                        cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                   <  1:10319/69‑1 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:402579/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:2466910/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:2320863/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:2911731/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:2266500/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:1442231/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:3143102/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:849153/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:404390/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:555006/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:1762985/1‑69 (MQ=255)
                         aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                  >  1:1717678/1‑69 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGc                                                                 >  1:546982/1‑35 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:3005774/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:534425/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:3264554/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:3277111/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:2845550/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:2794604/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:3418075/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:2657429/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:1718341/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:1372707/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:107721/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              >  1:2108906/1‑54 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:2357954/54‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:2184681/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:1426276/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:3108588/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:2209225/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:966132/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:2726469/69‑1 (MQ=255)
                            gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                               <  1:2569471/69‑1 (MQ=255)
                                    atatTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                              <  1:2044144/46‑1 (MQ=255)
                                     tatTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      <  1:1245686/53‑1 (MQ=255)
                                     tatTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                      <  1:1832175/53‑1 (MQ=255)
                                         tCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGgag                   <  1:2936930/68‑1 (MQ=255)
                                         tCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGgag                   <  1:2854057/68‑1 (MQ=255)
                                         tCAGGTTGATGCCCTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGgag                   <  1:3145955/68‑1 (MQ=255)
                                             gTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGAtct               <  1:1057133/68‑1 (MQ=255)
                                             gTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGAtct               <  1:1053174/68‑1 (MQ=255)
                                              ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg             >  1:2726201/1‑69 (MQ=255)
                                              ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg             >  1:1565207/1‑69 (MQ=255)
                                              ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg             >  1:1750821/1‑69 (MQ=255)
                                              ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg             >  1:1452155/1‑69 (MQ=255)
                                              ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg             >  1:2237150/1‑69 (MQ=255)
                                               tGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcgc            >  1:2981947/1‑69 (MQ=255)
                                               tGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcgc            >  1:2268195/1‑69 (MQ=255)
                                               tGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGAGGCAGGCACGGAGATCTcgc            >  1:1455829/1‑69 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:1449906/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:867190/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:414113/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:1615740/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:1786529/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3431678/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3423736/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3421443/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3371411/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3364368/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3347156/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:2686109/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:1876738/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:1990295/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:2071896/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:3068068/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:2963934/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:2501151/69‑1 (MQ=255)
                                                       tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg    <  1:2704103/69‑1 (MQ=255)
                                                         cATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTTATtgtg  <  1:440542/69‑1 (MQ=255)
                                                         cATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtgtg  <  1:1378048/69‑1 (MQ=255)
                                                            |                                                                 
ACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTG  >  minE/2906344‑2906469

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: