Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,400,915 A→C 51.7% A172A (GCT→GCG yigB ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,400,9150AC51.7% 10.3 / 99.9 95A172A (GCT→GCGyigBpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (28/18);  new base C (29/20);  total (57/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCACGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAGCGCG  >  minE/2400848‑2400968
                                                                   |                                                     
ccacGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCa                                                       >  1:303112/1‑68 (MQ=255)
ccacGTCAGTGGTGAGATCGTCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAg                                                      >  1:2002894/1‑68 (MQ=255)
     tCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCa                                                   >  1:1841065/1‑67 (MQ=255)
          ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:3054535/69‑1 (MQ=255)
          ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:1974259/69‑1 (MQ=255)
          ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:2043056/69‑1 (MQ=255)
          ggtgAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTa                                            <  1:2592991/69‑1 (MQ=255)
             gagaTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACat                                         <  1:1146226/69‑1 (MQ=255)
             gagaTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACat                                         <  1:401171/69‑1 (MQ=255)
               gaTCGTCCCCTACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1409210/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:691020/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1702160/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:406804/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1659416/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1797347/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1478213/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:191353/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1409519/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2982314/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:34510/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:3404758/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2598814/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1284605/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1067277/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2100533/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2147927/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:1142637/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:3271970/60‑4 (MQ=255)
               gaTCGTCCCCGACATGTAAGATCTCACCGATCGGCACGTTAAGTTTTTCCGCCGCCagta                                                <  1:2481475/60‑4 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1956605/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:3184165/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:959562/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:952083/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1370992/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACACATc                                      >  1:1379013/1‑69 (MQ=255)
                aTCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTATCCGCAGCCAAAAAGTACATATc                                      >  1:3316488/1‑69 (MQ=255)
                aTCGACCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATc                                      >  1:1274250/1‑69 (MQ=255)
                  cGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCg                                     >  1:1350946/1‑68 (MQ=255)
                        cAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACg                              >  1:512223/1‑69 (MQ=255)
                        cAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACg                              >  1:2486762/1‑69 (MQ=255)
                              gTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:2484042/69‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1460458/68‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1883399/68‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:631514/68‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:2056456/68‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:1061653/68‑1 (MQ=255)
                               tAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTg                        <  1:854722/68‑1 (MQ=255)
                                 aGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgc                     >  1:3022510/1‑69 (MQ=255)
                                 aGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgc                     >  1:1044282/1‑69 (MQ=255)
                                   aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTgcgcgc                   >  1:3254618/1‑69 (MQ=255)
                                   aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:2166222/1‑69 (MQ=255)
                                   aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:800342/1‑69 (MQ=255)
                                   aTCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAgcgc                   >  1:1122722/1‑69 (MQ=255)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3236482/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3272100/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:440926/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:734785/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3301241/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3229626/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:711580/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:3222384/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:707668/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:683000/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:571602/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1037626/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1989534/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1344377/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1210763/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1523768/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1537530/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1081351/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1692197/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1782095/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:186064/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1997158/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2102633/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2360175/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2375132/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2473070/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1015/3‑50 (MQ=25)
                                       caccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:2992143/3‑50 (MQ=25)
                                       aaccGATCGGCACGTTAAGTTTTTCCGCCGCCAGTAAGTACATGTCGCTg                                  >  1:1075678/3‑50 (MQ=25)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCTCTGAACGGTATTGAg                      >  1:3386127/1‑55 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:96022/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:938017/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1040414/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1633524/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:499362/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:1989145/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCgg           >  1:2684652/1‑66 (MQ=255)
                                              cGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCACGTGCgg           >  1:647950/1‑66 (MQ=255)
                                                    gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1371458/69‑1 (MQ=255)
                                                    gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1443310/69‑1 (MQ=255)
                                                    gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:1227816/69‑1 (MQ=255)
                                                    gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:2213096/69‑1 (MQ=255)
                                                    gTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAgcgcg  <  1:3071822/69‑1 (MQ=255)
                                                                   |                                                     
CCACGTCAGTGGTGAGATCGTCCCCAACATGTAAGATCTCGCCGATCGGCACGTTGAGTTTTTCCGCAGCCAAAAAGTACATATCGCTGAACGGTTTTGAGCGCCCGTGCGGGCCAGCGCG  >  minE/2400848‑2400968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: