Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,157 T→C 5.4% 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.4% 2303.0 / 47.7 685pseudogene (1509/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (317/311);  new base C (32/4);  total (365/320)
Fisher's exact test for biased strand distribution p-value = 3.19e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 723 total aligned reads displayed.

ATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAA  >  CP009273/1461090‑1461224
                                                                   |                                                                   
 tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:39469884/67‑1 (MQ=255)
 tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:35415149/67‑1 (MQ=255)
 tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:25800824/67‑1 (MQ=255)
 tCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACt                                                                     <  1:12982079/67‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                    >  1:9670918/1‑67 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                    >  1:1754606/1‑67 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:29114543/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:7599265/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:43342704/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:38804926/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:26254139/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:17051197/68‑1 (MQ=255)
  ccGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                   <  1:12609467/68‑1 (MQ=255)
    ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                                                                  >  1:11603920/1‑67 (MQ=255)
    ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                 <  1:9128865/68‑1 (MQ=255)
    ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                 <  1:10669473/68‑1 (MQ=255)
    ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                 <  1:19671731/68‑1 (MQ=255)
    ggAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                 <  1:42133067/68‑1 (MQ=255)
      aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:39247609/1‑68 (MQ=255)
      aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:31831373/1‑68 (MQ=255)
      aGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                               >  1:41563769/1‑68 (MQ=255)
       gTCCATCGGTATCCAGTTTGACGGCGACCAGGCCATTGTAAAAAACGAAGCCGAGAGCCCAATCACCa                                                              >  1:1124688/1‑68 (MQ=11)
       gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGCGCACTATCACCa                                                              >  1:34479460/1‑68 (MQ=255)
       gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:29609135/1‑68 (MQ=255)
       gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:16883505/1‑68 (MQ=255)
       gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:16555341/1‑68 (MQ=255)
       gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              >  1:30896587/1‑68 (MQ=255)
       gTCCATCGGTATCCAGATTGACGGCGACCACCCCATCGTCAAAAACGAAGGCGAGCGCACTATCACCa                                                              >  1:17343806/1‑68 (MQ=21)
         ccATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                              <  1:1246578/66‑1 (MQ=255)
             cGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        <  1:26493703/68‑1 (MQ=255)
              ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        >  1:30602366/1‑67 (MQ=255)
              ggTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                                                        >  1:33226651/1‑67 (MQ=255)
                 aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                     <  1:14498852/67‑1 (MQ=255)
                  tCCAGTTTGACGGCGACCAGGCCATCGTAAATAACGACGGCGAGAGCACTATCACCAATGGcggcacc                                                   <  1:20583538/68‑1 (MQ=21)
                  tCCAGATTGACGGCGACCAGGCCATCGTCAATAACTAAGGCGAGAGCACTATCACCAATGGcggcacc                                                   <  1:41840480/68‑1 (MQ=255)
                  tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcacc                                                   <  1:10135886/68‑1 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggccccg                                                  >  1:42901273/1‑68 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:15005580/1‑68 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:27940586/1‑68 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:2759611/1‑68 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:6727545/1‑68 (MQ=255)
                   ccAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg                                                  >  1:23539273/1‑68 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:35057023/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:36872225/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:37981171/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:38526317/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:39644881/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                >  1:20270176/1‑67 (MQ=255)
                      gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaacggc                                                >  1:18774877/1‑67 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:3419582/1‑67 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                              >  1:19077108/1‑67 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:25050295/65‑1 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:22380978/65‑1 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:4827172/65‑1 (MQ=255)
                        ttGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                <  1:13826725/65‑1 (MQ=255)
                          gACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCa                                           <  1:28565789/68‑1 (MQ=255)
                            cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         <  1:29830699/68‑1 (MQ=255)
                            cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         <  1:32600038/68‑1 (MQ=255)
                             ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         >  1:21011380/1‑67 (MQ=255)
                             ggCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa                                         >  1:35834123/1‑67 (MQ=255)
                               tgACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAATGGAGGCACCGGCACTCAGATCa                                      >  1:17746458/2‑68 (MQ=11)
                                 aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                    >  1:13279/1‑68 (MQ=255)
                                 aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                    >  1:10400478/1‑68 (MQ=255)
                                  ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                   <  1:778004/68‑1 (MQ=255)
                                  ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                   >  1:5218834/1‑68 (MQ=255)
                                  ccAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACg                                   <  1:23830736/68‑1 (MQ=255)
                                  caaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGCGGCACCGGCACTCAGATCAACg                                   >  1:9498943/3‑68 (MQ=14)
                                   cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                  <  1:2184173/68‑1 (MQ=255)
                                       ccATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGt                                 <  1:8290569/65‑1 (MQ=255)
                                       ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACTCAGATAAACGGTAac                              >  1:2453356/1‑68 (MQ=11)
                                            gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGCTCAACGGTAacgacg                          >  1:24795833/1‑67 (MQ=255)
                                            gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                          >  1:33769256/1‑67 (MQ=255)
                                            gTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacgacg                          >  1:32307258/1‑67 (MQ=255)
                                            gTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACACAGATCAAAGGTAacgacg                          >  1:15686342/1‑67 (MQ=11)
                                             tCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGcc                        >  1:43518648/1‑68 (MQ=255)
                                             tCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTAACGACGcc                        >  1:21416170/1‑68 (MQ=11)
                                              cAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGcc                        <  1:31440904/67‑1 (MQ=255)
                                               aaTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                      <  1:28263897/68‑1 (MQ=255)
                                                aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                     <  1:25625282/68‑1 (MQ=255)
                                                aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                     >  1:32911150/1‑68 (MQ=255)
                                                aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAACGACGCCAcc                     >  1:37805199/1‑68 (MQ=11)
                                                  aaCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                   <  1:28755799/68‑1 (MQ=255)
                                                   aCGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                  >  1:6290266/1‑68 (MQ=11)
                                                   aCGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                  >  1:1199649/1‑68 (MQ=255)
                                                   aCGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACTCAGATCa                                      >  1:16103935/1‑48 (MQ=14)
                                                    cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 >  1:40299982/1‑68 (MQ=255)
                                                    cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGa                 >  1:10870235/1‑68 (MQ=17)
                                                     gAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGaa                >  1:8816165/1‑68 (MQ=255)
                                                     gAAGGCGAGAGCACAATCACCAACGGAGGCACCGGCACACAGATAAACGGTTAAGACGCCACCGCGaa                >  1:37446543/1‑68 (MQ=14)
                                                        ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATa              >  1:23189543/1‑67 (MQ=255)
                                                          cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc            <  1:22789129/67‑1 (MQ=255)
                                                          cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc            <  1:7246141/67‑1 (MQ=255)
                                                          cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa           >  1:26961552/1‑68 (MQ=255)
                                                          cGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa           >  1:41033820/1‑68 (MQ=255)
                                                               gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg       <  1:4557805/67‑1 (MQ=255)
                                                               gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTgg       <  1:18270512/67‑1 (MQ=255)
                                                               gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa      >  1:5828508/1‑68 (MQ=255)
                                                               gCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa      >  1:24257737/1‑68 (MQ=255)
                                                                   tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa  <  1:40119809/68‑1 (MQ=255)
                                                                   tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa  <  1:36440161/68‑1 (MQ=255)
                                                                   |                                                                   
ATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAA  >  CP009273/1461090‑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: