Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,132,044 T→A 6.8% I324I (ATA→ATT wcaD ← predicted colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,132,0440TA 6.8% 327.1 / 13.3 117I324I (ATA→ATTwcaDpredicted colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (44/65);  new base A (6/2);  total (50/67)
Fisher's exact test for biased strand distribution p-value = 7.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

Only 100 of 123 total aligned reads displayed.

CCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATC  >  AP009048/2131974‑2132114
                                                                      |                                                                      
ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCt                                                                        <  1:5874472/71‑1 (MQ=255)
  aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                      >  1:2380665/1‑71 (MQ=255)
  aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                      >  1:4674114/1‑71 (MQ=255)
      aCCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtt                                                                  >  1:7868667/1‑71 (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:4393406/1‑71 (MQ=255)
                 aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATc                                                       >  1:4903912/1‑71 (MQ=255)
                   aaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                             <  1:496102/63‑1 (MQ=255)
                    aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                                    <  1:2247013/71‑1 (MQ=255)
                    aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGc                                                    <  1:7792887/71‑1 (MQ=255)
                      cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                               <  1:7760597/58‑1 (MQ=255)
                      cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                                  >  1:1527612/1‑71 (MQ=255)
                        ccGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                 <  1:6704634/70‑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:719922/1‑41 (MQ=25)
                                   gAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGt                                                                   >  1:3871945/1‑41 (MQ=25)
                                     aaaaTAAATAATCAGCATATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTc                                      <  1:6973641/68‑1 (MQ=255)
                                     aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGtt                                              >  1:3257898/1‑60 (MQ=255)
                                     aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAata                                         <  1:1552683/65‑1 (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)
                                            ataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATACCGAATCATGc                            >  1:4985058/1‑71 (MQ=255)
                                             taatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                            >  1:156481/1‑70 (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)
                                                      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)
                                                        tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                           <  1:7500150/44‑1 (MQ=255)
                                                        tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                           <  1:1211683/44‑1 (MQ=255)
                                                        tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAatat                                        <  1:2907253/47‑1 (MQ=255)
                                                         aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                           <  1:1627961/43‑1 (MQ=255)
                                                         aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                           <  1:114641/43‑1 (MQ=255)
                                                         aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa                >  1:7709138/1‑70 (MQ=255)
                                                          cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAat                                          >  1:2190207/1‑43 (MQ=255)
                                                          cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa                >  1:7704147/1‑69 (MQ=255)
                                                          cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa                >  1:4978071/1‑69 (MQ=255)
                                                          cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa                >  1:1413066/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:6126730/71‑1 (MQ=255)
                                                            aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            <  1:5011014/71‑1 (MQ=255)
                                                            aaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            >  1:1234661/1‑71 (MQ=255)
                                                            aaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            >  1:2943038/1‑71 (MQ=255)
                                                             aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatga                               >  1:3905560/1‑51 (MQ=255)
                                                             aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            >  1:2179585/1‑70 (MQ=255)
                                                             aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg           <  1:7831205/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)
                                                             aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTg           <  1:4213394/71‑1 (MQ=255)
                                                              ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            >  1:1001865/1‑69 (MQ=255)
                                                              ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt            >  1:4473159/1‑69 (MQ=255)
                                                              ccATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa          <  1:4609703/71‑1 (MQ=255)
                                                               cATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                           <  1:3417241/37‑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:3789834/45‑1 (MQ=255)
                                                                 ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                                 <  1:3250845/45‑1 (MQ=255)
                                                                 ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACatat                     >  1:6215813/1‑57 (MQ=255)
                                                                    tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatg                                <  1:2632786/43‑1 (MQ=255)
                                                                    tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa    >  1:4656689/1‑71 (MQ=255)
                                                                     cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaa     <  1:4277059/69‑1 (MQ=255)
                                                                      tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatg                                <  1:3214572/41‑1 (MQ=255)
                                                                      tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatg                                <  1:686224/41‑1 (MQ=255)
                                                                      tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatg                                <  1:1461846/41‑1 (MQ=255)
                                                                      tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATc  >  1:3430441/1‑71 (MQ=255)
                                                                      |                                                                      
CCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATC  >  AP009048/2131974‑2132114

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: