Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,132,042 T→A 6.8% D325V (GAC→GTC)  wcaD ← predicted colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,132,0420TA 6.8% 322.1 / 13.2 117D325V (GAC→GTC) wcaDpredicted colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (43/66);  new base A (6/2);  total (49/68)
Fisher's exact test for biased strand distribution p-value = 6.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 119 total aligned reads displayed.

CCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAA  >  AP009048/2131973‑2132112
                                                                     |                                                                      
cccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTc                                                                       <  1:4375625/71‑1 (MQ=255)
  cATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTa                                                                     <  1:3636939/71‑1 (MQ=255)
   aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                    >  1:3312240/1‑71 (MQ=255)
   aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                    >  1:4674114/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)
          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)
                    aaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                           <  1:496102/63‑1 (MQ=255)
                     aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                                  <  1:2247013/71‑1 (MQ=255)
                       cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                             <  1:7760597/58‑1 (MQ=255)
                       cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                                >  1:1527612/1‑71 (MQ=255)
                        acCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                               >  1:7390172/1‑71 (MQ=255)
                         ccGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                               <  1:6704634/70‑1 (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)
                              aaCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                                <  1:7189474/64‑1 (MQ=255)
                               aCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                         >  1:3623947/1‑70 (MQ=255)
                                 cAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAatat                                      <  1:2688317/71‑1 (MQ=255)
                                  aGGAAAAATAAATAATCAGCAGATACAAACCATTGTCt                                                                      >  1:884881/1‑38 (MQ=255)
                                  aGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                     <  1:7504274/71‑1 (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:2551503/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)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGtt                                            >  1:3257898/1‑60 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                     <  1:7081431/67‑1 (MQ=255)
                                      aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGa                                 <  1:5846635/71‑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)
                                              taatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                          >  1:156481/1‑70 (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:2314345/36‑1 (MQ=255)
                                                        aTACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                                  <  1:2336324/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)
                                                         tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                         <  1:7500150/44‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                         <  1:114641/43‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                         <  1:1627961/43‑1 (MQ=255)
                                                          aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa              >  1:7709138/1‑70 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAat                                        >  1:2190207/1‑43 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa              >  1:1440963/1‑69 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa              >  1:1413066/1‑69 (MQ=255)
                                                           cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa              >  1:4978071/1‑69 (MQ=255)
                                                           aaaaCCATAGACAATGGTTTTTCCGACATCCGAACCGTTAAATATTCCGAATGATGCGACATACTCATaa             <  1:1012166/69‑1 (MQ=37)
                                                            aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                               <  1:5520678/51‑1 (MQ=255)
                                                            aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat                   <  1:7146331/63‑1 (MQ=255)
                                                             aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt          <  1:6126730/71‑1 (MQ=255)
                                                             aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt          <  1:5011014/71‑1 (MQ=255)
                                                             aaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt          >  1:2943038/1‑71 (MQ=255)
                                                             aaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt          >  1:1234661/1‑71 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatga                             >  1:3905560/1‑51 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt          >  1:7126892/1‑70 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg         <  1:7831205/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg         <  1:4213394/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg         <  1:4513645/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg         <  1:3224193/71‑1 (MQ=255)
                                                              aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg         <  1:3551191/71‑1 (MQ=255)
                                                               ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAg           >  1:6895349/1‑68 (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)
                                                                cATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAtt                                     <  1:2990989/41‑1 (MQ=255)
                                                                 aTTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGc      <  1:7738393/71‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                               <  1:3250845/45‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                          <  1:5630100/50‑1 (MQ=255)
                                                                  ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat                   >  1:6215813/1‑57 (MQ=255)
                                                                     tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatg                              <  1:2632786/43‑1 (MQ=255)
                                                                     tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa  >  1:4656689/1‑71 (MQ=255)
                                                                     |                                                                      
CCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAA  >  AP009048/2131973‑2132112

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: