Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,356 T→C 11.0% pseudogene (708/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,3560TC11.0% 1733.5 / 124.7 597pseudogene (708/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (362/155);  new base C (2/62);  total (379/218)
Fisher's exact test for biased strand distribution p-value = 5.58e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 635 total aligned reads displayed.

ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422
                                                                   |                                                                  
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:25032569/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:2253190/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:31160199/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:6737959/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:4301450/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGaaa                                                                  >  1:6423389/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:17310898/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:22332885/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:26807407/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:32330595/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:34832336/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:4795525/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:10520161/1‑68 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               >  1:10915015/1‑68 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCaagg                                                                <  1:13154077/67‑3 (MQ=18)
     gAGCTTGATATTACGGGTGATAAAGCCACTGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                               <  1:34585114/68‑1 (MQ=14)
       gCTTGTTATTTCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:18020417/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:18557343/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:28841443/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:18443977/1‑67 (MQ=255)
        cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                             >  1:13814928/1‑67 (MQ=255)
        cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                             >  1:32134511/1‑67 (MQ=255)
         gtGTTATTACGGGTGATTATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                             <  1:28247051/65‑1 (MQ=18)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:14465278/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:26648798/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:1785088/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:21258062/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:5016192/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                      <  1:7354707/68‑1 (MQ=255)
               ttaCGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:8776173/68‑1 (MQ=255)
               ttaAGGGTGATAATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:14224836/68‑1 (MQ=14)
                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     >  1:15205328/1‑67 (MQ=255)
                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     >  1:15405705/1‑67 (MQ=255)
                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     >  1:8893177/1‑67 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:30135952/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:7604203/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:5432411/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:16749722/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCCCCGGTACg                                                  >  1:11278371/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:29132927/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:20976569/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:25692936/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:26922769/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:25381504/1‑68 (MQ=255)
                     gTGCTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                <  1:8467277/67‑1 (MQ=14)
                     gTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                <  1:26363554/67‑1 (MQ=11)
                     gTGATAAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCaagg                                                                <  1:19202079/51‑3 (MQ=14)
                       gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                               <  1:3070973/50‑1 (MQ=17)
                         tAACGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                           <  1:10196593/68‑1 (MQ=11)
                          aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                           >  1:33330911/1‑67 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:18925498/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:7020126/1‑68 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                        <  1:33124697/68‑1 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGAACGCAAATCAAggg                                       >  1:16077361/1‑68 (MQ=18)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACgg                                       <  1:25255791/68‑1 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:24821974/1‑68 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:1809150/67‑1 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGta                                     <  1:12342748/68‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    >  1:13252761/1‑68 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:21726269/1‑67 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:23070577/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               >  1:12013406/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACTGTAAAGATTACACCGGTACGTAAATCTAGGGGAATAACggg                             >  1:12744656/1‑68 (MQ=17)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:5649464/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             <  1:453099/68‑1 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTTACGGTAAAGATTCCACCGGTACGGAAATCCAGGGTAATAACggca                            >  1:23451328/1‑66 (MQ=14)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:20192498/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCTACGGAAATAACGGCaaa                          >  1:9740260/1‑68 (MQ=14)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAGCGGTAATAACGGGaaa                          >  1:17102075/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          >  1:33699743/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:15063118/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:23862581/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:9264631/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:11086790/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:11442582/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:24256720/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:7987642/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:22839385/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:20609003/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:15664239/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:29540058/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:3529887/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:3786921/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGAAAATCGCGGGGAATAACGGCaat                          >  1:31579231/1‑67 (MQ=14)
                                            aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                        >  1:24650818/1‑68 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:30437297/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:194246/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:29880907/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGAAATAACGGGAAAGTAAt                     >  1:34040936/1‑68 (MQ=255)
                                                 cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                   >  1:27926634/1‑68 (MQ=255)
                                                    aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                <  1:6256490/68‑1 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:17633564/1‑68 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:34350101/1‑68 (MQ=255)
                                                         ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc           <  1:14762054/68‑1 (MQ=255)
                                                         ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc           <  1:326618/68‑1 (MQ=255)
                                                             tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt        <  1:27555958/67‑1 (MQ=255)
                                                             tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt        <  1:21967462/67‑1 (MQ=255)
                                                             tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt        <  1:19808920/67‑1 (MQ=255)
                                                              gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt      <  1:5948986/68‑1 (MQ=255)
                                                              gACGGTAAAGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt      <  1:32833679/68‑1 (MQ=255)
                                                                   |                                                                  
ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422

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: