Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,356 T→C 13.7% 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,3560TC13.7% 1666.5 / 166.2 607pseudogene (708/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (356/148);  new base C (3/78);  total (380/227)
Fisher's exact test for biased strand distribution p-value = 1.07e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 644 total aligned reads displayed.

ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422
                                                                   |                                                                  
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:19211782/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:24935296/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:2608966/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:26585926/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:13516935/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:12505164/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    <  1:12205855/68‑1 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:18420835/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:19992032/1‑68 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   <  1:2764757/68‑1 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   >  1:32794124/1‑68 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   <  1:16416889/68‑1 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:2878509/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:25832618/1‑68 (MQ=255)
     gAGCTTGTTATTACTGGTGATTAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGa                                                               <  1:21108110/68‑1 (MQ=14)
     gAGCTTGTTATTACTGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:28488618/67‑3 (MQ=18)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACATCAATGGCAAAACCACCGTTGACGGTAAAg                                                                <  1:20826466/67‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                <  1:24546770/67‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               >  1:32445061/1‑68 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGCCAAAACCACCGTTGACGGTAAAGa                                                               <  1:1681224/68‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCaagg                                                                <  1:5164121/67‑3 (MQ=21)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                           >  1:14444006/1‑68 (MQ=255)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                             <  1:6947949/66‑1 (MQ=255)
          tGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCa                                                          <  1:33026916/68‑1 (MQ=255)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:14046151/1‑68 (MQ=255)
             tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        <  1:15793556/67‑1 (MQ=255)
             tattaAGGGTGATTATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                       <  1:31000280/68‑1 (MQ=11)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:1987000/1‑68 (MQ=255)
               ttaCTGGTGATTATGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:3045044/68‑1 (MQ=11)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:24225478/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:23479570/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:8072880/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:3122101/68‑1 (MQ=255)
               ttaCGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:1166935/68‑1 (MQ=17)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:9031902/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:33591205/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:12907782/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:28978828/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:10473551/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAccaa                                                >  1:7695700/1‑65 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:7637103/1‑68 (MQ=255)
                    ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                >  1:29521313/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                               >  1:5615392/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                               >  1:14256522/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:11943307/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:12606739/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:28220394/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:13382451/67‑1 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                <  1:3545961/67‑1 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                        <  1:31436197/68‑1 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                        <  1:7318349/68‑1 (MQ=255)
                            tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                        <  1:25555544/68‑1 (MQ=255)
                             ggcacaGTAATCAACAATGGCAAAACCACCGTTGACGGTAAAGAATCCACCGGTACGGAAATCAACgg                                       <  1:34225696/66‑1 (MQ=18)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                       <  1:915702/68‑1 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                       >  1:27193586/1‑68 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                       <  1:10840175/68‑1 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:9848923/43‑3 (MQ=14)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                               >  1:26975433/1‑59 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:26709810/1‑68 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:11321994/67‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGATATCAACGGtaa                                    >  1:29187799/1‑68 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCCCGGtaa                                    >  1:13147135/1‑68 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    <  1:9579258/68‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    >  1:4851531/1‑68 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAACGGtaa                                    >  1:5950240/1‑68 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                         >  1:3443606/1‑47 (MQ=11)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:4289134/1‑68 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:33391744/1‑67 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:22585980/1‑67 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:10029800/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:22776820/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               <  1:32080648/68‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               >  1:17641539/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGTCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:21895599/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:14765649/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:20320919/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:24785312/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:5390120/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCaa                                          >  1:17101390/1‑52 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGaaa                          >  1:21344310/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               >  1:2151409/1‑63 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGTGaaa                          >  1:14973206/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          >  1:22319440/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:6658856/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:10267324/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:32502380/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:32248151/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:31943642/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:31642919/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:29990628/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGCaat                          >  1:26263432/1‑67 (MQ=14)
                                          acaaTGGCAAAACCACCGATGACGGAAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGaaa                          >  1:16024455/1‑68 (MQ=14)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     >  1:19469776/1‑68 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:33821839/1‑68 (MQ=255)
                                                        accGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg             >  1:18899135/1‑67 (MQ=255)
                                                          cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg          >  1:17463471/1‑68 (MQ=255)
                                                           gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa         >  1:3817850/1‑68 (MQ=255)
                                                              gagggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt      <  1:632030/65‑1 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:8665156/1‑68 (MQ=255)
                                                                   tAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          <  1:16643072/43‑1 (MQ=12)
                                                                   |                                                                  
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: