Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,357,757 A→T 45.7% V455V (GTT→GTA yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,357,7570AT45.7% 21.1 / 106.0 94V455V (GTT→GTAyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/21);  new base T (36/7);  total (66/28)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGA  >  minE/1357693‑1357825
                                                                |                                                                    
aataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCaaa                                                                     <  1:3299561/66‑1 (MQ=255)
 ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                                                                 <  1:1993475/69‑1 (MQ=255)
 ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                                                                 <  1:1494450/69‑1 (MQ=255)
 ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                                                                 <  1:3149362/69‑1 (MQ=255)
 ataatCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTg                                                                 <  1:3333752/69‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:2546507/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:257050/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:1589436/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:2260476/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:2514337/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:888173/68‑1 (MQ=255)
   aatCGGCAATGTCATAGATATCGTTTTCATAACGGCGATAATTATCTTTATTGTTTGCCCATACGTGa                                                                <  1:229277/68‑1 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACg                                                                   >  1:3246537/1‑64 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2482539/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2620685/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2521419/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1005576/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2036950/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2294010/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2217821/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2284952/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2272103/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2538497/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2499493/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2909104/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:2925923/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:3199891/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:3235754/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:3324299/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:3435005/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:355597/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:860386/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:965058/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1939048/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1441863/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1418098/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:131130/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:144747/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1294988/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1499353/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1225090/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1211199/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1116903/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:104223/1‑69 (MQ=255)
    atCGGCAATGTCATAGATATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCATACGTGATc                                                              >  1:1009725/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)
             gTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGt                                                     >  1:792135/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:3339207/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:1174116/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:355770/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:3258781/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:77995/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:2064948/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTc                                                    >  1:2763998/1‑69 (MQ=255)
              tCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCAAAACGTg                                                                 >  1:3241623/1‑56 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGTTCGTTa                                                          <  1:3381877/61‑1 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa                                                          <  1:3175851/61‑1 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa                                                          <  1:2808984/61‑1 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa                                                          <  1:1664257/61‑1 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa                                                          <  1:1632649/61‑1 (MQ=255)
                aTAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTa                                                          <  1:751005/61‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGt                                                >  1:3039693/1‑68 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGt                                                >  1:3193636/1‑68 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:2177174/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               <  1:926338/69‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1116968/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               <  1:1704472/69‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               <  1:761241/69‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1833056/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1870769/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:2910083/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               <  1:870463/69‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1018319/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               <  1:2943730/69‑1 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:238934/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1405258/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1630985/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:258333/1‑69 (MQ=255)
                   gACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTa                                               >  1:1262618/1‑69 (MQ=255)
                      aTCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAAt                                             >  1:368067/1‑68 (MQ=255)
                         gTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCAc                                          >  1:3309974/1‑68 (MQ=255)
                               atcaCGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGcgacg                                    <  1:2911681/68‑1 (MQ=255)
                               atcaCGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGcgacg                                    <  1:1237190/68‑1 (MQ=255)
                                          aaTTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAgg                        >  1:113075/1‑69 (MQ=255)
                                          aaTTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAgg                        >  1:1588354/1‑69 (MQ=255)
                                              aTCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGgccgc                    <  1:1749414/69‑1 (MQ=255)
                                              aTCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGgccgc                    <  1:497451/69‑1 (MQ=255)
                                              aTCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGgccgc                    <  1:1675106/69‑1 (MQ=255)
                                              aTCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGgccgc                    <  1:935240/69‑1 (MQ=255)
                                                                aaCGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCg   >  1:1071688/1‑68 (MQ=255)
                                                                aaCGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGa  >  1:3379587/1‑69 (MQ=255)
                                                                |                                                                    
AATAATCGGCAATGTCATAGACATCGTTTTCATCACGGCGATAATTATCTTTATTGTTTGCCCAAACGTGATCGTTAAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGA  >  minE/1357693‑1357825

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: