Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,356 T→C 12.4% 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,3560TC12.4% 1527.4 / 134.9 545pseudogene (708/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (312/145);  new base C (7/58);  total (337/208)
Fisher's exact test for biased strand distribution p-value = 2.69e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 577 total aligned reads displayed.

ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422
                                                                   |                                                                  
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:26744928/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:10051723/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:22429165/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    <  1:1514455/68‑1 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:7528233/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGaaa                                                                  >  1:17906606/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:16472965/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:25989568/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:23453177/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:13373838/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:16976964/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaa                                                                  >  1:3020263/1‑68 (MQ=255)
    tGAGCTTTTTATTACGGGTTATAATGCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:27950319/67‑3 (MQ=18)
     gAGCTTGTTATTACGGGTGATTATGCCACAGCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                               <  1:7133783/68‑1 (MQ=18)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               <  1:10888094/68‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:5436080/67‑3 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:5970794/67‑3 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:17259740/1‑67 (MQ=255)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                           >  1:21720717/1‑68 (MQ=255)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:8015412/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:6791830/1‑68 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:28257069/68‑1 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACGGGTAc                                                   >  1:1914188/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:20149213/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:22672357/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:18448706/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:25066699/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:8505213/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:12178296/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:8175928/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  <  1:28935411/68‑1 (MQ=255)
                   gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                 >  1:27365381/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:5245164/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:24649261/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:27604966/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:17739790/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:18225854/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGa                                                <  1:4018100/67‑1 (MQ=255)
                         tAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                           <  1:25702664/68‑1 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCTAc                                         >  1:21828277/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:11825455/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         <  1:13076915/68‑1 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:14272884/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:15670940/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:16286295/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                               >  1:29115250/1‑62 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                        <  1:21981757/68‑1 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                <  1:28554273/60‑1 (MQ=11)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                       >  1:3640503/1‑68 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                       <  1:16675901/68‑1 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATATACgg                                       >  1:453643/1‑68 (MQ=17)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:19682125/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:3209742/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:27681639/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:12747066/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:10328338/1‑68 (MQ=255)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:13906887/1‑68 (MQ=255)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:12285759/1‑68 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCaa                                          >  1:5666978/1‑59 (MQ=255)
                                    tAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGt                                                     <  1:902692/47‑1 (MQ=11)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:21481138/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:21025556/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:26582903/65‑1 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCCGGgaataa                                 >  1:9734886/1‑64 (MQ=14)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:23979035/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:14149242/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:4651420/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGAAGGTAAAGATTCCACGGGTACGAAAATAGCGGGAAATAACggca                            >  1:17557463/1‑66 (MQ=2)
                                          acaaTGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGTaaa                          >  1:6365284/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCCGGGAATAACGGCaat                          >  1:27277342/1‑67 (MQ=14)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACTGTAATAACGGGaaa                          >  1:6583995/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          >  1:10751864/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          >  1:1605944/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTaaa                          >  1:936534/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:23079464/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:11343988/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:12417483/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:27073071/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:14481133/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:4312341/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:22993016/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:29490131/68‑1 (MQ=255)
                                            aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                        >  1:11008585/1‑68 (MQ=255)
                                              tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGa                      <  1:20846599/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACTGTAAAGATTCCTCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:10587908/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGAt                     >  1:15373825/1‑68 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:29285726/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:1933815/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:26265709/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:18598068/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:13570228/68‑1 (MQ=255)
                                                    aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                <  1:24476767/68‑1 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              <  1:20704990/68‑1 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGGAAAGTGATTCAGGAc              >  1:23668166/1‑68 (MQ=255)
                                                           gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGCAAAGATATTCAGGACGGCGa         >  1:18805900/1‑68 (MQ=11)
                                                             tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt        <  1:4830591/67‑1 (MQ=255)
                                                              gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt      <  1:25319276/68‑1 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    <  1:12594810/68‑1 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:14892390/1‑68 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:3796619/1‑68 (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: