Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,222,892 T→C 43.5% G489G (GGT→GGC yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,222,8920TC43.5% 15.9 / 66.7 85G489G (GGT→GGCyebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/36);  new base C (14/23);  total (26/59)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

TTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTT  >  minE/1222827‑1222952
                                                                 |                                                            
tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTgag                                                           >  1:2124943/1‑69 (MQ=255)
tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTgag                                                           >  1:1692588/1‑69 (MQ=255)
tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTga                                                            >  1:2008975/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:3327919/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:2619194/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:2393636/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:2912283/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:2927065/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:2988393/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:1905203/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:157121/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:144889/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:1441117/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:732136/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:1075892/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:1170015/1‑68 (MQ=255)
       gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                                                     >  1:804513/1‑68 (MQ=255)
            aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc                                               >  1:2502480/1‑69 (MQ=255)
            aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc                                               >  1:388387/1‑69 (MQ=255)
            aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc                                               >  1:1289915/1‑69 (MQ=255)
                   cataaTGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1060533/57‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGTCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1990151/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1307341/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1442937/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:2418095/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1553681/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:2971013/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:2014779/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:1340909/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:2992223/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:3182979/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:3319001/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAACGCGAGATGAAAGCCa                                                 <  1:122040/60‑1 (MQ=255)
                   cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCCATAAAGGCGAGATGAAAGCCa                                                 <  1:2620653/60‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAGGCCAGCTATCCACTg                                      <  1:1398691/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:3390389/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:308074/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2954673/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:844736/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:992538/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2821099/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2621188/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:556369/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1075718/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:587139/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1897147/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1574888/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1685291/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1547393/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:1396870/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2053302/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2100379/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2129709/68‑1 (MQ=255)
                      aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg                                      <  1:2257479/68‑1 (MQ=255)
                         ggATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTAt                                   >  1:1759572/1‑68 (MQ=255)
                             cAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                                                 <  1:2158091/50‑1 (MQ=39)
                                tGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCAAGCTATCCACTGTAt                                   >  1:1865580/1‑61 (MQ=255)
                                            tATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcgc                >  1:1275555/1‑68 (MQ=255)
                                            tATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAgcgc                >  1:1835017/1‑68 (MQ=255)
                                                cTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgag           <  1:3169455/69‑1 (MQ=255)
                                                cTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgag           >  1:1341745/1‑69 (MQ=255)
                                                cTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgag           >  1:822310/1‑69 (MQ=255)
                                                 tGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGgag           <  1:1154799/68‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:744312/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:746964/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:997172/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:964728/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:494278/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:461018/67‑1 (MQ=255)
                                                     gggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAAc        <  1:2693277/67‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:3414160/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:337841/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:2895222/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:2762083/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:1164580/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:912731/62‑1 (MQ=255)
                                                     gggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:2321169/62‑1 (MQ=255)
                                                      ggCGATATAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGc     <  1:3329795/69‑1 (MQ=255)
                                                      ggCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGc     <  1:414700/69‑1 (MQ=255)
                                                      ggCGATAAAGGCGAGATGAAAGCCAGCTATCCACTGTATGCCAATCAGGAAAAAGCGCTgg             <  1:807216/61‑1 (MQ=255)
                                                         gATAAAGGTGATATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACATCCtt  <  1:2311612/69‑1 (MQ=255)
                                                         gATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCtt  <  1:3412889/69‑1 (MQ=255)
                                                         gATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCtt  <  1:1691657/69‑1 (MQ=255)
                                                         gATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCtt  <  1:2238134/69‑1 (MQ=255)
                                                         gATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCtt  <  1:2866755/69‑1 (MQ=255)
                                                                 |                                                            
TTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTATGCCAATTTGGAAAAAGCGCTGGAGAACAGCCTT  >  minE/1222827‑1222952

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: