Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,653 G→T 8.9% pseudogene (1005/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,6530GT 8.9% 1309.1 / 8.9 488pseudogene (1005/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (203/195);  new base T (4/49);  total (212/276)
Fisher's exact test for biased strand distribution p-value = 2.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

CCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460586‑1460720
                                                                   |                                                                   
ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                     >  1:40386679/1‑68 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:36252185/68‑1 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:34979058/68‑1 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:50531235/68‑1 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:39201967/68‑1 (MQ=40)
      cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               <  1:34690328/68‑1 (MQ=39)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:44847143/1‑67 (MQ=37)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:31915454/1‑67 (MQ=37)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:53874514/1‑67 (MQ=37)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:7993693/1‑67 (MQ=37)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:11891467/1‑67 (MQ=37)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              <  1:12882684/68‑1 (MQ=37)
          gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                           >  1:54680752/1‑68 (MQ=35)
          gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                           >  1:51671922/1‑68 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                            <  1:41213850/65‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         >  1:44738419/1‑68 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:30737407/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:26136801/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:25735714/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:28303412/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:48595657/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:56795809/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:33432069/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         >  1:53337785/1‑68 (MQ=35)
              aaTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                       >  1:47777283/1‑68 (MQ=35)
                tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     <  1:2011529/68‑1 (MQ=31)
                tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     >  1:55166191/1‑68 (MQ=31)
                tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     >  1:52201542/1‑68 (MQ=31)
                        aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAAc                                             >  1:21347481/1‑68 (MQ=14)
                         cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:5951428/68‑1 (MQ=39)
                         agCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:51060142/67‑1 (MQ=11)
                          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACTGAAATCAACgg                                           <  1:15771491/68‑1 (MQ=255)
                           ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                          >  1:42636032/1‑68 (MQ=255)
                             aCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaa                                        >  1:50858540/1‑68 (MQ=255)
                              cGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaat                                       <  1:39743365/68‑1 (MQ=255)
                               ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:57714010/68‑1 (MQ=255)
                                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:42230764/1‑67 (MQ=255)
                                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGACGGtaata                                      >  1:3776837/1‑67 (MQ=255)
                                    acaacaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                 <  1:19672954/68‑1 (MQ=255)
                                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAACGGAAATAACGGa                                 >  1:45622784/1‑68 (MQ=255)
                                       tcaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTTCGGAAATCATCGGTAATAACGGaaaa                              <  1:6502145/67‑1 (MQ=14)
                                       acaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:17204923/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:10312943/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:2178679/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:14142187/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:9275023/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:3680746/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:48940897/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:10080969/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:38395649/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:27123735/68‑1 (MQ=255)
                                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:18619586/68‑1 (MQ=255)
                                          atggCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                           <  1:7425531/66‑1 (MQ=255)
                                          acGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                           >  1:16776895/1‑68 (MQ=255)
                                           cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTa                          <  1:25032380/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGGAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:45208126/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCTAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:47398920/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                 <  1:8412725/60‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:42851471/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         >  1:24613533/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         >  1:55580493/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         >  1:29306724/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:2080146/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:49296076/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:18435607/68‑1 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCTACGGTAATAACGGAAATGTTATTCAgg                    >  1:23928841/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGTAATAACGGCAAAGTTATTCAgg                    >  1:44284028/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:28889349/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:16539579/1‑67 (MQ=255)
                                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                  >  1:11120822/1‑68 (MQ=255)
                                                     accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                <  1:22822078/68‑1 (MQ=255)
                                                      ccGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:53045385/68‑1 (MQ=255)
                                                      ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:26495554/68‑1 (MQ=255)
                                                      ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:57065440/68‑1 (MQ=255)
                                                      ccGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:9583490/68‑1 (MQ=255)
                                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:17762611/1‑68 (MQ=255)
                                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:16392726/1‑68 (MQ=255)
                                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:49708041/1‑68 (MQ=255)
                                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             <  1:47002341/68‑1 (MQ=255)
                                                          tGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:35981101/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:32371058/1‑66 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:22308608/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:26926704/67‑1 (MQ=255)
                                                          tGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:31337614/67‑1 (MQ=255)
                                                          tGACGGCAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:5616896/67‑1 (MQ=255)
                                                           gACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt          <  1:1211125/68‑1 (MQ=255)
                                                           gACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt          <  1:13463944/68‑1 (MQ=255)
                                                           gACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt          <  1:38684245/68‑1 (MQ=255)
                                                            aCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg         >  1:3423539/1‑68 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCACCGGGAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:41711380/1‑68 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:52417992/1‑68 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:23473900/1‑68 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        <  1:19082410/68‑1 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:50115186/1‑68 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        <  1:14840319/68‑1 (MQ=255)
                                                               gCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAt      >  1:46536369/1‑68 (MQ=255)
                                                                 aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg     >  1:48226551/1‑67 (MQ=255)
                                                                 aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg     <  1:40018357/67‑1 (MQ=255)
                                                                 aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg     <  1:32890524/67‑1 (MQ=255)
                                                                   ggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa  >  1:43656978/1‑68 (MQ=255)
                                                                   |                                                                   
CCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460586‑1460720

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: