Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 9.6% 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 9.6% 2133.7 / 119.7 734pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (455/186);  new base G (8/60);  total (473/261)
Fisher's exact test for biased strand distribution p-value = 6.80e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 747 total aligned reads displayed.

GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426
                                                                   |                                                                   
  gAGCTTGTTATTACTGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                    <  1:17194610/67‑1 (MQ=255)
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:24630641/1‑68 (MQ=255)
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   <  1:38089833/68‑1 (MQ=255)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAATGAt                                                                  >  1:16054215/1‑67 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 <  1:7036046/66‑1 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:36002057/1‑68 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:16571991/1‑68 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTcc                                                               >  1:26072277/1‑68 (MQ=255)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:36566899/1‑67 (MQ=255)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCa                                                              <  1:870072/68‑1 (MQ=255)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:29000517/1‑67 (MQ=255)
         ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:39075987/1‑68 (MQ=255)
          tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                            <  1:5883923/67‑1 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:34821029/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:9421393/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:33649747/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:37342265/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:16983244/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:11743961/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:23344715/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:20031608/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:25669977/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:11328638/68‑1 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAACCCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:37860276/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:15072360/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:5235832/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       <  1:35167985/68‑1 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:1519468/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTAc                                                       >  1:13263697/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGGTGACGGTAAAGATTCCACCGGTAc                                                       >  1:2357274/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:17772573/1‑68 (MQ=255)
                gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                     >  1:15522951/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:32491895/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:20301365/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:6396552/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:4306058/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:24337005/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                                   >  1:18524043/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                    <  1:35845246/67‑1 (MQ=255)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                  <  1:7996423/67‑1 (MQ=255)
                      tAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               <  1:12685095/68‑1 (MQ=11)
                      tAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                >  1:3467984/1‑67 (MQ=255)
                      tAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                <  1:30141544/67‑1 (MQ=11)
                       aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                               >  1:14815076/1‑67 (MQ=255)
                        aTGCCACAGTAACCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:24048360/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:33371485/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:34292290/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:33023754/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTACACCGGTACGGAAATCAACgg                                           >  1:32209499/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGt                                          >  1:25376689/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:37576364/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:13955100/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:10006793/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:30474093/1‑68 (MQ=255)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCCGGta                                         >  1:29742391/1‑68 (MQ=255)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGta                                         >  1:38392924/1‑68 (MQ=255)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGta                                         <  1:35469765/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtga                                        >  1:2935826/1‑66 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:27783928/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:24917705/1‑68 (MQ=255)
                                gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:21431495/1‑67 (MQ=255)
                                gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACAGAAATCAACGGtaata                                      >  1:14113322/1‑67 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:11052906/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:10574089/65‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:9649172/1‑68 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:723932/1‑68 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:22436508/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:18246209/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:19646065/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACTGAAATCAACGGTAATAACGGGaaa                              >  1:30690515/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATATCGGGaaa                              >  1:28287660/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:33941921/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:31979695/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:31300725/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:31056368/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:4562502/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:23214454/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:26424865/68‑1 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGAAATAACGGGAAAGt                            >  1:5701019/1‑68 (MQ=255)
                                          aTGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                           <  1:12228725/68‑2 (MQ=11)
                                          aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTg                           <  1:31795639/68‑1 (MQ=255)
                                           tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGa                          <  1:6095761/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACTGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:26621337/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:27116469/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:3267158/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:13653441/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:25799823/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:38660260/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:20905714/1‑68 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:19207892/1‑68 (MQ=255)
                                                      ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:4874499/68‑1 (MQ=255)
                                                      ccGTTGACGGTAAAGAATCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:28477865/68‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:21983015/1‑66 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:14427513/67‑1 (MQ=255)
                                                           gACTGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGTAAGTGATTCAGGACGGCGATCt          <  1:37000524/68‑1 (MQ=255)
                                                            aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTg         >  1:18640758/1‑68 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        >  1:9143282/1‑68 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:22925533/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  >  1:27450019/1‑68 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  <  1:750624/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: