Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,356 T→C 9.6% 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,3560TC 9.6% 3212.4 / 196.4 1069pseudogene (708/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (640/301);  new base C (8/92);  total (673/396)
Fisher's exact test for biased strand distribution p-value = 6.32e-33
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 1142 total aligned reads displayed.

ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422
                                                                    |                                                                  
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:9395975/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:8707497/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:36695641/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:20548595/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:49585625/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    <  1:51861916/68‑1 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:1672859/1‑68 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   >  1:46308985/1‑68 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   <  1:4676355/68‑1 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAGAACCACCGTTGACGGTaa                                                                  >  1:19399759/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:32778735/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:19935736/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:18853037/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:32158382/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGACAAAACCACCGTTGACGGTaa                                                                  >  1:56328792/1‑68 (MQ=255)
     gAGCTTGTTATTACGGGTGATTATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:28733257/67‑3 (MQ=21)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                <  1:25715228/67‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGa                                                               <  1:31450956/68‑1 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:50850315/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:17628381/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:48871409/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:55747409/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:38873916/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:24311562/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGAtt                                                             <  1:13972867/68‑1 (MQ=18)
         ttGTTATTACGGGTGATTAAGCCACTGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGAtt                                                             <  1:48366431/66‑1 (MQ=17)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                             <  1:22714611/66‑1 (MQ=255)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:37340015/1‑68 (MQ=255)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:3482484/1‑68 (MQ=255)
             tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        <  1:41055291/67‑1 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:55305181/1‑68 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                      <  1:27825134/68‑1 (MQ=14)
               ttaCGGGTGATTATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:10044789/68‑1 (MQ=17)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:42298913/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:57320230/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:30832658/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:15929349/68‑1 (MQ=255)
               ttaCGGGTGATAAAGCCACAGCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:43820249/68‑1 (MQ=14)
                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     >  1:6267404/1‑67 (MQ=255)
                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACAACCGTTGACGGTAAAGATTCCACCGGt                                                     >  1:1198336/1‑67 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:45830526/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:30041593/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:47313809/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:22375135/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACAACCGTTGACGGTAAAGATTCGACCGGTa                                                    >  1:23429638/1‑67 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:29669785/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                               >  1:35409408/1‑68 (MQ=255)
                       gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaaa                                              <  1:1043382/67‑1 (MQ=255)
                       gATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              <  1:21433056/67‑1 (MQ=14)
                         tAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                           >  1:52768884/1‑68 (MQ=255)
                          aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                           >  1:53190458/1‑67 (MQ=255)
                          aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                           >  1:40174292/1‑67 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACGGTTGACGGAAAAGATTCCACCGGTACCCAAATccag                                         >  1:7252980/1‑65 (MQ=17)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         <  1:2184774/68‑1 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:27172635/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCAGAGTGTTATGGAAATATTCCACCCCTCCGGAAATCAAAgga                                      >  1:28219113/1‑67 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACTGTAAAGATTCCACCGGTACGGAATTCAACGGt                                      >  1:5411300/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCaa                                          >  1:41997845/1‑64 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:25414481/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:44299491/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:17250624/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACAACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGt                                      >  1:47768753/1‑68 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:7662917/67‑1 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:4273115/67‑1 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:2082959/67‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAGGGgaa                                    >  1:1087155/1‑68 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    <  1:40759340/68‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    >  1:43277602/1‑68 (MQ=255)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:39705151/1‑68 (MQ=255)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:26135919/1‑68 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCTACCGGTACGGAAATCAACGGtaata                                  >  1:33484410/1‑67 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACAGGTACGGAAATCAACGGtaata                                  >  1:5239063/1‑67 (MQ=255)
                                    tAAACAACAATGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCa                                           <  1:35967581/57‑1 (MQ=11)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               >  1:45477428/1‑68 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               <  1:51325147/68‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:4146352/64‑1 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:38122857/1‑68 (MQ=255)
                                          tcaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:12697573/2‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                          <  1:31794716/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGACGGAAATAACGGCaaa                          >  1:54756847/1‑68 (MQ=14)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaat                          >  1:49234342/1‑67 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:6946659/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:16333429/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:13073422/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAAAGGGAATAACGGGaaa                          >  1:45148804/1‑68 (MQ=17)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGCAAATCAACGGTAATAACGGCaaa                          >  1:342754/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACAGTTGACGGTAATGATTAAACCCGTACGCAAATCAACGGTAATAACGGGaaa                          >  1:8130842/1‑68 (MQ=14)
                                          acaaCGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:57822109/1‑68 (MQ=255)
                                              tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGa                      <  1:14638337/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:15394655/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:1230833/68‑1 (MQ=255)
                                                 cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCaa                                          >  1:11367576/1‑45 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:50070979/1‑68 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:24874290/1‑68 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGGACGGAAATCTCCGGTAATATCGGGATAGTGATTCAGGAc              >  1:44628040/1‑68 (MQ=17)
                                                       caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg             >  1:26757897/1‑68 (MQ=255)
                                                               gACGGTAAAGATTCCACCGGTACGGAAATCACCGGGAATTACGGGATTGTGATTCAGGAAGgaga         >  1:1818625/1‑65 (MQ=17)
                                                                aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg     >  1:13464749/1‑68 (MQ=255)
                                                                 cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:20376988/1‑68 (MQ=255)
                                                                  ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGa   <  1:5370528/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: