Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,357,714 C→T 34.0% V470I (GTC→ATC)  yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,357,7140CT34.0% 104.4 / 71.5 129V470I (GTC→ATC) yehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (54/31);  new base T (37/7);  total (91/38)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 129 total aligned reads displayed.

TTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTA  >  minE/1357653‑1357780
                                                             |                                                                  
ttGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtt                                                              >  1:2670122/1‑68 (MQ=255)
ttGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtt                                                              >  1:175884/1‑68 (MQ=255)
ttGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtt                                                              >  1:1630075/1‑68 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGTCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            <  1:140828/68‑1 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            <  1:1549830/68‑1 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            <  1:957230/68‑1 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            >  1:1051144/1‑68 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            <  1:2555571/68‑1 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGtttt                                                            <  1:2726034/68‑1 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTtc                                                           >  1:1503321/1‑69 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTtc                                                           >  1:1581392/1‑69 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTtc                                                           >  1:2099141/1‑69 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTtc                                                           >  1:3368981/1‑69 (MQ=255)
  ggCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTtc                                                           >  1:831744/1‑69 (MQ=255)
        aaaGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACg                                                     >  1:114496/1‑69 (MQ=255)
        aaaGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACg                                                     >  1:2304447/1‑69 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:49512/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:2924930/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:2909996/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:2733615/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:2440774/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:2554164/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:814875/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:1053926/69‑1 (MQ=255)
                           aaaTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTg                                  <  1:832059/69‑1 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGc                                                   >  1:3008073/1‑50 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:53392/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:2475916/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:1078257/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:1208295/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:3416309/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:2758780/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:2682210/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:829353/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:101689/1‑66 (MQ=255)
                             aTCGTTCTGGTAATAATCGGCAATGTCATAAACATCGTTTTCATCACGGCGATAATTATCTTTAtt                                   >  1:1107698/1‑66 (MQ=255)
                                         ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                    <  1:1993475/69‑1 (MQ=255)
                                         ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                    <  1:3149362/69‑1 (MQ=255)
                                         ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                    <  1:1494450/69‑1 (MQ=255)
                                         ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                    <  1:3333752/69‑1 (MQ=255)
                                           aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                   <  1:888173/68‑1 (MQ=255)
                                           aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                   <  1:2260476/68‑1 (MQ=255)
                                           aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                   <  1:2514337/68‑1 (MQ=255)
                                           aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                   <  1:257050/68‑1 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAAt                                             >  1:869868/1‑41 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACg                      >  1:3246537/1‑64 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:355597/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:3199891/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:3324299/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:3435005/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2521419/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:471511/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:691348/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:761320/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:965058/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2538497/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1009725/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1116903/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1211199/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1225090/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1294988/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:131130/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1441863/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:144747/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1939048/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2036950/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2217821/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2272103/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2294010/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2284952/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:2482539/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                 >  1:1005576/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATc                 >  1:2407509/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATc                 >  1:1661378/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATc                 >  1:1619331/1‑69 (MQ=255)
                                            atCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATc                 >  1:1569833/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc       >  1:2064948/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc       >  1:3258781/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc       >  1:2763998/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc       >  1:3339207/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc       >  1:77995/1‑69 (MQ=255)
                                                      tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCAAAACGTg                    >  1:3241623/1‑56 (MQ=255)
                                                        aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa             <  1:3175851/61‑1 (MQ=255)
                                                        aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa             <  1:1664257/61‑1 (MQ=255)
                                                        aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa             <  1:2808984/61‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGt   >  1:3039693/1‑68 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGt   >  1:3193636/1‑68 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:1630985/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  <  1:761241/69‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:1405258/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  <  1:1704472/69‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:1870769/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:238934/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:2910083/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:1116968/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  <  1:870463/69‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  <  1:2943730/69‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  <  1:926338/69‑1 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:1018319/1‑69 (MQ=255)
                                                           gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa  >  1:258333/1‑69 (MQ=255)
                                                             |                                                                  
TTGGCGGAAAAGCTATTTTTGCGGCCAAAATCGTTCTGGTAATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTA  >  minE/1357653‑1357780

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: