Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,839 T→G 5.8% pseudogene (1191/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,8390TG 5.8% 2103.5 / 12.1 776pseudogene (1191/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (406/273);  new base G (49/6);  total (494/282)
Fisher's exact test for biased strand distribution p-value = 4.83e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 895 total aligned reads displayed.

ACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG  >  CP009273/1460772‑1460906
                                                                   |                                                                   
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:24831239/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:48984377/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:47993151/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:50291399/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:22157615/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:26778088/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:26355104/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:38307577/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:43348492/1‑68 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCggtggt                                                                   >  1:11304328/1‑68 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:36841136/68‑1 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:8867780/68‑1 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:6021951/68‑1 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtagt                                                                   >  1:18671362/1‑68 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCGgtagt                                                                   >  1:24184560/1‑68 (MQ=255)
       cACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                              >  1:41027057/1‑68 (MQ=255)
        aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGAGACCAGGCGGTTGTTaaca                                                              >  1:38001068/1‑67 (MQ=17)
        aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGaaca                                                              >  1:16823932/1‑67 (MQ=255)
        aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGaaca                                                              >  1:845345/1‑67 (MQ=255)
         ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                            <  1:44324344/68‑1 (MQ=255)
         ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                            <  1:7041316/68‑1 (MQ=255)
           aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGa                                                          <  1:51143928/68‑1 (MQ=255)
           aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGa                                                          >  1:15363425/1‑68 (MQ=255)
            tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGaa                                                         >  1:23502324/1‑68 (MQ=255)
              aCCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTAAAAAACGAAg                                                        >  1:17272256/1‑67 (MQ=14)
              aCCGTCAACGATAAGGAGTCAATCGGGATCCAGGTTTACGGCGACCAGGCGGTTGTTAACAACGAAg                                                        >  1:2955520/1‑67 (MQ=18)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:45448856/1‑53 (MQ=255)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAgg                                                       <  1:28311026/67‑1 (MQ=255)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAgg                                                       <  1:29833956/67‑1 (MQ=255)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGGGACCAGGCGGTTGTTaaaaa                                                             >  1:40474889/1‑60 (MQ=17)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCg                                                     >  1:20934275/1‑68 (MQ=255)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAAGAACGAAGGCg                                                     >  1:49883865/1‑68 (MQ=255)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGTGAACAACGAAGGCg                                                     >  1:27757706/1‑68 (MQ=255)
                   cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                  <  1:47150275/68‑1 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                  >  1:34916035/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                  >  1:18444612/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                  >  1:30804179/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAACAACGAAGGCgaga                                                  >  1:14815476/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTAGTTAAAAACGAAGGCgaga                                                  >  1:19859346/1‑67 (MQ=14)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCCTGGTAAACAACGAAGGCgaga                                                  >  1:15991511/1‑67 (MQ=11)
                    aCCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCGTTGTGAAGAACGAAGGCgagc                                                  >  1:31947106/1‑66 (MQ=11)
                    aCCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCGTTGTGAAGAACGAAGGCgaga                                                  >  1:39582590/1‑67 (MQ=11)
                      cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                               <  1:22628491/68‑1 (MQ=255)
                      cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                               <  1:34447729/68‑1 (MQ=255)
                      cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgc                                                >  1:42051582/1‑67 (MQ=255)
                       gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgc                                              <  1:45927129/68‑1 (MQ=255)
                       gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgc                                              <  1:46881459/68‑1 (MQ=255)
                         tCCGGAGTCCATCGGTATCCAGGTTGACGGCGCCCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCa                                            >  1:34005546/1‑68 (MQ=255)
                          ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:5431233/68‑1 (MQ=255)
                          ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:11622687/68‑1 (MQ=255)
                          ccGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:36419807/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATc                                          >  1:41800550/1‑67 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                         <  1:23050219/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                         <  1:46509536/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                         <  1:28917163/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAAGAACGAAGGCGAGAGCCCAATc                                          >  1:19100251/1‑67 (MQ=14)
                              aGTCCCTCGGTATCCAGTTTGACGGCGACCAGGCTGGTGTTAACAACTAAGGCGAGAGCGCCATCAc                                        <  1:32293976/67‑1 (MQ=18)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                       >  1:49460908/1‑68 (MQ=255)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                       >  1:12567974/1‑68 (MQ=255)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAc                                        <  1:42802606/67‑1 (MQ=255)
                               gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCa                                      >  1:12980223/1‑68 (MQ=255)
                               gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAGCGCCATCACCa                                      >  1:38675850/1‑68 (MQ=255)
                                 ccATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCaa                                     <  1:12081578/67‑1 (MQ=255)
                                  cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                   <  1:48580541/68‑1 (MQ=255)
                                  cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                   <  1:29187447/68‑1 (MQ=255)
                                   aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACgg                                  <  1:39824305/68‑1 (MQ=255)
                                     cGGTATCCAGCTGGACGGGGACCAGGCGGTGGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTg                                <  1:803128/68‑1 (MQ=18)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTg                                >  1:6290089/1‑67 (MQ=255)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAggag                                >  1:35264962/1‑65 (MQ=255)
                                      ggTATCCAGGTTGACGGCGACCAGGCCGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACggag                                >  1:33575031/1‑65 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                          >  1:44685406/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCg                          >  1:43636144/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGTGGCACCg                          >  1:10004912/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTAGTTAAGAACGAAGGCGAGAGCGCAATCACCAAAGGCGGCACCg                          >  1:10326623/1‑68 (MQ=17)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCGATCACCAAAGGGGGCACCGGc                        >  1:1365209/1‑67 (MQ=21)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:17708334/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:6709351/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGc                        >  1:20415028/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACTATCACCAATGGCGGCa                             >  1:33187208/1‑62 (MQ=14)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      <  1:47549555/67‑1 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:14320378/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:25303703/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:50556770/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:18212089/1‑67 (MQ=255)
                                                  gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGTGGCACCGGCACGCa                   >  1:21845592/1‑68 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAg                  <  1:30285646/67‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:9583967/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:39211842/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:9276701/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:10862608/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:1944261/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:21367590/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:2737909/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:32770656/68‑1 (MQ=255)
                                                        gACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa             <  1:3728966/68‑1 (MQ=255)
                                                           aaGGCGTTTGTTAATAACGAAGGCGAGAGCGGCATCACCAACGGTGGCACCGGCACGCAGATTAACgg          <  1:37665283/67‑1 (MQ=17)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgat      >  1:3592132/1‑68 (MQ=255)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACTCCGATAAACGGtgaa      >  1:7626180/1‑67 (MQ=18)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACAAACGGTGGCACCGGAACGCAGATTAACGGtgat      >  1:25959419/1‑68 (MQ=255)
                                                                   tgttAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACg  <  1:31536768/68‑1 (MQ=255)
                                                                   |                                                                   
ACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG  >  CP009273/1460772‑1460906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: