Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 9.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 9.8% 1801.4 / 109.3 618pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (395/151);  new base G (5/54);  total (403/215)
Fisher's exact test for biased strand distribution p-value = 2.28e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 630 total aligned reads displayed.

GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426
                                                                   |                                                                   
  gtgCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:21816651/66‑1 (MQ=21)
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACATTGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:25649704/1‑68 (MQ=255)
    gcgtGTTATTACGGGTGATTATGCCACTGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:20808472/65‑1 (MQ=17)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:29691244/1‑67 (MQ=255)
     cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGAtt                                                                 >  1:19221209/1‑67 (MQ=255)
     cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGAAGGTAAAGAtt                                                                 >  1:18935581/1‑67 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:8141036/66‑1 (MQ=18)
      ttGTTATTACGGGTGATAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:22353583/66‑1 (MQ=18)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:34878837/1‑67 (MQ=255)
       tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCa                                                              <  1:35558360/68‑1 (MQ=255)
         ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:28101936/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:43586/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:23269470/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:15860161/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:7154973/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:27181917/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:12377623/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:25378921/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:26742861/68‑1 (MQ=255)
             taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         >  1:18551498/1‑67 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:36754167/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:24830834/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:6070522/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:24706939/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:26254192/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:17479846/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:36411372/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:14230340/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                       <  1:33569135/68‑1 (MQ=14)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACg                                                      >  1:14441223/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:3945086/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:23646791/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:34407825/1‑68 (MQ=255)
                gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     <  1:14887852/68‑1 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:5274042/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAGGATTCCAc                                                             >  1:3094572/1‑59 (MQ=17)
                 ggTGATAATGACACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACAGGTACGGa                                                    >  1:6278924/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:16081347/67‑1 (MQ=255)
                   tGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  >  1:22950149/1‑68 (MQ=255)
                   tGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:10479575/53‑1 (MQ=17)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAGAGATTCCACCGGTACGGaaa                                                  <  1:11766705/67‑1 (MQ=255)
                    gATAAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAAGATTCCACCGGTACGGAAAt                                                 <  1:15505569/68‑1 (MQ=14)
                        aTGCCACAGTAAACAACAATGACAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:32321235/1‑68 (MQ=255)
                        aTGCCACAGAAAACAACAAAGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCAAc                                             <  1:17188200/68‑1 (MQ=14)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:570501/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:34294914/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           <  1:35212805/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                           <  1:7823949/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:11904142/1‑68 (MQ=255)
                            cacaGTAAACAACAAAGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGta                                         <  1:28452346/68‑1 (MQ=11)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:26569548/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:8354864/1‑68 (MQ=255)
                                gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:3032612/1‑67 (MQ=255)
                                  ccacaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:34696710/63‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:13448670/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:13897870/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:20132542/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:4560469/68‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:2900224/1‑68 (MQ=255)
                                    acaacaATGCCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:36131922/68‑1 (MQ=255)
                                       acaaTGGCAACACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACTGGaaa                              <  1:24578473/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGAAATAACGGGaac                              >  1:16730079/1‑67 (MQ=14)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:11169306/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:9107364/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:6828140/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:37060161/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:3797206/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:11476115/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:24436161/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:32037432/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:28687644/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:17926804/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:1968540/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:25544504/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGGaaa                              >  1:12885804/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATATCCGGAAATAACGGGaaa                              >  1:34021888/1‑68 (MQ=17)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGTaaa                              >  1:19424253/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAg                             >  1:35989906/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:4888658/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTGGACGGTAAAGATTCCAACGGTACGGAAATCAACGGTAAGAACGGGaaa                              >  1:2224130/1‑68 (MQ=17)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:10873758/1‑68 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            >  1:8600000/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:23787001/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:32733041/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:27795036/1‑68 (MQ=255)
                                               aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       <  1:21126044/67‑1 (MQ=255)
                                               aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCa                      >  1:31053544/1‑68 (MQ=255)
                                                 aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:33310400/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:10007935/68‑1 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:16528589/1‑68 (MQ=255)
                                                   ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                  >  1:15026736/1‑68 (MQ=255)
                                                    caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACg                 >  1:962277/1‑68 (MQ=255)
                                                       cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg              <  1:22606916/68‑1 (MQ=255)
                                                             cGGTAAATATTCCACCGTTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:24120393/68‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        >  1:20470759/1‑68 (MQ=255)
                                                              ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGa       <  1:29904540/68‑1 (MQ=255)
                                                               gTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:26925400/1‑36 (MQ=255)
                                                                 aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg     >  1:29353357/1‑67 (MQ=255)
                                                                 aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGt    <  1:21310061/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  >  1:12573767/1‑68 (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: