Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 7.7% 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 7.7% 2709.9 / 108.2 896pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (550/257);  new base G (18/50);  total (572/324)
Fisher's exact test for biased strand distribution p-value = 2.02e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 907 total aligned reads displayed.

GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426
                                                                   |                                                                   
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:347227/1‑68 (MQ=255)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:27964940/1‑67 (MQ=255)
     cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                                 >  1:45862141/1‑67 (MQ=255)
         ttattaCGGGTGATAATGCAACAGTAAACAACAATGTAAAAACCACCGTTGACGGAAAAGATTACAcc                                                            >  1:18898078/1‑68 (MQ=25)
          tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            <  1:47197337/67‑1 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:41726452/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:22319923/68‑1 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:9429375/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:79577/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:13134001/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:23936722/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:49391661/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTAc                                                       >  1:6776192/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:30312092/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:23633079/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:1659305/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAAGAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:29200762/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                    <  1:35478408/68‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:31523454/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:3759033/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:41066171/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:48918953/67‑1 (MQ=255)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:38210890/67‑1 (MQ=255)
                     aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                >  1:19592614/1‑68 (MQ=255)
                       aaTGCCACAGTAAACAACCATGGCAAAACCACCGTCGACGGCAAGGATTCCACCGGTACGGAAATCaa                                              <  1:49893459/68‑1 (MQ=14)
                       aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCaa                                              <  1:48240986/68‑1 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:46659005/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:30918899/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:41964946/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                                   >  1:47632237/1‑46 (MQ=14)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:20679993/1‑68 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:20905991/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           >  1:32116928/1‑68 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                           <  1:18401855/68‑1 (MQ=255)
                          gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                           <  1:20113194/68‑1 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:39063296/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:28741719/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAGGATTCCAc                                                             >  1:23320577/1‑49 (MQ=14)
                            cacaGTAAACAACAATGTCAAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          <  1:47980224/67‑1 (MQ=255)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              >  1:48645499/1‑47 (MQ=11)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGCTTCCACCGGTACGGAAATCAACGGtaa                                        <  1:40722458/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:10482175/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:7670126/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        <  1:33722825/68‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:644688/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:2301610/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGaa                                                   >  1:36712985/1‑57 (MQ=255)
                                 tAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCaa                                              <  1:50182615/58‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:11560847/68‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:14389668/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:25634247/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:37329074/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:50729524/68‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:50982341/68‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 <  1:5638223/68‑1 (MQ=255)
                                    acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                 >  1:33046026/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:2508745/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:17220593/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                                >  1:5875785/1‑68 (MQ=255)
                                     caacaaTGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggaa                                >  1:45003034/1‑66 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCGGGAAATAACGGCaaa                              >  1:22007038/1‑68 (MQ=14)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGta                                         >  1:8726457/1‑57 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTATTAACGGGaaa                              >  1:26071088/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:24261403/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:17938323/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaat                              >  1:33425536/1‑67 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:46155953/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:17492969/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:26433727/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:45538337/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:1704689/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:13731265/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:43175223/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:21665360/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:8439742/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:40113700/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:32664169/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:11155659/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGGaaa                              >  1:15508780/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGaaa                              >  1:36141727/1‑68 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            <  1:42408765/68‑1 (MQ=255)
                                          aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTg                           <  1:34143/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:12583182/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:27653360/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGaa                                                   >  1:26868395/1‑42 (MQ=255)
                                                 aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:14730333/68‑1 (MQ=255)
                                                      ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:16101900/68‑1 (MQ=255)
                                                      ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:43852464/68‑1 (MQ=255)
                                                      ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc               <  1:24671629/68‑1 (MQ=255)
                                                        gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             <  1:39650908/68‑1 (MQ=255)
                                                          ttaCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:46396632/65‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:12183272/67‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:9095999/67‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:1000708/67‑1 (MQ=255)
                                                           gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt          <  1:29759229/68‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTgg        <  1:7329231/68‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:38535415/68‑1 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGAGATCTgg        >  1:15242822/1‑68 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  <  1:44445267/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  >  1:21344489/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: