Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,166 T→C 5.4% pseudogene (1518/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,461,1660TC 5.4% 2346.8 / 53.2 691pseudogene (1518/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (344/288);  new base C (30/6);  total (397/294)
Fisher's exact test for biased strand distribution p-value = 4.85e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 693 total aligned reads displayed.

CATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTG  >  CP009273/1461100‑1461233
                                                                  |                                                                   
cATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                                                                    <  1:40263959/68‑1 (MQ=255)
   cGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 <  1:26493703/68‑1 (MQ=255)
    ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:6646072/1‑67 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:39639/67‑1 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:11368190/67‑1 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:5304676/67‑1 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:14498852/67‑1 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:16973138/67‑1 (MQ=255)
        tCCAGTTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                            <  1:17220749/68‑1 (MQ=255)
        tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                            <  1:29781172/68‑1 (MQ=255)
         ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                           >  1:4702228/1‑68 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:9729478/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:35057023/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:7457454/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:16324647/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:10573644/1‑67 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:35859693/67‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:2625281/67‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:41665377/67‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:33948253/67‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:21136907/65‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:9216510/65‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:32127004/68‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:28565789/68‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:22767248/68‑1 (MQ=255)
                  cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  <  1:24244228/68‑1 (MQ=255)
                   ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  >  1:24885145/1‑67 (MQ=255)
                       aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                             >  1:13279/1‑68 (MQ=255)
                       aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAGGGCGGCGCCGGCACTGAGaacaac                                             >  1:18170571/1‑68 (MQ=17)
                        ccAGGCCTTCGTAAAAAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            <  1:27036408/68‑1 (MQ=18)
                        ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            >  1:32287134/1‑68 (MQ=255)
                        caaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGCGGCACCGGCACTCAGATCAACg                                            >  1:9498943/3‑68 (MQ=14)
                         cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                           <  1:20745687/68‑1 (MQ=255)
                         cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                           <  1:2184173/68‑1 (MQ=255)
                             ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                       >  1:3652046/1‑68 (MQ=255)
                                 cGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   <  1:8076729/68‑1 (MQ=255)
                                  gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:6296609/1‑67 (MQ=255)
                                  gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:32671166/1‑67 (MQ=255)
                                   tCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGcc                                 >  1:43518648/1‑68 (MQ=255)
                                    cAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGcc                                 <  1:31440904/67‑1 (MQ=255)
                                      aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                              <  1:15308928/68‑1 (MQ=255)
                                      aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                              >  1:12291076/1‑68 (MQ=255)
                                      aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                              >  1:11069652/1‑68 (MQ=255)
                                      aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCAcc                              >  1:37805199/1‑68 (MQ=11)
                                        aaCGAAGGCGAGAGCACTATCACCAGTGGCGGCACCGGCACTCAGGTCAACGGTAACGACGCCACcgc                            >  1:18531451/1‑68 (MQ=25)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:15937575/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:38549890/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:39170870/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:25292469/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:40746806/68‑1 (MQ=255)
                                         aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                           >  1:11892490/1‑68 (MQ=255)
                                          cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          <  1:19792988/68‑1 (MQ=255)
                                          cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          <  1:2038005/68‑1 (MQ=255)
                                          cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          >  1:8401069/1‑68 (MQ=255)
                                          cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGAACTCAGATCAACGGTAACGACGCCACCGCGa                          >  1:7670918/1‑68 (MQ=255)
                                           gAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGaa                         <  1:43343922/68‑1 (MQ=255)
                                           gAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGaa                         >  1:8816165/1‑68 (MQ=255)
                                             aGGCGAGAGCTCAATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       <  1:28137736/68‑1 (MQ=255)
                                              ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       >  1:17394820/1‑67 (MQ=255)
                                              ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       >  1:32975460/1‑67 (MQ=255)
                                              ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       >  1:41248815/1‑67 (MQ=255)
                                                cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:1385174/67‑1 (MQ=255)
                                                cGAGAGCACAATCACCAACGGTGGCACCGGCACGCAGATGAACGGTAACGACGCCACCGCGAAAAACa                    >  1:30728209/1‑68 (MQ=14)
                                                 gagaGCTCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                    <  1:24514041/67‑1 (MQ=255)
                                                 gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                    <  1:1028351/67‑1 (MQ=255)
                                                 gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                    <  1:37511862/67‑1 (MQ=255)
                                                     gctcTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:35432733/63‑1 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCCCCGCGAATAACAGTGGa               >  1:12899065/1‑68 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:4834171/1‑68 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:23786770/1‑68 (MQ=255)
                                                     gCACTATCACCAATGGAGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:25712149/1‑68 (MQ=255)
                                                     gCACCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:4262280/1‑68 (MQ=18)
                                                       aCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGCCGCCCCCGCGAATAACAGTGGaaa             >  1:34365263/1‑68 (MQ=255)
                                                         tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa           <  1:7163598/68‑1 (MQ=255)
                                                          aTCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa           >  1:2773551/1‑67 (MQ=255)
                                                          aTCACCAACGGTGGCACCGGCACGCAGATGAACGGTAAAGACGCCACCGCGAAAAACAGTGGaaaaa           >  1:30458025/1‑67 (MQ=14)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:39392952/68‑1 (MQ=255)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:585366/68‑1 (MQ=255)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:4515299/68‑1 (MQ=255)
                                                            cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:10045690/1‑67 (MQ=255)
                                                            cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:14998602/1‑67 (MQ=255)
                                                            cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:23398670/1‑67 (MQ=255)
                                                            cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:30971669/1‑67 (MQ=255)
                                                             aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:18845767/1‑66 (MQ=255)
                                                             aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:22263511/1‑66 (MQ=255)
                                                             aCCAATGGCGGCACCGCCACTCAGATCAACGGTAACGACGCCAACGCGAATAACAGTGGAAACAcc         >  1:27303068/1‑66 (MQ=255)
                                                             aCCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                    <  1:12597014/55‑1 (MQ=18)
                                                              ccAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          <  1:17758865/48‑1 (MQ=18)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:29354477/1‑64 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:31840858/1‑64 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:23157059/1‑64 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:42906412/1‑64 (MQ=255)
                                                               cAACGGTGGCACCGGCACGCAGATGAACGGTGACGACGCCACGGCGAATAACAGTGGAAAAACCAcgg     >  1:41937114/1‑66 (MQ=11)
                                                                aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:21490066/68‑1 (MQ=255)
                                                                aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:19102139/68‑1 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:28804437/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:38079755/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:40039620/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:33336231/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:37276740/1‑68 (MQ=255)
                                                                  |                                                                   
CATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTG  >  CP009273/1461100‑1461233

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: