Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,169 C→T 6.1% pseudogene (1521/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,1690CT 6.1% 1712.2 / 49.7 547pseudogene (1521/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (270/220);  new base T (28/4);  total (320/227)
Fisher's exact test for biased strand distribution p-value = 3.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 719 total aligned reads displayed.

TCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATG  >  CP009273/1461102‑1461236
                                                                   |                                                                   
tCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggc                                                                     >  1:2294664/1‑68 (MQ=255)
 cGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    <  1:41920801/68‑1 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:4599968/1‑67 (MQ=255)
    tATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                                                                  <  1:29730137/67‑1 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 >  1:25563933/1‑67 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 <  1:11368190/67‑1 (MQ=255)
      tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                               <  1:27179615/68‑1 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggccccg                                                              >  1:42901273/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:20782282/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:34615930/1‑68 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:9496803/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:26654668/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:33233053/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:3563950/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:37339218/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:37981171/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:10573644/1‑67 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:40172251/1‑67 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          <  1:38749618/67‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          >  1:20087640/1‑67 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          <  1:43261033/67‑1 (MQ=255)
               aCGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAg                                                      >  1:25152093/1‑68 (MQ=255)
                cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAg                                                      <  1:12714984/67‑1 (MQ=255)
                cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAg                                                      <  1:16983598/67‑1 (MQ=255)
                 ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     >  1:6126286/1‑67 (MQ=255)
                 ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     >  1:23564904/1‑67 (MQ=255)
                     aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                                >  1:13279/1‑68 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                                <  1:21371455/67‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               >  1:30203112/1‑68 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:29220478/68‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:28349954/68‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:778004/68‑1 (MQ=255)
                      caaGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:39903016/66‑1 (MQ=255)
                      caaGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGGGGCACCGGCACTCAGATGAACg                                               >  1:39326686/3‑68 (MQ=11)
                         ggCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAACGGAGGCACCGGCACGCAGATCAACGGTa                                            >  1:6877126/1‑68 (MQ=11)
                          gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCaa                                                 >  1:37622312/1‑62 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCTCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:9693001/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                             <  1:39523701/65‑1 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTa                                            <  1:41636289/66‑1 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:23342577/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:22437018/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:27664679/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:1891444/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCCAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:8750724/1‑68 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGCTCAACGGTAacgacg                                      >  1:24795833/1‑67 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                      >  1:16549750/1‑67 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                      >  1:482934/1‑67 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACTCAGATCAACGGTAacgacg                                      >  1:22060713/1‑67 (MQ=14)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:14213047/1‑68 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 <  1:4367093/68‑1 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:15957868/1‑68 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:32911150/1‑68 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:11069652/1‑68 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGCCGCCACcgc                               >  1:240066/1‑68 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:40599943/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:15588720/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:25825349/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:8368266/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:14876369/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               >  1:287546/1‑68 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:35319522/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:17908593/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:2053596/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:32358141/68‑1 (MQ=255)
                                       aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                              >  1:11892490/1‑68 (MQ=255)
                                       aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                              >  1:17231877/1‑68 (MQ=255)
                                       aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGAACAACGGTAACGCCGCCACcgcg                              >  1:32249755/1‑68 (MQ=255)
                                        cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATCAACGGTAACGACGCCACCGCGa                             >  1:31629233/1‑68 (MQ=14)
                                           aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                          <  1:26352899/68‑1 (MQ=255)
                                            ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                          >  1:34167801/1‑67 (MQ=255)
                                            ggCGAGAGCACCATCACCAAAGGTGGCACCGGCACTCAGATCAAAGGTAACGACGCCACCGCGAATa                          >  1:12338925/1‑67 (MQ=17)
                                              cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                        <  1:24375628/67‑1 (MQ=255)
                                              cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACAGCGAATAACa                       >  1:43407068/1‑68 (MQ=255)
                                               gagaGCTCCATCACCAACGGGGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:7634535/59‑1 (MQ=17)
                                               gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       <  1:1028351/67‑1 (MQ=255)
                                                   gctcTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                    <  1:38118163/63‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCCCCGCGAATAACAGTGGa                  >  1:12899065/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                   <  1:18270512/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:5299356/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:5828508/1‑68 (MQ=255)
                                                       tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa              <  1:17648188/68‑1 (MQ=255)
                                                       tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa              <  1:31195755/68‑1 (MQ=255)
                                                          cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:10045690/1‑67 (MQ=255)
                                                            ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:30870731/1‑65 (MQ=255)
                                                            ccAACGGTGGCACCGGCACGCAGATAAACGGTTACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:19360794/1‑65 (MQ=17)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAaacaacgg        >  1:13825255/1‑66 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:35945546/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:32017089/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTg        <  1:9200059/68‑1 (MQ=255)
                                                              aaTGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCTAATAACAGTGGAAAAACCACTGt       <  1:15113742/68‑1 (MQ=255)
                                                              aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt       <  1:36586411/68‑1 (MQ=255)
                                                              aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt       <  1:27962136/68‑1 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:18576781/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:6462630/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:7431846/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:41116151/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:21694432/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:42461336/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:336433/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:30633270/1‑68 (MQ=255)
                                                                   |                                                                   
TCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATG  >  CP009273/1461102‑1461236

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: