Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,198,249 C→T 41.6% R256R (CGC→CGT degS → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,198,2490CT41.6% 27.8 / 49.8 72R256R (CGC→CGTdegSserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (35/7);  new base T (20/10);  total (55/17)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

ATTCCTTTCCAGTTAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCAA  >  minE/2198187‑2198315
                                                              |                                                                  
aTTCCTTTCCAGTTAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTa                                                               <  1:26263/68‑1 (MQ=255)
 ttCCTTTCCAGTTAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTa                                                               <  1:77644/67‑1 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:2737701/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:1622527/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:87457/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:2082366/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:2425672/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:2494395/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:2709895/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:1012405/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:757694/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:3061288/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:752967/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:37240/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:422234/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:438323/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:549232/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:1348589/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:1493182/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:1275114/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:11703/1‑68 (MQ=255)
            ttAGCGACCAAAATTATGGATAAACTGATCCGTGATGGCCGTGTGATCCGTGGCTAcattggcattgg                                                   >  1:113637/1‑68 (MQ=255)
             tAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcg                                                 >  1:2486039/1‑69 (MQ=255)
             tAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcg                                                 >  1:2184773/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcgg                                                >  1:2671102/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcgg                                                >  1:2507529/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcgg                                                >  1:3009184/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcgg                                                >  1:1800168/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcgg                                                >  1:1645328/1‑69 (MQ=255)
              aGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggcg                                                 >  1:2428023/1‑68 (MQ=255)
                 aaCCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGAc                                              >  1:3083623/1‑68 (MQ=255)
                    cAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGAc                                              >  1:2768212/1‑65 (MQ=255)
                    cAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGAc                                              >  1:736118/1‑65 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcg                                                    >  1:1329095/1‑55 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:3231873/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:3237292/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:1126991/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:2551247/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:91353/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:466885/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:783800/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:204287/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:1581708/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:1432305/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTgag                                         >  1:2285079/1‑66 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTga                                          >  1:103645/1‑65 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTga                                          >  1:894326/1‑65 (MQ=255)
                        aTTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTga                                          >  1:1473696/1‑65 (MQ=255)
                        aTTATGGATAAGCTGATACGCGAAGGTCGCGTGATCCGCGGCTACATTgg                                                         >  1:21370/1‑50 (MQ=255)
                             ggATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcac                                  <  1:819860/68‑1 (MQ=255)
                             ggATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcac                                  <  1:3127900/68‑1 (MQ=255)
                             ggATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcac                                  <  1:2349419/68‑1 (MQ=255)
                             ggATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcac                                  <  1:1102579/68‑1 (MQ=255)
                             ggATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcac                                  <  1:962678/68‑1 (MQ=255)
                              gATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcacca                                >  1:202966/1‑69 (MQ=255)
                              gATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGcacca                                >  1:3340877/1‑69 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATcggc                                                  >  1:680751/1‑48 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTg                             >  1:3336146/1‑69 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTg                             >  1:2819922/1‑69 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTg                             >  1:2706315/1‑69 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTg                             >  1:1432970/1‑69 (MQ=255)
                                 aaGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTg                             >  1:1481799/1‑69 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:1032583/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:641597/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:3213427/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:3045935/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:2715422/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:2325445/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:1934950/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:1532048/57‑1 (MQ=255)
                                     tGATCCGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:1318796/57‑1 (MQ=255)
                                     tGATACGTGATGGCCGTGTGATCCGTGGCTACATTGGCATTGGTGGACGTGAGATCg                                     <  1:982044/57‑1 (MQ=255)
                                                            cgcgGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCaa  <  1:3426088/69‑1 (MQ=255)
                                                            cgcgGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCaa  <  1:758788/69‑1 (MQ=255)
                                                            cgcgGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCaa  <  1:2489312/69‑1 (MQ=255)
                                                            cgcgGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCaa  <  1:1682039/69‑1 (MQ=255)
                                                              |                                                                  
ATTCCTTTCCAGTTAGCAACCAAAATTATGGATAAGCTGATCCGCGATGGTCGCGTGATCCGCGGCTACATTGGTATCGGCGGACGTGAGATCGCACCACTGCACGCGCAGGGCGGTGGTATAGATCAA  >  minE/2198187‑2198315

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: