Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 11.8% 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,3590AG11.8% 1609.8 / 121.7 570pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (340/148);  new base G (7/58);  total (354/216)
Fisher's exact test for biased strand distribution p-value = 2.43e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 581 total aligned reads displayed.

TGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460293‑1460426
                                                                  |                                                                   
 gAGCTTGTTATTACGGGTGATCATGCCACTGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:14426498/68‑1 (MQ=18)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   <  1:20801844/68‑1 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:14073565/1‑68 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:26850920/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGAtt                                                                 >  1:25355505/1‑68 (MQ=255)
    cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGAtt                                                                 >  1:34060391/1‑67 (MQ=255)
     ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:11206683/1‑68 (MQ=255)
     ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:13294003/1‑68 (MQ=255)
     ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               <  1:16702148/68‑1 (MQ=255)
        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:29567341/1‑68 (MQ=255)
         tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCg                                                           >  1:9388212/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:4757522/68‑1 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:14828674/68‑1 (MQ=17)
           ttaCTGGTGATAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:17250680/68‑1 (MQ=14)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:5533081/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:13427449/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:2148899/68‑1 (MQ=255)
            taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         >  1:17523428/1‑67 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTAc                                                       >  1:19533973/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:12907782/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:5089109/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:33591205/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:30618557/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:31506002/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:8501632/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:21067635/1‑68 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:662190/1‑68 (MQ=255)
               gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                     >  1:19879944/1‑68 (MQ=255)
                ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:31321076/1‑68 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:3538349/1‑68 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:28220394/67‑1 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:26059994/67‑1 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:12606739/67‑1 (MQ=255)
                 gTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                    <  1:24451130/67‑1 (MQ=14)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:29009901/67‑1 (MQ=255)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:8962456/67‑1 (MQ=255)
                      aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                               >  1:22542138/1‑67 (MQ=255)
                      aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGGACGGAAATCa                                               >  1:10462071/1‑67 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAtt                                                >  1:32871274/1‑64 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:16488989/1‑68 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:1153826/1‑68 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              >  1:18944475/1‑51 (MQ=14)
                        tGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:34245147/68‑1 (MQ=11)
                        tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:28803314/68‑1 (MQ=255)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:19968981/1‑68 (MQ=255)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           <  1:32407718/68‑1 (MQ=255)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:27193586/1‑68 (MQ=255)
                          ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:26975433/1‑59 (MQ=255)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          <  1:13113577/67‑1 (MQ=255)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGt                                          <  1:7489614/67‑1 (MQ=255)
                            acaGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                             >  1:3443606/1‑47 (MQ=11)
                              aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:4289134/1‑68 (MQ=255)
                               gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGACGGaaata                                      >  1:26297434/1‑67 (MQ=255)
                               gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:24654305/1‑67 (MQ=255)
                               gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:33391744/1‑67 (MQ=255)
                                tAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACg                                            <  1:5747220/60‑1 (MQ=255)
                                tAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGa                                                    <  1:10805367/52‑1 (MQ=11)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:29336180/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:10029800/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:22776820/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:21335158/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:7293826/65‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:30171273/67‑1 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAGAAGGGAATAACg                                   >  1:7672604/1‑68 (MQ=14)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGAATCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:16057998/68‑1 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGTCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:21895599/1‑68 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:11856914/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:26638731/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:12052543/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:32502380/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:19236568/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:32053722/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:31954478/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:3100782/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:23791709/68‑1 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:28449322/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGGaat                              >  1:22160509/1‑67 (MQ=255)
                                      acaaTGGCAAAAACACCGTTGACGGTAAAGATTCCACCGGTACGGAAAGCAACGGTAATAAAGGGaaa                              >  1:32212297/1‑68 (MQ=17)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTaaa                              >  1:9813702/1‑66 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:27565746/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            <  1:14271914/68‑1 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:7940757/1‑68 (MQ=255)
                                           tgCAAAACCACCGGTGTAGGTAAATTTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:31798607/67‑1 (MQ=17)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:18511893/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:12494073/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:15197114/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:15888870/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:525695/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:25692095/68‑1 (MQ=255)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATATACGGTAATAACGGGAAAGTGATTc                       >  1:9117216/1‑68 (MQ=255)
                                                      cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAGGGAAATAACGGCAAAGTTATTCAGGAGGGCg              >  1:16851821/1‑68 (MQ=11)
                                                       gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:6936391/1‑68 (MQ=255)
                                                       gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:3817850/1‑68 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:61895/67‑1 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:2480893/67‑1 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:10250787/67‑1 (MQ=255)
                                                          taCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATGCAGGACGGCGATCt          <  1:17739057/67‑1 (MQ=255)
                                                                aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGt    <  1:330335/68‑1 (MQ=255)
                                                                  aGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   >  1:23390809/1‑35 (MQ=33)
                                                                  aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  >  1:20277395/1‑68 (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: