Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 8.8% 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 8.8% 1628.1 / 65.7 531pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (267/205);  new base G (25/17);  total (309/222)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 596 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCA  >  CP009273/1460526‑1460660
                                                                   |                                                                   
tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                     <  1:46198339/68‑1 (MQ=255)
  gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa                                                                   <  1:35085308/68‑1 (MQ=255)
   aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAAAGGAGGCACCGGCACACAg                                                                  >  1:19395204/1‑68 (MQ=14)
    cGGTGACAAGGCCATCGTCAATAACGTAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:5903239/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:895600/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:36627140/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:41880154/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:43343477/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:2538997/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:24271040/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:21678431/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:2091282/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:20051578/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:19723681/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:18423740/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:2794756/68‑1 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGAt                                                                >  1:48331767/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                                >  1:8645746/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCAACCAGAt                                                                >  1:15780510/1‑68 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAACGGTGGCACCGGCACGCAGAt                                                                >  1:27120114/1‑68 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                             >  1:6955530/1‑67 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGTGGCACCGGCACTCAGATaaac                                                            >  1:7276378/1‑67 (MQ=11)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            <  1:51084916/67‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                           >  1:50654873/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                           >  1:33306767/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                           >  1:11627172/1‑68 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:50101581/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:45197224/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:7607417/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:15366236/68‑1 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:46781155/1‑68 (MQ=255)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       >  1:49023962/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACGATCACCAAGGGTGGCACCGGCACTCAGATGAAGGGtgac                                                      >  1:39093991/1‑67 (MQ=11)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACgggtg                                                       >  1:18919708/1‑64 (MQ=14)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:40481415/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:32217028/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:37019059/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:33619623/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:14574696/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        <  1:12597252/66‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        <  1:24030745/66‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGCCACGCAGATAAACGGtgaa                                                      >  1:42789857/1‑67 (MQ=14)
                 aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatg                                                     >  1:11295069/1‑67 (MQ=255)
                  tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                   <  1:22698871/68‑1 (MQ=255)
                  tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                   <  1:38468488/68‑1 (MQ=255)
                   cGTCAATAACGAAGGCGATAGCACCATCACCAACGGGGTCACCGGCACCCAGATTAATGGTGATGACg                                                  <  1:39678753/68‑1 (MQ=18)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                  <  1:45080531/68‑1 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGACg                                                  <  1:30026314/68‑1 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                  >  1:3736226/1‑67 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                  >  1:10866670/1‑67 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:38080173/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             <  1:11982321/68‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACGCAGATGAAAGGTGAAGACGCCACg                                             >  1:8088513/1‑68 (MQ=11)
                          aaCGAAGGCGAGAGCACCATCACCCACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:954146/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:36036526/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:17556693/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           >  1:41168591/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:25991717/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:49430022/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:33912578/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:13576312/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          <  1:14293349/68‑1 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACCCAGATAAACGGTGACGCCGCCACGGcg                                          >  1:47518971/1‑67 (MQ=11)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGcta                                         <  1:29069358/68‑3 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         <  1:14907672/68‑1 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         <  1:23122732/68‑1 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:24914494/1‑68 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                      <  1:28312198/68‑1 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                      >  1:1812393/1‑67 (MQ=255)
                                  cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                     <  1:6594144/66‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                               <  1:16883736/67‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                <  1:40833052/66‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                <  1:19976502/66‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                              >  1:47036751/1‑68 (MQ=255)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                          >  1:17188245/1‑67 (MQ=40)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                          >  1:4453039/1‑67 (MQ=40)
                                              cACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacc                        >  1:13053455/1‑67 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                       >  1:38779276/1‑67 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                       >  1:34117259/1‑67 (MQ=40)
                                                ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                     >  1:32703133/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                    >  1:39397433/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:101352/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:45501510/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:43943655/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:31136549/1‑68 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCGAATAACAGAGGAAAAACCACGGTTGa                >  1:5488030/1‑68 (MQ=14)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAACCCACCGTTg                 >  1:34552115/1‑67 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:26693707/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:42395244/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:1598315/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:31871427/1‑67 (MQ=40)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc               >  1:13889621/1‑68 (MQ=40)
                                                       tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              >  1:51294246/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:3571514/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:27538993/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:47744033/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:30622624/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGACAAGGCAAACAACAACGGCAAAACCACCGTGGACGGCAAgg        >  1:4134232/1‑68 (MQ=14)
                                                              ggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa       <  1:2219974/68‑1 (MQ=40)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc   >  1:29496445/1‑67 (MQ=37)
                                                                   |                                                                   
TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCA  >  CP009273/1460526‑1460660

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: