Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,169 C→T 5.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 5.1% 1606.9 / 30.8 504pseudogene (1521/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (215/240);  new base T (17/8);  total (252/252)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 682 total aligned reads displayed.

TCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATG  >  CP009273/1461102‑1461236
                                                                   |                                                                   
  ggTATCCAGATTGACGGCGCCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATggcg                                                                    >  1:12483534/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:8652634/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:6413043/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:36919569/1‑67 (MQ=255)
  ggTATCCAGATTGACGGCGACCAGGCCATCGACAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                                    >  1:32545405/1‑67 (MQ=255)
     aTCCAGTTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 <  1:5121836/67‑1 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 <  1:27489538/67‑1 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 >  1:24423431/1‑67 (MQ=255)
     aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                 <  1:13579374/67‑1 (MQ=255)
      tCCAGATTGACGGCGACCAGGGCTTGGTCATTAACGAAGTCGAGAGCACTATCACCAATGGcggcacc                                                               <  1:15995585/68‑1 (MQ=21)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:3600286/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:23277626/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:16272656/1‑68 (MQ=255)
       ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                              >  1:29592138/1‑68 (MQ=255)
          gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            >  1:36193650/1‑67 (MQ=255)
            ttGACGGCGACCAGGCGTTTGTTAATAACGAAGGCGAGAGCTCTATCACCAATGGcggcaccggcac                                                          <  1:14669668/67‑1 (MQ=11)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCTATGGcggcaccggc                                                            <  1:19438613/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                          >  1:27996430/1‑67 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:17365834/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:81029/65‑1 (MQ=255)
            ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:23783735/65‑1 (MQ=255)
            gtGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                            <  1:23054369/64‑1 (MQ=255)
              gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                                       <  1:33666118/68‑1 (MQ=255)
              gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGCGCACTATCACCAATGGCGGCACCGGCACTCa                                                       >  1:31730354/1‑68 (MQ=255)
              gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGGGGCACCGGCACTCa                                                       >  1:8205973/1‑68 (MQ=255)
                cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     <  1:7183905/68‑1 (MQ=255)
                cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     <  1:17605511/68‑1 (MQ=255)
                 ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                                     >  1:25227376/1‑67 (MQ=255)
                     aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCatc                                                >  1:17938955/1‑66 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:15793123/68‑1 (MQ=255)
                      ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:25914811/68‑1 (MQ=255)
                      caaGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                               <  1:26721767/66‑1 (MQ=255)
                       cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                              <  1:27086829/68‑1 (MQ=255)
                         ggCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCCCCGGCACTCAGAGCAACGGTa                                            >  1:10312381/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                             <  1:30863042/65‑1 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                             <  1:10045181/65‑1 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTaaa                                          >  1:11631502/1‑67 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:9711657/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAac                                          >  1:8952123/1‑68 (MQ=255)
                           ccATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACTCAGATCAACGGTAac                                          >  1:3995558/1‑68 (MQ=17)
                            cATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacg                                         <  1:18491900/68‑1 (MQ=255)
                              tCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgac                                       <  1:19792492/68‑1 (MQ=255)
                                  cAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCACCGGTAACGACGcc                                    >  1:13299329/1‑67 (MQ=255)
                                   aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                                  <  1:32838115/68‑1 (MQ=255)
                                    aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                 >  1:20188774/1‑68 (MQ=255)
                                      aaCGAAGGCGATAGCTCCATCACCAATGGAGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:33294616/68‑1 (MQ=18)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:21849630/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:21986144/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:27337999/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               <  1:16433093/68‑1 (MQ=255)
                                      aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                               >  1:682780/1‑68 (MQ=14)
                                        cGAAGGCGAGAGCTCAATCACCAATGGGGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:34704691/68‑1 (MQ=21)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:30196235/68‑1 (MQ=255)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             >  1:28412128/1‑68 (MQ=255)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             <  1:27700359/68‑1 (MQ=255)
                                        cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                             >  1:34228404/1‑68 (MQ=255)
                                        cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGa                             >  1:20835553/1‑68 (MQ=14)
                                        cGAAGGCGAGAGCACAATCACCAAAGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGa                             >  1:30463162/1‑68 (MQ=17)
                                           aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa                          <  1:24098650/68‑1 (MQ=255)
                                              cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAAt                           >  1:28826350/1‑64 (MQ=255)
                                              cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       >  1:746152/1‑68 (MQ=255)
                                               gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       <  1:24744121/67‑1 (MQ=255)
                                               gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                       <  1:5970029/67‑1 (MQ=255)
                                                 gagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGt                     <  1:35116211/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                   <  1:4823335/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg                   <  1:5246328/67‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTg                    <  1:21589572/66‑1 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:10657356/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:4568649/1‑68 (MQ=255)
                                                   gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  >  1:30009145/1‑68 (MQ=255)
                                                     aCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                  <  1:13899516/66‑1 (MQ=255)
                                                       caTCACCAACGGTGGCACCGGCACTCAGATCAACGGTaa                                           <  1:17180563/38‑1 (MQ=11)
                                                        aTCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa              >  1:32196937/1‑67 (MQ=255)
                                                        aTCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa              >  1:2052437/1‑67 (MQ=255)
                                                         tCACCAACGGTGGCACCGGCACGCAGATAAAAGGTGAAGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:29192075/1‑68 (MQ=14)
                                                          cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:15362596/1‑67 (MQ=255)
                                                            ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:12767816/1‑65 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:22569376/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:10936404/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:3768315/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:5664708/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc            >  1:26489563/1‑64 (MQ=255)
                                                             cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTg        <  1:11171405/68‑1 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCCCTGtt      >  1:16824404/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:15252794/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:2550300/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:31398615/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      <  1:14260340/68‑1 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:11883372/1‑68 (MQ=255)
                                                               aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGtt      >  1:30631368/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGCCGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:33922093/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:31529477/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:1898726/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:35671860/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:18403517/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:14463837/1‑68 (MQ=255)
                                                                 ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGa    >  1:36595170/1‑68 (MQ=255)
                                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATg  >  1:11342860/1‑68 (MQ=255)
                                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATg  >  1:4201188/1‑68 (MQ=255)
                                                                   cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATg  >  1:22153570/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: