Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,327,185 C→A 15.0% intergenic (‑175/+101) 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,1850CA15.0% 144.5 / 10.0 53intergenic (‑175/+101)wcaK/wzxCpredicted pyruvyl transferase/colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base C (14/31);  new base A (2/6);  total (16/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCG  >  Exported/1327120‑1327255
                                                                 |                                                                      
caaaaCCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGt                                                                    >  1:1589576/4‑70 (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:1220100/71‑1 (MQ=11)
     ccGTAGGACGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATc                                                              <  1:1226928/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)
       ccaGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAg                                                            >  1:5340866/3‑71 (MQ=9)
        tAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                           >  1:4931444/1‑71 (MQ=255)
        tAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAg                                                            <  1:3563313/70‑1 (MQ=255)
        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:3192198/1‑71 (MQ=255)
                  aaGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTc                                                 >  1:3858079/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)
                                 cATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCg                                  <  1:4400215/71‑1 (MQ=255)
                                   tCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATc                                                              <  1:478267/41‑1 (MQ=255)
                                   tCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATc                                                              <  1:4169605/41‑1 (MQ=255)
                                     cGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGaa                                     >  1:1200941/1‑64 (MQ=255)
                                     cGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAg                               >  1:2889671/1‑70 (MQ=255)
                                          ttCAGTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAg                               <  1:1815605/65‑1 (MQ=255)
                                          ttCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGca                                          <  1:1959955/54‑1 (MQ=255)
                                            cAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGAt                        <  1:608985/70‑1 (MQ=255)
                                               tGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAgg                    >  1:1385006/1‑71 (MQ=39)
                                                 ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCg                  <  1:3449394/71‑1 (MQ=39)
                                                    gATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGttt               >  1:883617/1‑71 (MQ=40)
                                                     aTGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGttt               <  1:5063851/70‑1 (MQ=255)
                                                      tGCGACGCTGACGCGTCTTATCAGGCCTACAGGTccc                                               <  1:2329526/37‑1 (MQ=255)
                                                        cGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAcg           >  1:2645926/1‑71 (MQ=255)
                                                        cGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAc            <  1:4167895/70‑1 (MQ=255)
                                                        cGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAc            >  1:1559463/1‑70 (MQ=255)
                                                           cGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGa                             <  1:1834838/50‑1 (MQ=255)
                                                           cGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGa                             <  1:5540024/50‑1 (MQ=255)
                                                            gCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                           <  1:2689230/51‑1 (MQ=255)
                                                            gCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACg                           <  1:4302795/51‑1 (MQ=255)
                                                            gCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTAcgccg        >  1:3448418/1‑70 (MQ=255)
                                                             cTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                            <  1:421379/49‑1 (MQ=255)
                                                             cTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                            <  1:2202897/49‑1 (MQ=255)
                                                             cTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGAc                            <  1:2338759/49‑1 (MQ=255)
                                                               gACGCGTCTTATCAGGCCTAAAGGTCCCGAGCACAGAAc                                    >  1:5619618/1‑39 (MQ=255)
                                                                aCGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATc    <  1:1860875/70‑1 (MQ=255)
                                                                 cgcgTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCg  <  1:592782/71‑1 (MQ=255)
                                                                 |                                                                      
CAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAGGTCCCGAGCACAGAACCGTAGGACGGATAAGGCGTTTTACGCCGCATCCG  >  Exported/1327120‑1327255

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: