Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,132,040 T→A 7.1% N326Y (AAT→TAT)  wcaD ← predicted colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,132,0400TA 7.1% 315.6 / 13.5 113N326Y (AAT→TAT) wcaDpredicted colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (42/63);  new base A (6/2);  total (48/65)
Fisher's exact test for biased strand distribution p-value = 6.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 118 total aligned reads displayed.

CCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGC  >  W3110S.gb/2131973‑2132108
                                                                   |                                                                    
 ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCt                                                                  <  1:5874472/71‑1 (MQ=255)
  cATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCAtt                                                                      <  1:4403130/66‑1 (MQ=255)
   aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                >  1:2380665/1‑71 (MQ=255)
   aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                >  1:3312240/1‑71 (MQ=255)
    tGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                               >  1:2265734/1‑71 (MQ=255)
       aCCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTc                                                                   <  1:2112846/64‑1 (MQ=255)
       aCCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtt                                                            >  1:7868667/1‑71 (MQ=255)
        ccACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtt                                                            <  1:292320/70‑1 (MQ=255)
         cacaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                          <  1:672559/71‑1 (MQ=255)
          acaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAACCCATTGTCTATGGtttt                                                          <  1:7782127/70‑1 (MQ=255)
             gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCg                                                      <  1:1592781/71‑1 (MQ=255)
               gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                          <  1:3725786/65‑1 (MQ=255)
                aTAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACa                                                   <  1:1251562/71‑1 (MQ=255)
                aTAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACa                                                   <  1:6821477/71‑1 (MQ=255)
                  aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATc                                                 >  1:4903912/1‑71 (MQ=255)
                  aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATc                                                 >  1:4393406/1‑71 (MQ=255)
                     aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                              <  1:2247013/71‑1 (MQ=255)
                       cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                         <  1:4094845/58‑1 (MQ=255)
                       cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                         <  1:7760597/58‑1 (MQ=255)
                       cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                            >  1:1527612/1‑71 (MQ=255)
                        acCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                           >  1:7390172/1‑71 (MQ=255)
                         ccGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCg                                          <  1:4873065/71‑1 (MQ=255)
                          cgcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                    >  1:891136/1‑60 (MQ=255)
                          cgcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACa                                                   <  1:259865/61‑1 (MQ=255)
                          cgcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGt                                         <  1:6829592/71‑1 (MQ=255)
                           gcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATTGTTTTTCCGACATCCGCACCGtt                                        <  1:6409211/71‑1 (MQ=255)
                           gcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGtt                                        <  1:3768782/71‑1 (MQ=255)
                             gAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTc                                                        <  1:4950666/53‑1 (MQ=255)
                               aCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     >  1:3623947/1‑70 (MQ=255)
                                 cAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAatat                                  <  1:2688317/71‑1 (MQ=255)
                                  aGGAAAAATAAATAATCAGCAGATACAAACCATTGTCt                                                                  >  1:884881/1‑38 (MQ=255)
                                    gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaa                                      <  1:232927/64‑1 (MQ=255)
                                    gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                               >  1:5899788/1‑71 (MQ=255)
                                    gAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGt                                                             >  1:6849481/1‑41 (MQ=25)
                                    gAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGt                                                             >  1:719922/1‑41 (MQ=25)
                                    gAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGt                                                             >  1:3871945/1‑41 (MQ=25)
                                     aaaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCg                              <  1:1183524/71‑1 (MQ=255)
                                     aaaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCg                              <  1:7813022/71‑1 (MQ=255)
                                      aaaaTAAATAATCAGCATATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTc                                <  1:6973641/68‑1 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTAt                                                                >  1:1595862/1‑36 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAata                                   <  1:4488095/65‑1 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAata                                   <  1:1552683/65‑1 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                 <  1:7081431/67‑1 (MQ=255)
                                       aaataaatAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGaa                            >  1:4209092/1‑71 (MQ=255)
                                         ataaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATc                                                 <  1:3640064/48‑1 (MQ=255)
                                         ataaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTa                                       <  1:3069071/58‑1 (MQ=255)
                                          taaataaTCAGCAGATACAAACCATTGTCTTTGGTTTTTCCGACATCCGCACCGTTa                                       <  1:3324015/57‑1 (MQ=255)
                                             ataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                      >  1:1082683/1‑71 (MQ=255)
                                             ataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATACCGAATCATGc                      >  1:4985058/1‑71 (MQ=255)
                                                   agcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGa                    >  1:7711657/1‑67 (MQ=255)
                                                    gcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat               <  1:7388549/71‑1 (MQ=255)
                                                      agATACAAACCATTGTCTATGGTTTTTCCGACATcc                                                <  1:7400778/36‑1 (MQ=255)
                                                      agATACAAACCATTTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCa            >  1:7013265/1‑71 (MQ=255)
                                                       gATACAGACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                 >  1:1020373/1‑50 (MQ=255)
                                                       gATACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                              >  1:4906023/1‑37 (MQ=255)
                                                       gATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaa                                      <  1:3844485/45‑1 (MQ=255)
                                                       gATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                 >  1:2813243/1‑50 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGCCATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCa            >  1:3159121/1‑70 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                              <  1:2336324/36‑1 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                              <  1:2314345/36‑1 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                            >  1:4032211/1‑38 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                            >  1:5865505/1‑38 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGa                    <  1:3941636/62‑1 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCa            >  1:4897633/1‑70 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCa            >  1:3168239/1‑70 (MQ=255)
                                                         tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     <  1:7500150/44‑1 (MQ=255)
                                                         tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAatat                                  <  1:2907253/47‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     <  1:114641/43‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     <  1:1627961/43‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa          >  1:7709138/1‑70 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa          >  1:1413066/1‑69 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa          >  1:4978071/1‑69 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa          >  1:1440963/1‑69 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa          >  1:7704147/1‑69 (MQ=255)
                                                           aaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaa         <  1:6946266/69‑1 (MQ=38)
                                                           aaaaCCATAGACAATGGTTTTTCCGACATCCGAACCGTTAAATATTCCGAATGATGCGACATACTCATaa         <  1:1012166/69‑1 (MQ=37)
                                                            aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                           <  1:5520678/51‑1 (MQ=255)
                                                            aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat               <  1:7146331/63‑1 (MQ=255)
                                                             aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      <  1:5011014/71‑1 (MQ=255)
                                                             aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      <  1:6126730/71‑1 (MQ=255)
                                                             aaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:1234661/1‑71 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:7126892/1‑70 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:2179585/1‑70 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg     <  1:4513645/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg     <  1:3551191/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg     <  1:4213394/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg     <  1:3224193/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg     <  1:7831205/71‑1 (MQ=255)
                                                               ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:1001865/1‑69 (MQ=255)
                                                               ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:4960397/1‑69 (MQ=255)
                                                               ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt      >  1:4473159/1‑69 (MQ=255)
                                                               ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa    <  1:4609703/71‑1 (MQ=255)
                                                                cATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     <  1:6878025/37‑1 (MQ=255)
                                                                cATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                     <  1:3417241/37‑1 (MQ=255)
                                                                cATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                 <  1:2990989/41‑1 (MQ=255)
                                                                 aTTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGc  <  1:7738393/71‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                           <  1:3789834/45‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                           <  1:3250845/45‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                      <  1:5630100/50‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat               >  1:6215813/1‑57 (MQ=255)
                                                                   |                                                                    
CCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGC  >  W3110S.gb/2131973‑2132108

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: