Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 99,964 G→A 50.3% I167I (ATC→ATT hofC ← assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE99,9640GA50.3% ‑4.1 / 142.3 147I167I (ATC→ATThofCassembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base A (60/14);  total (97/50)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-02

Only 100 of 159 total aligned reads displayed.

GCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCC  >  minE/99909‑100030
                                                       |                                                                  
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:92532/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 >  1:56620/1‑69 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:3107085/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:2876586/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1219437/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1776730/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1335206/69‑1 (MQ=255)
      cAAAATGCAGCATTGCCACAACCACCAAGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:462758/69‑1 (MQ=255)
       aaaaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1058807/68‑1 (MQ=255)
       aaaaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1499986/68‑1 (MQ=255)
       aaaaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1496545/68‑1 (MQ=255)
       aaaaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTg                                                 <  1:1354886/68‑1 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:508807/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:1403349/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:2350284/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:1320483/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:3305115/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:2772812/1‑68 (MQ=255)
         aaTGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGAt                                               >  1:1680499/1‑68 (MQ=255)
             cagcaTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCAc                                          <  1:814048/69‑1 (MQ=255)
             cagcaTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCAc                                          <  1:1470803/69‑1 (MQ=255)
             cagcaTTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCAc                                          <  1:2513336/69‑1 (MQ=255)
                    gcaacaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:1778304/4‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:1392883/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:1836982/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTc                                    >  1:2732516/1‑68 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:3252790/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:2573752/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:2809695/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:3011808/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:3029555/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:3033607/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:376836/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:334655/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:803061/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:113317/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:445240/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:538027/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:1558276/1‑69 (MQ=255)
                    gCCACAACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCg                                   >  1:1728982/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:3062919/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCTGt                                 >  1:1291584/1‑68 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGt                                 >  1:772654/1‑68 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGt                                 >  1:2774649/1‑68 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:1624452/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:280999/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:1147819/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:292434/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:3086989/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:1433025/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:3224683/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:617705/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:2656997/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:1130280/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:2360449/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:2310690/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:2175647/1‑69 (MQ=255)
                       acaACCACCATGATTGCCATCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTc                                >  1:428869/1‑69 (MQ=255)
                          accaccATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCaa                              >  1:1672672/1‑68 (MQ=255)
                          accaccATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCaa                              >  1:2447024/1‑68 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:3300753/69‑1 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:550606/69‑1 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:1154516/69‑1 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:1696270/69‑1 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:2650911/69‑1 (MQ=255)
                           ccaccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:2183332/69‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:1026088/68‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:757897/68‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:309571/68‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:1800781/68‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:551943/68‑1 (MQ=255)
                            caccaTGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACt                            <  1:546078/68‑1 (MQ=255)
                                         aTCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              >  1:306332/1‑69 (MQ=255)
                                         aTCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              >  1:1673104/1‑69 (MQ=255)
                                         aTCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              >  1:1878646/1‑69 (MQ=255)
                                         aTCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              >  1:226280/1‑69 (MQ=255)
                                         aTCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              >  1:3018853/1‑69 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:2488194/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:3384671/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:1448250/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:1323899/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:2080634/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:643294/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:595084/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:2580555/69‑1 (MQ=255)
                                         aTCGCTAAAATGATAATGGCATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCtttt              <  1:3191454/69‑1 (MQ=255)
                                            gCTAAAATGATGATGTGATAACGTAACGCTGATTTCa                                           >  1:2914500/1‑37 (MQ=255)
                                            gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTg            >  1:2438599/1‑68 (MQ=255)
                                            gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTg            >  1:587084/1‑68 (MQ=255)
                                            gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTg            >  1:2067109/1‑68 (MQ=255)
                                            gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTg            >  1:308038/1‑68 (MQ=255)
                                            gCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTg            >  1:877963/1‑68 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACgcgc   >  1:3145875/1‑68 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:3223341/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:1814211/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:1885471/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:465552/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:2126332/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:3256240/1‑69 (MQ=255)
                                                     atgatgGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGcc  >  1:3043410/1‑69 (MQ=255)
                                                       |                                                                  
GCAGAACAAAATGCAGCATTGCCACAACCACCATGATTGCCATCGCTAAAATGATGATGGGATAACGTAACGCTGATTTCACTTTGTCGGTCAACTGACGCTGGGCTTTTTGCTGACGCGCC  >  minE/99909‑100030

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: