Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 10.2% pseudogene (711/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,460,3590AG10.2% 1433.9 / 87.2 503pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (292/143);  new base G (10/39);  total (309/194)
Fisher's exact test for biased strand distribution p-value = 3.17e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 516 total aligned reads displayed.

GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426
                                                                   |                                                                   
gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaaa                                                                     >  1:12394231/1‑68 (MQ=255)
  gAGCTTGTTATTACGGGTGATTATGCCACAGCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:7133783/68‑1 (MQ=18)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:2972935/1‑67 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:10617635/1‑68 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               <  1:16167758/68‑1 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               <  1:4281676/68‑1 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               <  1:26431775/68‑1 (MQ=255)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAGGATTcc                                                               >  1:552468/1‑67 (MQ=255)
          tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCg                                                           <  1:10193970/68‑1 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:5904908/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:6791830/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:19148563/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:10925007/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:24887050/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAACATTCCACCgg                                                          >  1:18605341/1‑68 (MQ=255)
            ttaCGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:2287345/68‑1 (MQ=17)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGtcc                                                       >  1:27265769/1‑66 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:18448706/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:20149213/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:28220326/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTAc                                                       >  1:23437627/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                       >  1:28962850/1‑68 (MQ=14)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACg                                                      >  1:18877673/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:27591696/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:27349734/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:21851857/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:306740/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:18225854/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGa                                                    <  1:4018100/67‑1 (MQ=255)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:6443330/67‑1 (MQ=255)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:17921827/67‑1 (MQ=255)
                     aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                >  1:10757690/1‑68 (MQ=255)
                      ttaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               <  1:10384028/66‑1 (MQ=14)
                       aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAttgc                                              <  1:12924756/68‑4 (MQ=11)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:11825455/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:14272884/1‑68 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACg                                            <  1:4769294/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:9553557/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:8085417/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:21981757/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACTGGTACGGAAATCAACgg                                           <  1:15507856/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:3640503/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:13905471/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                           <  1:21734461/68‑1 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:11022248/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:19682125/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:3209742/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAATGATTCCACCGGTACGGa                                                    <  1:1751033/58‑1 (MQ=17)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:15973376/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:902017/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:27656830/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:25794461/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:17676791/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:12285759/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGGACGGAAATCAACGGtaata                                      >  1:12071086/1‑68 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:21481138/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:21025556/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:22994846/65‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:26185073/1‑68 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:23979035/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:5250705/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:26685161/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGACGGGAATAACGGGaaa                              >  1:14064350/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:3837389/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:4511867/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:4937146/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:23079464/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:2370482/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:27073071/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:14481133/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:25602212/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:11343988/1‑68 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:8290391/1‑68 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:16340068/1‑68 (MQ=255)
                                          aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTg                           <  1:6123793/68‑1 (MQ=255)
                                           tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGa                          <  1:20846599/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACTTTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:7890414/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCTACGGTAATAACGGGAAAGTGAt                         >  1:2251615/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:1933815/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:25267515/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:28426698/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:2788739/68‑1 (MQ=255)
                                              cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       >  1:15006609/1‑68 (MQ=255)
                                                 aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:24476767/68‑1 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  <  1:20704990/68‑1 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:12514719/1‑68 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:29396625/1‑68 (MQ=255)
                                                    caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg                 <  1:17145799/68‑1 (MQ=255)
                                                      ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:1268078/68‑1 (MQ=255)
                                                       cGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGGAATAACGGGAAAGTGATTCAGGACGGCg              >  1:15355691/1‑68 (MQ=255)
                                                       cGTTGACGGTAAAGATTCCACCGGTACGGAAATACCAGGAAATAACGGGAATGTGATTCAGGACGGCg              >  1:13121812/1‑68 (MQ=11)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:13538922/67‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:29171234/67‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:9350002/67‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:16644203/68‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:14825380/68‑1 (MQ=255)
                                                                 aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg     >  1:13643528/1‑67 (MQ=255)
                                                                 aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGt    <  1:28749237/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  <  1:16778331/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  <  1:4594240/68‑1 (MQ=255)
                                                                   |                                                                   
GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426

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: