Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 8.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,3590AG 8.2% 1594.9 / 68.6 514pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (309/148);  new base G (6/35);  total (322/192)
Fisher's exact test for biased strand distribution p-value = 3.03e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 522 total aligned reads displayed.

TGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460293‑1460426
                                                                  |                                                                   
tGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                    >  1:8042649/1‑68 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                    <  1:28844640/67‑1 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                    <  1:20417441/67‑1 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:16531497/1‑68 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:1231320/1‑68 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:31629036/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:4807572/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:21899579/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGAt                                                                  >  1:29877553/1‑67 (MQ=255)
    cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 >  1:26968113/1‑67 (MQ=255)
       gTTATTACGGGTGATTATGCCACAGTAAACAACAAGGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                             <  1:8262475/68‑1 (MQ=17)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:21211972/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:7239387/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:22308390/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:10500136/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:30643665/1‑68 (MQ=255)
           ttaCTGGTGATTATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:2318303/68‑1 (MQ=11)
           ttaCTGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:1860619/68‑1 (MQ=14)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:28615935/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:9243422/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:15438300/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:18749815/68‑1 (MQ=14)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACGGGTa                                                        >  1:10050123/1‑67 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:11373954/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:14859011/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:28375291/1‑68 (MQ=255)
               gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACgg                                                     >  1:28126272/1‑68 (MQ=255)
               gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                     >  1:12730792/1‑68 (MQ=255)
               gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                     >  1:26407762/1‑68 (MQ=255)
                ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGa                                                    >  1:3724121/1‑68 (MQ=255)
                ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:20976764/1‑68 (MQ=255)
                ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGa                                                    >  1:31100896/1‑68 (MQ=255)
                 gttaTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:24382451/64‑1 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:7598309/67‑1 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGaa                                                   >  1:17993971/1‑68 (MQ=255)
                 gTGATAATGCAACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:27381045/1‑68 (MQ=255)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:29478506/67‑1 (MQ=255)
                   gATAAAGCCACAGGAAACAACAAAGGCAAAACCACCGTTGACGGGAATGATTCCACCGGTACGGaaa                                                  <  1:27233836/67‑1 (MQ=14)
                    aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                >  1:25910863/1‑68 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:30875022/1‑68 (MQ=255)
                        tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:10848878/68‑1 (MQ=255)
                        tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:24854409/68‑1 (MQ=255)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           <  1:23067826/68‑1 (MQ=255)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           <  1:25126223/68‑1 (MQ=255)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGta                                         <  1:1620023/68‑1 (MQ=255)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          <  1:24289109/67‑1 (MQ=255)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGta                                         <  1:10236047/68‑1 (MQ=255)
                            acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:23819248/68‑1 (MQ=255)
                            acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:21454646/1‑68 (MQ=255)
                            acaGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     >  1:24058030/1‑55 (MQ=11)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAATATTCCACCGGTACGGAAATCAACGGTAATAACg                                   >  1:7117591/1‑68 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:26859013/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:9593824/67‑1 (MQ=255)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:4585094/1‑68 (MQ=255)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:11691321/68‑1 (MQ=255)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:8500592/68‑1 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:19675394/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:10145924/1‑47 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCGGGAAATATCGGGaat                              >  1:20809035/1‑67 (MQ=14)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:13109741/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:22746618/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:8905862/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:23468946/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:23632498/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:29330012/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:19256216/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:29750324/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:18517956/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:13878900/68‑1 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:7307629/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:25341807/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:4152800/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAAGGTAATAACGGGAAAGt                            >  1:864002/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:1371809/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:17472546/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:30541487/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:6517406/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:6057623/68‑1 (MQ=255)
                                           ggCAAAACCACAGTTGACGCAAAAGATTCCACCGGAACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:8777954/1‑68 (MQ=14)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       >  1:31735362/1‑68 (MQ=255)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       >  1:514935/1‑68 (MQ=255)
                                                aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:31855576/68‑1 (MQ=255)
                                                  ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGAc                  >  1:8579940/1‑68 (MQ=255)
                                                  ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:16882104/1‑68 (MQ=255)
                                                   caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg                 >  1:19850043/1‑68 (MQ=255)
                                                     ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:15057365/68‑1 (MQ=255)
                                                     ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:3889059/68‑1 (MQ=255)
                                                     ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:2759833/68‑1 (MQ=255)
                                                      cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg              >  1:17011252/1‑68 (MQ=255)
                                                       gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             <  1:14203446/68‑1 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:5463603/1‑66 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:26210788/1‑66 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:30020559/67‑1 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:24606385/67‑1 (MQ=255)
                                                          gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt          <  1:620404/68‑1 (MQ=255)
                                                           aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg         >  1:11956739/1‑68 (MQ=255)
                                                           aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg         >  1:32129226/1‑68 (MQ=255)
                                                             ggTAAAGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGa       <  1:18221016/68‑1 (MQ=255)
                                                               tAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:19404977/43‑1 (MQ=12)
                                                                aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg     <  1:15838527/67‑1 (MQ=255)
                                                                  |                                                                   
TGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460293‑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: