Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 7.6% 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.6% 1758.3 / 59.7 563pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (306/205);  new base G (22/16);  total (342/221)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 656 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCA  >  CP009273/1460526‑1460660
                                                                   |                                                                   
    cGGTGACAAGGCCTTTGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:1980333/68‑1 (MQ=14)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:23781715/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                  <  1:23960387/67‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:21487748/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:20611923/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:57443039/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:31247746/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:55890419/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:2943418/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:51123829/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                 <  1:47857033/68‑1 (MQ=255)
    cGGGGACAAGGCCTTCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                 <  1:35651507/68‑1 (MQ=17)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                                >  1:27245882/1‑68 (MQ=255)
        gtcAAGGCCATCGTCATTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaa                                                             <  1:50983965/66‑1 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            >  1:51373667/1‑68 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            >  1:14653430/1‑68 (MQ=255)
         aCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            >  1:50147298/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTaagg                                                           >  1:44114185/1‑66 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                            <  1:3234192/67‑1 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                           >  1:46408624/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                           >  1:24764087/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACGCAGATCaacg                                                           >  1:12018302/1‑66 (MQ=11)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGTGGCACCGGCACGCAGATTaacg                                                           >  1:8066967/1‑66 (MQ=14)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:13332814/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:52063077/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                          <  1:34287157/68‑1 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGtg                                                        >  1:3787591/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:50481957/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:1905935/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:28455720/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        >  1:22043141/1‑68 (MQ=255)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       >  1:49315619/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                      >  1:34772714/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:30173908/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                       <  1:45429473/67‑1 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                        <  1:1652461/66‑1 (MQ=255)
                  tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                   <  1:6858696/68‑1 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                   >  1:22610027/1‑67 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                  <  1:38951384/68‑1 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                  >  1:56826392/1‑67 (MQ=255)
                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                <  1:3728072/67‑1 (MQ=255)
                      cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGcc                                                <  1:45006492/67‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:5458257/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:41637378/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             >  1:42503781/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                             <  1:53229887/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:1766581/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                           <  1:23426269/68‑1 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCa                                          >  1:2589874/1‑68 (MQ=25)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:40293430/1‑68 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                          >  1:9297218/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:1002195/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         <  1:18102097/68‑1 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:15253847/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:28973333/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                         >  1:37675619/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACAAACTGTTGCACCGGAACCAAGAATAAAGGAGATTACGCCACGGCaa                                         >  1:26389500/1‑68 (MQ=9)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                        <  1:52482567/68‑1 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                      <  1:19540986/68‑1 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                      >  1:32990098/1‑67 (MQ=255)
                                 gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaac                                    >  1:51574807/1‑68 (MQ=255)
                                 gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaac                                    >  1:41061643/1‑68 (MQ=255)
                                  cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                   >  1:13902363/1‑68 (MQ=255)
                                       tcACCATCACAAACGGCTGCACCGTCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                               <  1:4281775/66‑1 (MQ=14)
                                       gCACCATCACCAACGGTGGCACCGGCACTCAGATAAAAGGTGAGGACGCCCAGGCAATAAAAAACGGc                              >  1:43561566/1‑68 (MQ=2)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                <  1:53977916/66‑1 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                              >  1:33344507/1‑68 (MQ=255)
                                         aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                            >  1:49676394/1‑68 (MQ=255)
                                         aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                            >  1:18413634/1‑68 (MQ=255)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                          >  1:39899429/1‑67 (MQ=40)
                                               aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                       >  1:6358088/1‑67 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                    >  1:30887764/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                    >  1:16173461/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:35195674/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:25358748/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:36162410/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                  >  1:7781768/1‑68 (MQ=40)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGAAGACGCCACCGCGaaaaa                                     >  1:8593378/1‑49 (MQ=14)
                                                   aCGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCggt                  >  1:5902334/1‑66 (MQ=255)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:7077119/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:22809539/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:36841146/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:48507777/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                 >  1:28085515/1‑67 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAAGGGTGAGGACGCCACGGCGAAAAACAGAGGAAAAACCACTGTTGa                >  1:39413981/1‑68 (MQ=11)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGAAGCCAAGGAAAACAAAAACGGCAAAACCACCGTTGAc               >  1:10558975/1‑68 (MQ=14)
                                                       tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                >  1:20165417/1‑50 (MQ=255)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              <  1:11615735/67‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg              <  1:33691031/67‑1 (MQ=40)
                                                          caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc            <  1:41964535/67‑1 (MQ=40)
                                                          caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAAAGGCAAAACCACCGTTGACGGCa           >  1:31837020/1‑68 (MQ=255)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                       >  1:40492851/1‑55 (MQ=255)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa          >  1:12543979/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:46642181/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:55097145/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg        <  1:16838671/68‑1 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt     <  1:47364962/67‑1 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAAAAACATCGGAAAAACCACGGTTGACGGAAAGGAtt     >  1:31581553/1‑67 (MQ=14)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc   >  1:17205775/1‑67 (MQ=37)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa  <  1:12882684/68‑1 (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: