Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,166 T→C 5.7% 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.7% 1562.2 / 37.6 461pseudogene (1518/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (208/219);  new base C (18/8);  total (234/227)
Fisher's exact test for biased strand distribution p-value = 4.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 463 total aligned reads displayed.

CATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTG  >  CP009273/1461100‑1461233
                                                                  |                                                                   
cATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                                                                    <  1:18941525/68‑1 (MQ=255)
cATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                                                                    <  1:2629189/68‑1 (MQ=255)
  tCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGACGGCGAGAGCACTATCACCAATggc                                                                  >  1:16194121/1‑68 (MQ=255)
    ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:10569751/1‑67 (MQ=255)
    ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:5011134/1‑67 (MQ=255)
    ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:471893/1‑67 (MQ=255)
    ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:23245058/1‑67 (MQ=255)
      tATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                                                               >  1:3002936/1‑67 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              >  1:16303843/1‑67 (MQ=255)
       aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              >  1:19296658/1‑67 (MQ=255)
        tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                            <  1:15955654/68‑1 (MQ=255)
        tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAACGGtggcacc                                                            <  1:23111472/68‑1 (MQ=255)
         ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                           >  1:6920292/1‑68 (MQ=255)
         ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                           >  1:22745908/1‑68 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:25485548/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:27985809/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:29251723/1‑67 (MQ=255)
            gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:4548625/1‑67 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:12686234/67‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:23378490/65‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:25973968/65‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:22362886/65‑1 (MQ=255)
              ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:10361396/65‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:7508827/68‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:28863453/68‑1 (MQ=255)
                gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                    <  1:19540658/68‑1 (MQ=255)
                  cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  <  1:15594023/68‑1 (MQ=255)
                   ggCGCCCAGGCCATCGTCATTAACGAAGGCGAGAGCTCTATCACCAATGGCGGCACCGGCACTCAGa                                                  >  1:7267032/1‑67 (MQ=18)
                   ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  >  1:22147512/1‑67 (MQ=255)
                   ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  >  1:28302238/1‑67 (MQ=255)
                        ccAGGCGGTTGTTAACAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            <  1:14007008/68‑1 (MQ=11)
                        ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            >  1:4177663/1‑68 (MQ=255)
                         cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                           <  1:24173188/68‑1 (MQ=255)
                             ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                       >  1:18517749/1‑68 (MQ=255)
                             ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGAAac                                       >  1:7915940/1‑68 (MQ=255)
                              cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                      <  1:28632437/68‑1 (MQ=255)
                                tCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgac                                    <  1:20518492/68‑1 (MQ=255)
                                  gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:10022769/1‑67 (MQ=255)
                                  gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:1725708/1‑67 (MQ=255)
                                     aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                               <  1:27460742/68‑1 (MQ=255)
                                      aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                              >  1:24568016/1‑68 (MQ=255)
                                       tAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                              >  1:14993974/1‑67 (MQ=255)
                                        aaCGACGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACCCCACcgc                            <  1:6581448/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTTTCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:13407775/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            >  1:8906642/1‑68 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:25334432/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:16878869/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:26999541/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:21968781/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:7725095/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:17852714/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:18725732/68‑1 (MQ=255)
                                        aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:1127727/68‑1 (MQ=255)
                                         aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                           >  1:25897963/1‑68 (MQ=255)
                                         aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                           >  1:10853538/1‑68 (MQ=255)
                                          cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          >  1:9815908/1‑68 (MQ=255)
                                          cGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACACAGATAAACGGTAACGACGCCACCGCGa                          >  1:24375031/1‑68 (MQ=14)
                                           gAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGaa                         >  1:2655010/1‑68 (MQ=255)
                                            aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGAAt                        >  1:17619822/1‑68 (MQ=11)
                                             aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       <  1:250893/68‑1 (MQ=255)
                                             aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       <  1:11468603/68‑1 (MQ=255)
                                              ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                       >  1:25649961/1‑67 (MQ=255)
                                                cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:28501332/67‑1 (MQ=255)
                                                cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:17118369/67‑1 (MQ=255)
                                                cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:5295795/67‑1 (MQ=255)
                                                 gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                    <  1:8198615/67‑1 (MQ=255)
                                                 gagaGCACCATCACCAACGGGGGCACCGGCACTCAGATAAACGGTTAAGACGCCACCGCGAATAACAg                   >  1:24984420/1‑68 (MQ=17)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:20235634/66‑1 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:23713102/66‑1 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:13901712/66‑1 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:21674634/66‑1 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:20010038/1‑68 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:14456382/1‑68 (MQ=255)
                                                     gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:9378950/1‑68 (MQ=255)
                                                     gCACCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:22539200/1‑68 (MQ=17)
                                                       accaTCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaa              <  1:26133850/64‑1 (MQ=18)
                                                       aCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaa             >  1:16499443/1‑68 (MQ=255)
                                                         tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa           <  1:352415/68‑1 (MQ=255)
                                                         tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa           <  1:24822808/68‑1 (MQ=255)
                                                          aTCACCAATGGCGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa           >  1:21407687/1‑67 (MQ=255)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:12411620/68‑1 (MQ=255)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:379139/68‑1 (MQ=255)
                                                           tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         <  1:13039671/68‑1 (MQ=255)
                                                             aCCAACGGTGGCACCGGCACGCAGATTAAAGGTAAAGACGCCACCGCGAAAAACAGAGGAAAAACCAc       >  1:20911027/1‑68 (MQ=11)
                                                              ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:12914000/1‑65 (MQ=255)
                                                              ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:1513160/1‑65 (MQ=255)
                                                              ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACt      <  1:4078936/68‑1 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:27811983/1‑64 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:20685727/1‑64 (MQ=255)
                                                               cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:10668495/1‑64 (MQ=255)
                                                                aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:15646136/68‑1 (MQ=255)
                                                                aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:18230847/68‑1 (MQ=255)
                                                                aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:640121/68‑1 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:12073564/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:17481951/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:5550712/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:21131570/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:2727813/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:14847473/1‑68 (MQ=255)
                                                                 aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:23966305/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: