Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 7.2% pseudogene (945/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,460,5930CG 7.2% 1070.6 / 35.5 338pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (191/117);  new base G (16/9);  total (212/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 387 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCC  >  CP009273/1460526‑1460659
                                                                   |                                                                  
tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                    <  1:2455999/68‑1 (MQ=255)
  gACGGTGACAAGGCGATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCa                                                                  <  1:21329514/68‑1 (MQ=17)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCa                                                                  >  1:17972588/1‑68 (MQ=255)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                  <  1:7343172/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                 <  1:34215975/67‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:11403261/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:23389245/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:4390807/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:24031872/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:8063395/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:31841487/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:31181451/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:15184126/68‑1 (MQ=255)
       tGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                              <  1:28035375/67‑1 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                           >  1:1556865/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                          <  1:32047992/68‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                          <  1:28551705/68‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                          >  1:27139361/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                          >  1:14085047/1‑68 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                         <  1:6371273/68‑1 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:19311782/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:28692401/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:33699895/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:33948501/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:8405882/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:24851543/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:3660355/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:6273340/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:22601142/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:10239241/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGtgat                                                     >  1:20412317/1‑68 (MQ=255)
                   cctCAATAACGAAGGATAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 <  1:32261543/66‑1 (MQ=17)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 <  1:20728353/68‑1 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 >  1:22221432/1‑67 (MQ=255)
                     tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTAACGACGcc                                               >  1:32410150/1‑68 (MQ=14)
                       aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAc                                             <  1:32680603/68‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                            >  1:16476134/1‑52 (MQ=16)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:14230633/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:34798155/1‑68 (MQ=255)
                         tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACgg                                           <  1:1061439/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:6803697/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:12104039/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:26076057/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:1987080/68‑1 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:25529009/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:23007162/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:12227146/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:4724144/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:16866931/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:25034484/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:21017249/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:13303242/1‑68 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:33193450/1‑67 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGAAGACGCCACGGCaa                                        >  1:21889615/1‑67 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     <  1:14747239/68‑1 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     <  1:10874190/68‑1 (MQ=255)
                                  cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                  >  1:9536423/1‑68 (MQ=255)
                                    agagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                >  1:27505138/1‑68 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                              <  1:9147389/67‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                             >  1:12913932/1‑68 (MQ=255)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                         >  1:17245340/1‑67 (MQ=40)
                                             tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacc                       <  1:21277745/68‑1 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                     <  1:13823997/68‑1 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                      >  1:29259513/1‑67 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                      >  1:3101296/1‑67 (MQ=40)
                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGc                                                >  1:11769011/1‑40 (MQ=255)
                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                    >  1:14625046/1‑68 (MQ=40)
                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                    >  1:23728120/1‑68 (MQ=40)
                                                  aaCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                  <  1:18857551/68‑1 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCCAAAAACACGGGCAAAACCCCGGtt                 >  1:19113316/1‑68 (MQ=11)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:26656796/1‑68 (MQ=40)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:33306934/1‑68 (MQ=40)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGAGGACGCCACGGCAAAAAACAGGGGAAAAACCACAGTTg                >  1:18356187/1‑68 (MQ=14)
                                                    cGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:30172258/1‑68 (MQ=25)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCGAAGAACAGGGGAAAAACCACGGTTGa               >  1:33114498/1‑68 (MQ=14)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:10004999/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:31599881/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:7662131/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:17552897/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:16143007/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:25233296/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAAGGGTGAAGACGCCACGGCAAATAACAAGGGAAAAACCACTGTTGa               >  1:30570246/1‑68 (MQ=11)
                                                     ggTGGCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAAAAACAGTGGAAAAACCACTGTTGa               >  1:19634315/1‑68 (MQ=11)
                                                       tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             >  1:11165629/1‑68 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             >  1:19401951/1‑67 (MQ=40)
                                                          caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc           <  1:27901693/67‑1 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                        >  1:3969902/1‑53 (MQ=255)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa         >  1:22406947/1‑68 (MQ=40)
                                                            ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg        >  1:7041590/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:31679284/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:4491324/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:5565468/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:34629006/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:24242391/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:2391343/68‑1 (MQ=40)
                                                             cggcaccCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:18226330/68‑1 (MQ=25)
                                                               gcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt     >  1:29712270/1‑68 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt    >  1:31286231/1‑67 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt    >  1:22774827/1‑67 (MQ=40)
                                                                  cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc  <  1:13592930/68‑1 (MQ=39)
                                                                   |                                                                  
TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCC  >  CP009273/1460526‑1460659

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: