Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,359 A→G 8.4% 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 8.4% 1896.1 / 84.0 624pseudogene (711/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (386/171);  new base G (10/41);  total (401/223)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 633 total aligned reads displayed.

GTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460292‑1460426
                                                                   |                                                                   
gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaaa                                                                     >  1:31349839/1‑68 (MQ=255)
gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaaa                                                                     >  1:25501077/1‑68 (MQ=255)
  gagaTTGTTATTACTGGTGATAATGCCACTGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:16899260/68‑1 (MQ=14)
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   <  1:7949532/68‑1 (MQ=255)
  gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                   >  1:4327116/1‑68 (MQ=255)
    gCTTGTTATTACTGGTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:13475607/68‑1 (MQ=17)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:8417036/1‑67 (MQ=255)
    gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                                  >  1:2776198/1‑67 (MQ=255)
      ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                               >  1:398143/1‑68 (MQ=255)
         ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                            >  1:4573208/1‑68 (MQ=255)
         ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              >  1:11188491/1‑66 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:9669571/1‑68 (MQ=255)
           attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                          >  1:12937576/1‑68 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:8844469/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                         <  1:15596702/68‑1 (MQ=255)
            ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:29857378/68‑1 (MQ=255)
            ttaCGGGTGATAAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:38937377/68‑1 (MQ=14)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:9120288/1‑68 (MQ=255)
              aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                       >  1:481938/1‑68 (MQ=255)
               cGGGTGATTATGCCACAGTAAACAACAATGGCAATACCACCGTTCACGGTAAAGATTCCACCGGTACg                                                      >  1:6408543/1‑68 (MQ=21)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:36507425/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:26785127/1‑68 (MQ=255)
               cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                      >  1:20380383/1‑68 (MQ=255)
                 ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    >  1:6966751/1‑68 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGa                                                    <  1:28344816/67‑1 (MQ=255)
                  gTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:24990796/54‑1 (MQ=14)
                    gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                                  <  1:10044218/67‑1 (MQ=255)
                      tAATGCGACAGTAATCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                               <  1:30052863/68‑1 (MQ=255)
                      tAATGCCACTGTAAACAACAAGGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               <  1:32625032/68‑1 (MQ=11)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:34756474/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                             >  1:19173588/1‑68 (MQ=255)
                        aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAACGATTCCACCGGt                                                         >  1:30942710/1‑56 (MQ=17)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:14103204/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                            <  1:30547022/68‑1 (MQ=255)
                         tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCg                                                           >  1:37992413/1‑53 (MQ=14)
                         agCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:2365598/67‑1 (MQ=11)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:22428436/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          >  1:14809840/1‑68 (MQ=255)
                           ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              >  1:12316054/1‑48 (MQ=11)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGta                                         >  1:33003494/1‑68 (MQ=255)
                            cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                          <  1:13657787/67‑1 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:21650531/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:26075573/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:37212544/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:21310337/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:15968378/1‑68 (MQ=255)
                             acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:20117284/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:36820616/1‑68 (MQ=255)
                               aGTAAACAACAATGGCAAAACCACCGTAGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:30361929/1‑68 (MQ=255)
                                gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:24037566/1‑67 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACGGGTACGGAAATCAACGGTAATAACg                                   >  1:37473122/1‑68 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:2930912/65‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:9473052/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:22978388/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAAc                                    <  1:5625914/67‑1 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAACGGTAATAACg                                   >  1:3375739/1‑68 (MQ=255)
                                  aaaCAACAATGGAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:10621326/67‑1 (MQ=14)
                                     caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAAGGTAATAACGGGa                                >  1:38363337/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCGGGGAATAACGGGaat                              >  1:10304855/1‑67 (MQ=14)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCAAGGTAATAACGGGaat                              >  1:12555371/1‑67 (MQ=17)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:1511112/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:16532687/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:17498241/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:9007630/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:8246955/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:1836362/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:7539199/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:19336950/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:19700233/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:21922223/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:24235514/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:38712687/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:27657610/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:31790303/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:14513012/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:34984156/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:35680093/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              <  1:28783841/68‑1 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:35894414/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                              >  1:36300382/1‑68 (MQ=255)
                                       acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGGAATAACGGGaaa                              >  1:15267375/1‑68 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACTGTAAAGATTCCCCCGGTACTGAAATCAACTGTAATAACTGCaaatt                            >  1:6118348/1‑66 (MQ=14)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            <  1:5355829/68‑1 (MQ=255)
                                         aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                            <  1:34256547/68‑1 (MQ=255)
                                          aTGGCAAAACCACCGTTGACGGTAAAGACTCCACCGGTACGGAAATCAACGGTAATAACCGGAAAGTg                           <  1:27154505/68‑1 (MQ=17)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:39343938/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:27974289/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         >  1:35328870/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                         <  1:18870349/68‑1 (MQ=255)
                                              cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                       >  1:862274/1‑68 (MQ=255)
                                                 aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    <  1:38084286/68‑1 (MQ=255)
                                                  accaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                    >  1:11761500/1‑67 (MQ=255)
                                                        gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa             >  1:17829920/1‑68 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:34168712/67‑1 (MQ=255)
                                                          tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt            <  1:26355319/67‑1 (MQ=255)
                                                            aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg         >  1:36134312/1‑68 (MQ=255)
                                                             cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        >  1:32283295/1‑68 (MQ=255)
                                                             cGGTAAAGATTCAACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg        <  1:20716151/68‑1 (MQ=255)
                                                                 aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGt    <  1:33224199/68‑1 (MQ=255)
                                                                   aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa  <  1:28619070/68‑1 (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: