Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,157 T→C 5.9% pseudogene (1509/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,1570TC 5.9% 1650.3 / 39.0 491pseudogene (1509/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (222/226);  new base C (25/3);  total (259/232)
Fisher's exact test for biased strand distribution p-value = 2.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 521 total aligned reads displayed.

TCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAA  >  CP009273/1461091‑1461224
                                                                  |                                                                   
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:17980105/67‑1 (MQ=255)
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:21947258/67‑1 (MQ=255)
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:27535682/67‑1 (MQ=255)
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                    >  1:3242129/1‑68 (MQ=255)
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                    >  1:12434656/1‑68 (MQ=255)
tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAACGCGAGAGCACTa                                                                    >  1:25287255/1‑68 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:8200494/68‑1 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:20702284/68‑1 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:6824537/68‑1 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:18382816/68‑1 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:28004218/68‑1 (MQ=255)
 ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:30731823/68‑1 (MQ=255)
   ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                                                                  >  1:26034459/1‑67 (MQ=255)
   ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                                                                  >  1:26197323/1‑67 (MQ=255)
     aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:19966912/1‑68 (MQ=255)
     aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:17198057/1‑68 (MQ=255)
     aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:4520327/1‑68 (MQ=255)
      gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               <  1:24071642/67‑1 (MQ=255)
      gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:9129917/1‑68 (MQ=255)
      gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:13848106/1‑68 (MQ=255)
       tCCATCGGTATCCAGATTGACGGCGACCAGGCCAACGTCAATAACGAAGGCGAGAGCACTATCACCaa                                                             >  1:29301848/1‑68 (MQ=255)
        ccATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAt                                                            >  1:10961821/1‑68 (MQ=255)
          aTCGGTATCCAGATTGACGGCGCCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATg                                                           >  1:2457469/1‑67 (MQ=255)
           tCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggc                                                         >  1:3016610/1‑68 (MQ=255)
           tCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggc                                                         >  1:11730986/1‑68 (MQ=255)
             ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        >  1:8783/1‑67 (MQ=255)
             ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        >  1:26420085/1‑67 (MQ=255)
             ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        >  1:25769280/1‑67 (MQ=255)
                aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                     >  1:26584557/1‑67 (MQ=255)
                 tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                   <  1:5966073/68‑1 (MQ=255)
                  ccAGGTTGACGGCGACCAGGCGTTGGTGAATAACGAAGGCGAGAGCACAATCACCAAAGGcggcaccg                                                  >  1:10027064/1‑68 (MQ=11)
                  ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcgggaccg                                                  >  1:29086933/1‑68 (MQ=255)
                  ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:1317663/1‑68 (MQ=255)
                  ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:28978555/1‑68 (MQ=255)
                  ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:27723879/1‑68 (MQ=255)
                   cAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccgg                                                 <  1:3004605/68‑1 (MQ=255)
                     gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:16817237/1‑67 (MQ=255)
                     gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:7017764/1‑67 (MQ=255)
                     gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:162648/1‑67 (MQ=255)
                     gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:22313716/1‑67 (MQ=255)
                     gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:31891387/1‑67 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:7636789/1‑67 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              <  1:1020719/67‑1 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:16455261/1‑67 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:20572317/1‑67 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:14345014/1‑67 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:14538414/65‑1 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:23708848/65‑1 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:20833726/65‑1 (MQ=255)
                       ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:25203935/65‑1 (MQ=255)
                       ttGACGGCGACCAGGCCATAGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:2158025/65‑1 (MQ=255)
                       ttGACGGCGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:19406068/65‑1 (MQ=255)
                         gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                           <  1:28094993/68‑1 (MQ=255)
                           cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         <  1:22998386/68‑1 (MQ=255)
                            ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTTTCACCAATGGCGGCACCGGCACTCAgc                                         >  1:21369536/1‑66 (MQ=255)
                            ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCTCCAATGGCGGCACCGGCACTCAGa                                         >  1:27146508/1‑67 (MQ=255)
                            ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         >  1:23491472/1‑67 (MQ=255)
                                 ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                   <  1:29548874/68‑1 (MQ=255)
                                    ggCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACTCAGATAAAAGGtt                                >  1:17634496/1‑67 (MQ=11)
                                     gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTa                                >  1:1480395/1‑67 (MQ=255)
                                      ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGCTCAACGGTAac                              >  1:4201359/1‑68 (MQ=255)
                                      ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                 <  1:24858526/65‑1 (MQ=255)
                                      ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                              >  1:11443419/1‑68 (MQ=255)
                                         tCGTCAATAACGAAGGCGAGAGCTCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgac                           <  1:22739662/68‑1 (MQ=255)
                                          cGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgac                           >  1:28485639/1‑67 (MQ=255)
                                           gTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              <  1:24184287/47‑1 (MQ=255)
                                              aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                           <  1:7291362/47‑1 (MQ=255)
                                              aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                      <  1:16660868/68‑1 (MQ=255)
                                               aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAg                                          >  1:9805159/1‑47 (MQ=255)
                                               aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                     >  1:31681265/1‑68 (MQ=255)
                                               aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                     >  1:15970560/1‑68 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:25057731/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:6074842/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   >  1:4897216/1‑68 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:1553509/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:15774222/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:21732805/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:30989616/68‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACCATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                 <  1:5144040/54‑1 (MQ=255)
                                                 aaCGAAGGCGAGAGCACCATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:4740752/68‑1 (MQ=255)
                                                  aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                  >  1:31517382/1‑68 (MQ=255)
                                                   cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 >  1:7849844/1‑68 (MQ=255)
                                                   cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 <  1:12280437/68‑1 (MQ=255)
                                                   cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 <  1:3901109/68‑1 (MQ=255)
                                                   cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 <  1:15267962/68‑1 (MQ=255)
                                                   cGAAGGCGAGAGCACTATAAACAATGTATGCACCGTCACTCAGATCAACGGTAACGACGCCACCGCGa                 <  1:28751593/68‑1 (MQ=25)
                                                      aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa              <  1:18791184/68‑1 (MQ=255)
                                                       ggCGAGAGCACCATCACCAAAGGGGGCACCGGCACGCAGATAAAAGGTTACGACGCCACCGCGAATa              >  1:14986042/1‑67 (MQ=14)
                                                         cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc            <  1:5779126/67‑1 (MQ=255)
                                                         cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa           >  1:5449262/1‑68 (MQ=255)
                                                           agagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGCCGCCACCGCGAATAACAGt         >  1:14167316/1‑68 (MQ=255)
                                                            gagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAga         <  1:10752225/67‑2 (MQ=255)
                                                            gagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg        >  1:10459312/1‑68 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGCCGCCACCGCGAATAACAGTGGa      >  1:13664777/1‑68 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg       <  1:13489255/67‑1 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg       <  1:11877884/67‑1 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg        <  1:19245696/66‑1 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg        <  1:18989394/66‑1 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa      >  1:22631331/1‑68 (MQ=255)
                                                              gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa      >  1:25376287/1‑68 (MQ=255)
                                                                  |                                                                   
TCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAA  >  CP009273/1461091‑1461224

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: