Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,166 T→C 5.2% 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.2% 3474.5 / 73.6 1039pseudogene (1518/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (474/477);  new base C (37/16);  total (546/493)
Fisher's exact test for biased strand distribution p-value = 4.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1045 total aligned reads displayed.

CCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTG  >  CP009273/1461099‑1461233
                                                                   |                                                                   
  aTCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                                                                    >  1:21175422/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:29280127/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:28350755/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:5994560/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:25436960/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:31567979/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:45706901/1‑67 (MQ=255)
     ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                 >  1:51861422/1‑67 (MQ=255)
        aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCTCCAATGGcggca                                                              <  1:30929499/67‑1 (MQ=255)
        aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              <  1:40580552/67‑1 (MQ=255)
        aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                              >  1:54234158/1‑67 (MQ=255)
         tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                            <  1:26081516/68‑1 (MQ=255)
          ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                           >  1:39935871/1‑68 (MQ=255)
            aGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:11600248/1‑68 (MQ=255)
             gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         >  1:34803860/1‑67 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       >  1:32103324/1‑67 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       >  1:29739928/1‑67 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       >  1:56988437/1‑67 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       >  1:43102277/1‑67 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:19001853/65‑1 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                         <  1:33495895/65‑1 (MQ=255)
               ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt                                                      <  1:23303756/68‑1 (MQ=255)
               tcgACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                       <  1:26286932/65‑1 (MQ=11)
                 gACGGCGACCAGCGCAGCGTCTATTTCGATGTCTTGAGCACTATACACACCGGCGGCACCGGCACTCa                                                    <  1:41269431/68‑1 (MQ=255)
                   cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  <  1:30275815/68‑1 (MQ=255)
                    ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                  >  1:46506378/1‑67 (MQ=255)
                         ccAGGCCTTCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            <  1:6765245/68‑1 (MQ=255)
                         ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                            >  1:37086297/1‑68 (MQ=255)
                            ggCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGCGGCACCGGCACTCAGATCAACGGTa                                         >  1:48921162/1‑68 (MQ=17)
                            ggCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACGCAGATCAACGGTa                                         >  1:19381332/1‑68 (MQ=14)
                              cgaTCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                          <  1:27627913/63‑1 (MQ=255)
                              ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                          <  1:54969640/65‑1 (MQ=255)
                              ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                       >  1:14150170/1‑68 (MQ=255)
                              ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                       >  1:42213075/1‑68 (MQ=255)
                              ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                       >  1:39217760/1‑68 (MQ=255)
                              ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGAGGCACCGGCACTCAGATCAACGGTAac                                       >  1:29023134/1‑68 (MQ=255)
                               cATCGTCAATAACGAAGTCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                      <  1:49753304/68‑1 (MQ=255)
                               cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                      <  1:43811852/68‑1 (MQ=255)
                               cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                      <  1:9370221/68‑1 (MQ=255)
                               cATCGTCAATAACGAAGGCGAGAGCACGATCACCAATGGCGGCACCGGCACTCa                                                    >  1:15467451/1‑54 (MQ=255)
                                  cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAacgac                                    >  1:52278000/1‑67 (MQ=11)
                                   gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:27201113/1‑67 (MQ=255)
                                   gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:14895905/1‑67 (MQ=255)
                                   gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                   >  1:16769584/1‑67 (MQ=255)
                                     cAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGcc                                 <  1:22857880/67‑1 (MQ=255)
                                      aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                               <  1:21468984/68‑1 (MQ=255)
                                       aTAACGAAGGCGAGAGCACCATCACCAAAGGGGGCACCGGCACTCAGATAAACGGTAACGACGCCAcc                              >  1:22266948/1‑68 (MQ=17)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:4699913/68‑1 (MQ=255)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            >  1:12532604/1‑68 (MQ=255)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            >  1:13223626/1‑68 (MQ=255)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:53198711/68‑1 (MQ=255)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:34233625/68‑1 (MQ=255)
                                         aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                            <  1:23704465/68‑1 (MQ=255)
                                         aaCGAAGGCGAGAGCACCATCACCAAAGGAGGCACCGGCACACAGATAAACGGTAACGACGCCACcgc                            >  1:48208242/1‑68 (MQ=14)
                                          aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATAAACGGTAACGACGCCACcccg                           >  1:10011125/1‑68 (MQ=255)
                                          aCGAAGGCGAGAGCACAATCACCAACGGAGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcc                           >  1:13615301/1‑67 (MQ=18)
                                           cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          <  1:36520315/68‑1 (MQ=255)
                                           cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                          <  1:12145088/68‑1 (MQ=255)
                                           cGAAGGCGAGAGCACTATCACCAATGGCCGCACCGGCACACAGATCAACGGTAACGACGCCACCGCGa                          >  1:48525733/1‑68 (MQ=25)
                                           cGAAGGCGAGAGCACTATCAACAAAGGCGGAACCGGCACTCAGATCAACGGTAACGACGACAACGCGa                          >  1:49312152/1‑68 (MQ=21)
                                                 cGAGAGCTCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:54456485/67‑1 (MQ=255)
                                                 cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                     <  1:47904179/67‑1 (MQ=255)
                                                 cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATGAACGGTAAAGACGCCACCGCGAATAACa                    >  1:7794975/1‑68 (MQ=11)
                                                   agagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGt                  >  1:3795207/1‑68 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:35850772/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:17441475/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:45976764/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:50945278/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:24726416/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:20819543/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                <  1:10672849/67‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:53610587/66‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                 <  1:38906028/66‑1 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:35448364/1‑68 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:56380632/1‑68 (MQ=255)
                                                      gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa               >  1:4531685/1‑68 (MQ=255)
                                                      gCACCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCGAATAACAGTGGa               >  1:49984868/1‑68 (MQ=11)
                                                           aTCACCAATGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAc          <  1:38253398/68‑1 (MQ=255)
                                                            tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:1988065/1‑68 (MQ=255)
                                                             cACCAATGGCGGCACCGGCACTCAGTTCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:12804206/1‑67 (MQ=255)
                                                             cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:20244152/1‑67 (MQ=255)
                                                             cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:44236132/1‑67 (MQ=255)
                                                             cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCa        <  1:51490133/68‑1 (MQ=255)
                                                              aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:47605681/1‑66 (MQ=255)
                                                              aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:50424894/1‑66 (MQ=255)
                                                              aCCAACGGTGGCACCGGCACGCAGATAAACGGTAAAGACGCCACCGCGAATAACAGTGGAAAAACCAc       >  1:33063766/1‑68 (MQ=17)
                                                               ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACt      <  1:44778670/68‑1 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAACCCAcgg     >  1:55629651/1‑66 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:40886601/1‑64 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:24238200/1‑64 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:42547607/1‑64 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:57616174/1‑64 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:33904674/1‑64 (MQ=255)
                                                                cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:52576972/1‑64 (MQ=255)
                                                                cAATGGCGGAACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc         >  1:8794801/1‑64 (MQ=255)
                                                                 aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:11224451/68‑1 (MQ=255)
                                                                 aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:7294080/68‑1 (MQ=255)
                                                                 aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt    <  1:26774413/68‑1 (MQ=255)
                                                                  aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:18007716/1‑68 (MQ=255)
                                                                  aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt   >  1:1858846/1‑68 (MQ=255)
                                                                   |                                                                   
CCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTG  >  CP009273/1461099‑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: