Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,169 C→T 5.5% 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 5.5% 2500.1 / 62.0 791pseudogene (1521/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (350/367);  new base T (34/8);  total (415/376)
Fisher's exact test for biased strand distribution p-value = 4.56e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1094 total aligned reads displayed.

TCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATG  >  CP009273/1461102‑1461236
                                                                   |                                                                   
tCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggc                                                                     >  1:51835917/1‑68 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGCGCACTATCACCAATggcg                                                                    >  1:36680620/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:5994560/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:31127967/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:25200785/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:49907698/1‑67 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 <  1:4011067/67‑1 (MQ=255)
      tCCAGTTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                               <  1:18134279/68‑1 (MQ=255)
      tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                               <  1:38433623/68‑1 (MQ=255)
      tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                               <  1:44039622/68‑1 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:27421973/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:4341308/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:48779720/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:29135710/1‑68 (MQ=255)
         aGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:48359858/1‑68 (MQ=255)
         aGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:19990054/1‑68 (MQ=255)
          gATTGCCGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:31687573/1‑67 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAGTATCAGCAATGGcggcaccggc                                                            <  1:14291295/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          <  1:25948810/67‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          <  1:23975750/67‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          <  1:12576983/67‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          >  1:56988437/1‑67 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggca                                                           <  1:32565930/66‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:16750304/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:3362325/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGACAGCACTATCACCAATGGcggcaccggcac                                                          >  1:29595888/1‑67 (MQ=255)
                cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     >  1:2057570/1‑68 (MQ=255)
                cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     <  1:47818174/68‑1 (MQ=255)
                 ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     >  1:21092695/1‑67 (MQ=255)
                 ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGAt                                                    >  1:4547099/1‑68 (MQ=255)
                    gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCaa                                                 <  1:53533127/68‑1 (MQ=255)
                      ccAGGCCTTCGTTAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:34368069/68‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:35291249/68‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:17723631/68‑1 (MQ=255)
                      caaGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAAAGGAGGCACCGGCACTCAGATAAACg                                               >  1:26214015/3‑68 (MQ=11)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTaa                                           <  1:40948669/67‑1 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:12873271/1‑68 (MQ=255)
                            cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                         <  1:9370221/68‑1 (MQ=255)
                            cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                         >  1:35999909/1‑68 (MQ=255)
                              tCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgac                                       <  1:46914323/68‑1 (MQ=255)
                               cGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                      <  1:37295255/68‑1 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                                      >  1:6813301/1‑67 (MQ=255)
                                gTCAATAACGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACTCAGATCAACGGTAacgacg                                      >  1:45858216/1‑67 (MQ=18)
                                 tCAATAACGAAGGCGAGAGCACCATCACCAAAGGTGGCACCGGCACGCAGATCAACGGTAACGACGcc                                    >  1:15498755/1‑68 (MQ=14)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:54130552/1‑68 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:53062808/1‑68 (MQ=255)
                                     tAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:14913084/1‑67 (MQ=255)
                                      aaCGAAGGCGAGAGCTCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:11645493/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:884701/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:49120753/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:49670853/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:21445244/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:13721227/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:39020198/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:13392106/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:8021410/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:28034612/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:23651941/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:12204256/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACAATCACCAACGGAGGCACCGGCACGCAGATCAACGGTAACGACGCCACcg                                >  1:38083569/1‑67 (MQ=14)
                                       aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                              >  1:46345813/1‑68 (MQ=255)
                                       aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATAAACGGTAACGACGCCACcccg                              >  1:10011125/1‑68 (MQ=255)
                                        ctaaGGCGAGATCAATATCACCAATTGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:7296740/66‑1 (MQ=255)
                                        cGAAGGCGAGAGCTCAATCACCAACGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:40733945/68‑1 (MQ=18)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:34067305/68‑1 (MQ=255)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:21979349/68‑1 (MQ=255)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:45065200/68‑1 (MQ=255)
                                              cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       >  1:48043332/1‑68 (MQ=255)
                                               gagaGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     <  1:17320177/37‑1 (MQ=255)
                                               gagaGCACCATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       <  1:18712849/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCCCCGGCACTCAGATCAACGGTAACGCCGCCACCGCGAATAACAGTGGa                  >  1:44492342/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                   <  1:10672849/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                   <  1:47325697/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:33199716/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:56380632/1‑68 (MQ=255)
                                                     accaTCACCAACGGTGGCACCGGCACTCAGATAAACGGTTACGACGCCACCGCGAATAACGGTGGaaa                >  1:24548526/4‑68 (MQ=17)
                                                      cTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaa               <  1:39415008/68‑1 (MQ=255)
                                                       tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa              <  1:4290869/68‑1 (MQ=255)
                                                         tCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            <  1:1913495/68‑1 (MQ=255)
                                                          cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATaa                         >  1:35529473/1‑54 (MQ=255)
                                                          cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:28657088/1‑67 (MQ=255)
                                                           aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:51383929/1‑66 (MQ=255)
                                                           aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:50744528/1‑66 (MQ=255)
                                                           aCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:5249088/1‑66 (MQ=255)
                                                            ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:23103031/1‑65 (MQ=255)
                                                            ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACt         <  1:50334782/68‑1 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAACCCAcgg        >  1:55629651/1‑66 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:26758370/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:57616174/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:11312893/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:2267850/1‑64 (MQ=255)
                                                              aaTGGCGGCACCGGGACGGAGATCAACGGTTAGGACTCCACCGCGAATAACAGTGGAAAAACCACTGt       <  1:54843646/68‑1 (MQ=255)
                                                              aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt       <  1:4228404/68‑1 (MQ=255)
                                                              aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt       <  1:29938600/68‑1 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:54405494/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:19479276/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:41464226/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:27070645/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:30551929/1‑68 (MQ=255)
                                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATg  >  1:52188363/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: