Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,418,890 G→T 47.2% P375T (CCT→ACT)  wzyE ← predicted Wzy protein involved in ECA polysaccharide chain elongation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,418,8900GT47.2% 12.4 / 56.0 55P375T (CCT→ACT) wzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base G (16/13);  new base T (18/8);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACGATCGCCCCG  >  minE/2418822‑2418957
                                                                    |                                                                   
cAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                                                                     <  1:2432164/69‑1 (MQ=255)
cAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                                                                     <  1:2717515/69‑1 (MQ=255)
cAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                                                                     <  1:3324369/69‑1 (MQ=255)
 aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                                                                     <  1:953096/68‑1 (MQ=255)
 aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGc                                                                    >  1:1217242/1‑69 (MQ=255)
 aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGc                                                                    >  1:2861628/1‑69 (MQ=255)
 aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGc                                                                    >  1:2460556/1‑69 (MQ=255)
            atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc                                                          <  1:1186567/68‑1 (MQ=255)
            atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc                                                          <  1:1317808/68‑1 (MQ=255)
            atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc                                                          <  1:1370078/68‑1 (MQ=255)
            atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc                                                          <  1:694908/68‑1 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTcg                                                                 >  1:1462487/1‑59 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:541787/1‑68 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:2905603/1‑68 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:1223014/1‑68 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:896083/1‑68 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:2671785/1‑68 (MQ=255)
              tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc                                                        >  1:202159/1‑68 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:101314/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:2378216/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:2498539/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:2672245/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:2701189/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:3149250/68‑1 (MQ=255)
                 aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:840325/68‑1 (MQ=255)
                   aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  <  1:2587716/69‑1 (MQ=255)
                   aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  <  1:3255101/69‑1 (MQ=255)
                   aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  <  1:1036896/69‑1 (MQ=255)
                   aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  <  1:244504/69‑1 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:320916/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1103389/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1258592/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:769780/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1463813/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:352404/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:3243062/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:3073821/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1898872/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1921868/1‑67 (MQ=255)
                     cGCCCCGAAGCAGAAACTGTGAAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa                                                  >  1:1473443/1‑67 (MQ=255)
                               cAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc                                                     <  1:2329470/54‑1 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACa                                              >  1:350332/1‑45 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCa                                 <  1:1178281/58‑1 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCa                                 <  1:501810/58‑1 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCa                                 <  1:2098932/58‑1 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:3430465/1‑69 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:2505234/1‑69 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:3280457/1‑69 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:1838767/1‑69 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:2687182/1‑69 (MQ=255)
                                               tCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCa                      >  1:1876685/1‑69 (MQ=255)
                                                 gCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  >  1:2163576/1‑39 (MQ=255)
                                                 gCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa                                                  >  1:2983203/1‑39 (MQ=255)
                                                    gCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGacca                 >  1:1800501/1‑69 (MQ=255)
                                                    gCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGacca                 >  1:2212473/1‑69 (MQ=255)
                                                    gCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGacca                 >  1:2889336/1‑69 (MQ=255)
                                                       ttATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGaccaac               >  1:1936404/1‑68 (MQ=255)
                                                          taGCGATTAGGCTCGCGGTTGCCCAGCTCAT‑‑AGCCAGTCGAACCATTTGATGATCAGACCAACCACg           <  1:562782/67‑1 (MQ=255)
                                                           aGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACg           <  1:2167724/68‑1 (MQ=255)
                                                                   ggCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACGATCGcccc   >  1:196023/1‑68 (MQ=255)
                                                                   ggCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACGATCGCCCCg  >  1:2734896/1‑69 (MQ=255)
                                                                   ggCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACGATCGCCACg  >  1:831931/1‑69 (MQ=255)
                                                                    |                                                                   
CAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCATACAGCCAGTCGAACCATTTGATGATCAGACCAACCACGATCGCCCCG  >  minE/2418822‑2418957

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: