Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 8.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,3590AG 8.8% 1660.0 / 83.8 556pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (347/147);  new base G (3/45);  total (351/205)
Fisher's exact test for biased strand distribution p-value = 7.78e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 568 total aligned reads displayed.

TGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460293‑1460426
                                                                  |                                                                   
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                    <  1:29546338/67‑1 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:10915015/1‑68 (MQ=255)
 gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:24689344/68‑1 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 >  1:25493961/1‑68 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:18557321/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:28841443/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:15214712/1‑67 (MQ=255)
   gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAt                                                                  >  1:3642240/1‑67 (MQ=255)
    cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 >  1:34848449/1‑67 (MQ=255)
    cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 >  1:32134511/1‑67 (MQ=255)
      tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCa                                                              <  1:19694145/68‑1 (MQ=255)
        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:32586277/1‑68 (MQ=255)
        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:13314628/1‑68 (MQ=255)
        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:30414428/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:28451104/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:16968118/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:5016192/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:26648798/1‑68 (MQ=255)
          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:7354707/68‑1 (MQ=255)
           ttaCGGGTGATAATGCGACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCTCCGGt                                                         <  1:11418321/68‑1 (MQ=255)
           ttaCGGGTGATAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGt                                                         <  1:20973464/68‑1 (MQ=17)
           ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:18045500/68‑1 (MQ=255)
           ttaCGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:8776173/68‑1 (MQ=255)
            taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGt                                                         >  1:15115114/1‑67 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:33206195/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:16354441/1‑68 (MQ=255)
             aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                             >  1:22429281/1‑62 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      <  1:21710509/68‑1 (MQ=255)
              cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:33826731/1‑68 (MQ=255)
                 gTGATAATGCCACTGAAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                    <  1:6181722/67‑1 (MQ=11)
                 gTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGa                                                    <  1:14575220/67‑1 (MQ=14)
                 gTGATAATGCCACAGAAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                    <  1:21793949/67‑1 (MQ=14)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:24837467/67‑1 (MQ=255)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                  <  1:6819594/67‑1 (MQ=255)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                  <  1:21243509/67‑1 (MQ=255)
                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGaaa                                                  <  1:13992016/67‑1 (MQ=255)
                   gATAAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                  <  1:20688628/67‑1 (MQ=11)
                     tAATGCCACAGTAACCAACAATGGCTAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCa                                               <  1:25093803/68‑1 (MQ=18)
                     tAACGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               <  1:10196593/68‑1 (MQ=11)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:34653583/1‑68 (MQ=255)
                       aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:25789102/1‑68 (MQ=255)
                        tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGa                                                    <  1:14269759/60‑1 (MQ=14)
                         gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGAACGCAAATCAAggg                                           >  1:16077361/1‑68 (MQ=18)
                           cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          <  1:29145257/67‑1 (MQ=255)
                            acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:22149495/1‑68 (MQ=255)
                              aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:24531821/1‑68 (MQ=255)
                              aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:33357784/1‑68 (MQ=255)
                                 aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:17973203/68‑1 (MQ=255)
                                   acaacaATGGCAAAACCACCGTTGACTGTAAAGATTACACCGGTACGTAAATCTAGGGGAATAACggg                                 >  1:12744656/1‑68 (MQ=17)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:26808249/1‑68 (MQ=255)
                                   acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:5649464/1‑68 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:1575247/1‑68 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:20152936/1‑68 (MQ=255)
                                    caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:18436804/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCTTTGACGGTAAAGAATCCACCGGTACGGAAATCAACGGTAATAACGGGaat                              >  1:27203189/1‑67 (MQ=17)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGTAAATCGCGGGCAATAACGGGatt                              >  1:32787166/1‑66 (MQ=11)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCTACGGAAATAACGGCaaa                              >  1:9740260/1‑68 (MQ=14)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:4434026/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:7987642/68‑1 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:25294732/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:22893267/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:17171312/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:27958204/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:17755771/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:14901936/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:11086790/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:1037780/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:26363666/1‑68 (MQ=255)
                                      acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGAAAATCGCGGGGAATAACGGCaat                              >  1:31579231/1‑67 (MQ=14)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:23743103/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            <  1:2266936/68‑1 (MQ=255)
                                        aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:5270079/1‑68 (MQ=255)
                                        aaTGGCAAAACCACCCTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:15430542/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:30109796/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:29880907/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:33998818/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:12297664/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGAATGTGAt                         >  1:32971108/1‑68 (MQ=255)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                        <  1:9259295/67‑1 (MQ=255)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       >  1:11986200/1‑68 (MQ=255)
                                             cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATATACGGGAATAACGGGAAAGTGATTc                       >  1:3840440/1‑68 (MQ=17)
                                                aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    >  1:4669413/1‑68 (MQ=255)
                                                aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:17445881/68‑1 (MQ=255)
                                                  ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCGGGGAATAACGGGAAAGTGATTCAGGAc                  >  1:6350804/1‑68 (MQ=14)
                                                  ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:34350101/1‑68 (MQ=255)
                                                  ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:30609757/1‑68 (MQ=255)
                                                     ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:7043947/68‑1 (MQ=255)
                                                     ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:11247747/68‑1 (MQ=255)
                                                      cGTTGACGGTAAAGATTCCGCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACgggg              >  1:23936056/1‑66 (MQ=255)
                                                      cGTTGACGGTAAAGATTCCACCGGTGCGGAAATCAACGGTAATAAGGGGAAAGTGATTCAGGACGGCg              >  1:23409743/1‑68 (MQ=17)
                                                      cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   >  1:13228422/1‑47 (MQ=255)
                                                      cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg              <  1:19081709/68‑1 (MQ=255)
                                                       gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGa             >  1:845346/1‑68 (MQ=255)
                                                       gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGa             >  1:19463709/1‑68 (MQ=255)
                                                        ttGACGGTAAAGATTCCACCGGTATGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:7661158/68‑1 (MQ=255)
                                                         tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:27555958/67‑1 (MQ=255)
                                                          gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt          <  1:5948986/68‑1 (MQ=255)
                                                           aCGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg         >  1:886074/1‑68 (MQ=255)
                                                            cGGTAAAGATTCCACCGGTACGGAAATCAACGGtaataa                                     >  1:18914897/1‑39 (MQ=255)
                                                                 aaGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   >  1:31067400/1‑36 (MQ=37)
                                                                  |                                                                   
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: