Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,327,151 C→T 15.6% intergenic (‑141/+135) wcaK ← / ← wzxC predicted pyruvyl transferase/colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,327,1510CT15.6% 119.8 / 12.1 51intergenic (‑141/+135)wcaK/wzxCpredicted pyruvyl transferase/colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base C (16/27);  new base T (2/6);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAC  >  Exported/1327082‑1327221
                                                                     |                                                                      
aCGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgcc                                                                        <  1:4415337/70‑1 (MQ=255)
  gCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa                                                                     <  1:3302846/71‑1 (MQ=255)
      aCGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATc                                                                   <  1:1118124/69‑1 (MQ=255)
      aCGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCg                                                                 >  1:3713051/1‑71 (MQ=255)
       cgcgTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                <  1:207690/71‑1 (MQ=255)
        gcgTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                               <  1:2292852/71‑1 (MQ=255)
          gTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAt                                                             >  1:3404179/1‑71 (MQ=255)
           tCTTATCAGCCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                            <  1:1924439/71‑1 (MQ=255)
             ttATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCAt                                                                    >  1:4148574/1‑61 (MQ=255)
             ttATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCa                                                              >  1:1181219/1‑67 (MQ=255)
             ttATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                            >  1:2375478/1‑69 (MQ=255)
                tCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTg                                                       >  1:4129314/1‑71 (MQ=255)
                 cAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCAcgcc                                                                        >  1:2802787/1‑53 (MQ=255)
                     ccTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa                                                  <  1:1270555/71‑1 (MQ=255)
                      cTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGAt                                                 >  1:2301945/1‑71 (MQ=255)
                       tACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCa                                                          >  1:2859225/1‑61 (MQ=255)
                          aaGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGa                                             <  1:23414/71‑1 (MQ=255)
                            gACCCGAGCACAGAACCGTAGGACGGATGCGGCGTGAAcgcctt                                                                      <  1:3899978/44‑3 (MQ=11)
                               ccGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                <  1:3648665/47‑1 (MQ=255)
                                 gAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCa                                                              >  1:4311658/1‑47 (MQ=255)
                                    cacaGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAt                                                             <  1:3433983/45‑1 (MQ=255)
                                     acaGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                                                                <  1:1782674/41‑1 (MQ=255)
                                      caaaaCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGt                                  >  1:1589576/4‑70 (MQ=255)
                                      caGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                            <  1:4478341/44‑1 (MQ=255)
                                       aGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                            <  1:3749507/43‑1 (MQ=255)
                                       aGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTc                                                           <  1:2551577/44‑1 (MQ=255)
                                        gAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTc                                                           <  1:1346503/43‑1 (MQ=255)
                                        gAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCtt                               <  1:3948198/71‑1 (MQ=255)
                                        accccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCtt                               <  1:347533/68‑1 (MQ=11)
                                        accccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCt                                >  1:895719/4‑70 (MQ=11)
                                          cccGTAGGACGTATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTAt                             <  1:2099847/70‑1 (MQ=11)
                                          cccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTa                              <  1:2000872/69‑1 (MQ=11)
                                          cccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTAt                             <  1:2093025/70‑1 (MQ=11)
                                          aCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTAt                             <  1:3576398/71‑1 (MQ=255)
                                          aCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTAt                             <  1:4782157/71‑1 (MQ=255)
                                           ccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATc                            <  1:1226928/71‑1 (MQ=11)
                                           ccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATc                            <  1:1220100/71‑1 (MQ=11)
                                            cGTAGTACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATc                            <  1:1832529/70‑1 (MQ=255)
                                            cGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCa                           >  1:3286212/1‑71 (MQ=11)
                                            cGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCa                           <  1:5316478/71‑1 (MQ=255)
                                            cGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCa                           <  1:2638605/71‑1 (MQ=255)
                                             gTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGc                                                      >  1:1305785/1‑43 (MQ=38)
                                             ccaGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAg                          >  1:5340866/3‑71 (MQ=9)
                                              tAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAgg                         >  1:4931444/1‑71 (MQ=255)
                                              tAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAg                          <  1:3563313/70‑1 (MQ=255)
                                              caGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCt                                         <  1:4084969/54‑1 (MQ=2)
                                              caGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAg                          <  1:950276/69‑1 (MQ=14)
                                              caGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAg                          <  1:1941816/69‑1 (MQ=14)
                                                ggACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                        <  1:3405046/70‑1 (MQ=255)
                                                  aCGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgc                                    <  1:5170376/56‑1 (MQ=9)
                                                        aaGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTc               >  1:3858079/1‑71 (MQ=255)
                                                        aaGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTc               >  1:3192198/1‑71 (MQ=255)
                                                              ttCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCa                           <  1:5477899/53‑1 (MQ=255)
                                                                cACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAg                          >  1:3742831/1‑52 (MQ=255)
                                                                     cgcATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTccc             <  1:2828115/60‑1 (MQ=255)
                                                                     cgcATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAc  <  1:3956110/71‑1 (MQ=255)
                                                                     |                                                                      
ACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAAC  >  Exported/1327082‑1327221

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: