Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 6.4% 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 6.4% 1301.5 / 35.4 412pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (156/221);  new base G (13/11);  total (180/232)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 458 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCC  >  CP009273/1460526‑1460659
                                                                   |                                                                  
tCGACGGTGACAAGGCCATCGTCAATAACGGAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                    <  1:22312119/68‑1 (MQ=255)
tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                    <  1:14483962/68‑1 (MQ=255)
  gACGGTGACATGGCCATCGTCAATAACGAAGTCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                  <  1:12104561/68‑1 (MQ=255)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAAAGGAGGCACCGGCACTCa                                                                  >  1:5388311/1‑68 (MQ=11)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                  <  1:160136/68‑1 (MQ=255)
    ctgTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCCTGACCAACGGTGGCACCGGCACCCAGa                                                                <  1:18549802/66‑1 (MQ=18)
    cGGTGACAAGGCCATGGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                <  1:35306321/68‑1 (MQ=18)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                <  1:16847578/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCTGCACCCAGa                                                                <  1:25548818/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCTCCCAGa                                                                <  1:32984688/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                 <  1:19492307/67‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:13607843/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:38177272/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:12585296/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:9039484/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:42537094/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:32452785/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:29809077/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:19108012/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:28765162/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGACCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:18023331/68‑1 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCGCAGAt                                                               >  1:3204390/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGAt                                                               >  1:839813/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                               >  1:9391980/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                               >  1:22548614/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                               >  1:18396648/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGAGGCACCGGCACACAGAt                                                               >  1:23269855/1‑68 (MQ=14)
       tGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                              <  1:27789437/67‑1 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGTGGCACCGGCACTCAGATCaac                                                           >  1:36149320/1‑67 (MQ=11)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                           <  1:29149917/67‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                           <  1:33245663/67‑1 (MQ=255)
            aGGCGATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                        <  1:16969534/68‑1 (MQ=11)
            aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                        <  1:42853227/68‑1 (MQ=255)
            aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGAGGCACCGGCACGCAGATTAAAGGt                                                        >  1:11150376/1‑68 (MQ=14)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:7824018/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACGCAGATGAACGGta                                                       >  1:11386930/1‑67 (MQ=11)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGTGGCACCGGCACCCAGATTAAAGGtg                                                       >  1:30834483/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:42783194/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:24765204/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:37474887/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:35638284/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACTCAGATTAAAGGtaac                                                     >  1:13273015/1‑65 (MQ=11)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 <  1:28979224/68‑1 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 <  1:6770834/68‑1 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 >  1:386077/1‑67 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 >  1:2413498/1‑67 (MQ=255)
                     tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATGAACGGTGAGGACGc                                                >  1:29202097/1‑67 (MQ=11)
                       aaTAACGAAGGCTAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAc                                             <  1:20495718/68‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGAt                                                               >  1:28132816/1‑49 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:35986901/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:14839191/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:33611728/1‑68 (MQ=255)
                         tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACgg                                           <  1:27397976/68‑1 (MQ=255)
                          caCGAAGGCGAGAGCACGATCACCTACGTTGGCACCGGGACCCATATTAATGGTGATGACGCCACTGc                                          <  1:10517623/67‑1 (MQ=14)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          >  1:31390345/1‑68 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:21710814/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:20881078/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          >  1:38996582/1‑68 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:37311732/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          >  1:33803790/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:23504725/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:26256261/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACAATCACCAACGGTGGCACCGGCACTCAGATaaa                                                            >  1:41572173/1‑49 (MQ=11)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTTATGGTGATGACGCCACGGCaa                                        >  1:5772449/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:15244976/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        <  1:6096040/68‑1 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                       <  1:3408987/68‑1 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        >  1:43092713/1‑67 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     >  1:36591741/1‑68 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     >  1:39964893/1‑67 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     >  1:7292006/1‑67 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     >  1:17242070/1‑67 (MQ=255)
                                       gCACGATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                              <  1:36793510/67‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                             >  1:12952724/1‑68 (MQ=255)
                                           cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                         <  1:26674631/68‑1 (MQ=40)
                                           cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                          >  1:40358044/1‑67 (MQ=40)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                         >  1:25849020/1‑67 (MQ=40)
                                             tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAcc                       >  1:12336488‑M2/1‑66 (MQ=255)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAAAAACACCGGAAAAAcca                      <  1:34604380‑M2/67‑4 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCGAATAACAGGGGAAAAAcc                       <  1:8874125‑M2/64‑3 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCCAAAAACAGGGGAAAAAccacgg                   <  1:39591925‑M2/68‑7 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccg                   >  1:21419262‑M2/1‑62 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccg                   <  1:9231028‑M2/68‑7 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAAAGGTGAAGACGCCACGGCGAAAAACAGGGGAAAAaccacgg                   >  1:2114903/1‑66 (MQ=11)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgtt                 <  1:11822541‑M2/68‑9 (MQ=255)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgtt                 <  1:30545798‑M2/68‑9 (MQ=255)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCTAAAAACACCGGAAAAAccacagttg                <  1:43536262‑M2/68‑10 (MQ=255)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgttg                <  1:17664857‑M2/68‑10 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgttg                <  1:25185834‑M2/67‑10 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACCGCAAAAccaccgttg                <  1:4973881‑M2/67‑10 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgttg                <  1:34862739‑M2/67‑10 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATAAACGGTGATGACGCCACGGCAAACAACAACGGCAAAAccaccgt                  >  1:1928293‑M2/1‑58 (MQ=255)
                                                      gTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAAAAACAAAGGAAAAACCACAGTTGag              >  1:16513426/1‑67 (MQ=11)
                                                       tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             >  1:2376503/1‑68 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             <  1:41369000/67‑1 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAATAACAGGGGAAAAACCCCGGTTGATGGCaa         >  1:40759901/1‑68 (MQ=14)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa         >  1:10147900/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCGAAAAACAGAGGAAAAACCACGGTTGATGGAAAgg       >  1:20006763/1‑68 (MQ=14)
                                                               gcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt     >  1:24334786/1‑68 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt    <  1:42111103/67‑1 (MQ=40)
                                                                   |                                                                  
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: