Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,839 T→G 5.3% 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.3% 1531.7 / 6.2 595pseudogene (1191/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (309/208);  new base G (38/2);  total (383/212)
Fisher's exact test for biased strand distribution p-value = 1.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 678 total aligned reads displayed.

ACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACG  >  CP009273/1460772‑1460906
                                                                   |                                                                   
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:19269235/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:23008700/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:27006496/1‑68 (MQ=255)
aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgtt                                                                     >  1:6328588/1‑68 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGCCGGCGACCAGGAGgtcgt                                                                   >  1:13020712/1‑68 (MQ=11)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:6908599/68‑1 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:13865223/68‑1 (MQ=255)
  aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                   <  1:13435479/68‑1 (MQ=255)
   aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgtt                                                                  >  1:29064003/1‑68 (MQ=255)
    gggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTa                                                                 >  1:22596761/1‑68 (MQ=255)
     ggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTAGACGGAGACAAGGCGgttgtaa                                                                 >  1:16177224/1‑65 (MQ=11)
       cACCATGATCCGCACCTTTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                              <  1:31702611/68‑1 (MQ=17)
       cACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCGGTTGTCaaca                                                              >  1:24761142/1‑68 (MQ=255)
        aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                              <  1:2132617/67‑1 (MQ=255)
        aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                              >  1:33327813/1‑67 (MQ=255)
        aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                              <  1:34477361/67‑1 (MQ=255)
         ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                            <  1:11499077/68‑1 (MQ=255)
           aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAAAAACGa                                                          >  1:14182060/1‑68 (MQ=14)
            tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTTAAGAACGaa                                                         >  1:10749912/1‑68 (MQ=255)
              aCCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGTGAAAAACGAAg                                                        >  1:36596754/1‑67 (MQ=14)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAgg                                                       <  1:29183216/67‑1 (MQ=255)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGc                                                      <  1:37310102/68‑1 (MQ=255)
               ccGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGc                                                      >  1:32513923/1‑68 (MQ=255)
                cGTCACCGATCCGGAGTCCATCGGTATCCATGTTGACGGCGACCAGGCGGTTGTGAACAACGAAGGCg                                                     >  1:15779609/1‑68 (MQ=18)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCg                                                     >  1:37032234/1‑68 (MQ=255)
                cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAACAACGAAGGCg                                                     >  1:11388850/1‑68 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAGAACGAAGGCgaga                                                  >  1:10372052/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                  >  1:35990691/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAAGAACGAAGGCgaga                                                  >  1:494013/1‑67 (MQ=255)
                    aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTAAAGAACGAAGGCgaga                                                  >  1:3972692/1‑67 (MQ=14)
                     ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAAAAACGAAGGCGAGAgc                                                >  1:13705583/1‑68 (MQ=255)
                      cGATCCGTAGTCCATCGTTATCCAGGTTGACGGCGAACAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                               <  1:34728556/68‑1 (MQ=18)
                      cGATCCGGAGTCCATCGGTATCCAGGGTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcg                                               <  1:22064692/68‑1 (MQ=255)
                       gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgc                                              <  1:33518261/68‑1 (MQ=255)
                       gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgc                                              <  1:38542211/68‑1 (MQ=255)
                       gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgc                                              <  1:12635270/68‑1 (MQ=255)
                         tCCGGAGTCCATGGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                             >  1:9564008/1‑51 (MQ=255)
                         tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAgcgcaa                                            >  1:27714489/1‑66 (MQ=255)
                         tCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCGTAGTTAAAAACGAAGGCGAGAgcgcaa                                            >  1:19007688/1‑66 (MQ=14)
                          ccGGAGTCCATCGGTATCCAGTTTGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:34996706/68‑1 (MQ=255)
                          ccGGAGTCCATCGGTATCCAGGTTGACGGCGCCCAGGCCGTAGTTAAAAACGAAGGCGAGAGCGCCa                                            >  1:29589645/1‑67 (MQ=17)
                          ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:10045411/68‑1 (MQ=255)
                          ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                           <  1:18645023/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                         <  1:4430461/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                         <  1:32012787/68‑1 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATc                                          >  1:25132187/1‑67 (MQ=255)
                            ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTCGTGAAAAACGAAGGCGAGCGCACGATc                                          >  1:38867733/1‑67 (MQ=11)
                              aGTCCATGGGTATCCAGGTAGACGGCGACAAGGCCCTTGTTAAAAACGAAGGc                                                      >  1:25721087/1‑53 (MQ=11)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                       >  1:30580507/1‑68 (MQ=255)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                       >  1:5795037/1‑68 (MQ=255)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAc                                        <  1:32826625/67‑1 (MQ=255)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTAGTTAAAAACGAAGGCGAGAGCCCAATCAcc                                       >  1:24327010/1‑68 (MQ=14)
                              aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCCTTGTTAAAAACGAAGGCGAGAGCGCCATCAcc                                       >  1:35551308/1‑68 (MQ=17)
                               gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCa                                      >  1:15416540/1‑68 (MQ=255)
                               gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCa                                      >  1:28133960/1‑68 (MQ=255)
                               gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTAAAGAACGAAGGCGAGAGCGCCATCACCa                                      >  1:29794339/1‑68 (MQ=255)
                                  cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                   <  1:26083376/68‑1 (MQ=255)
                                  cATCGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                   <  1:20693553/68‑1 (MQ=255)
                                    tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt                                 >  1:9243859/1‑68 (MQ=255)
                                    tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAGCGCAATCACCAAAggg                                 >  1:14850721/1‑67 (MQ=18)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCGATCACCAAAggcg                                >  1:21047724/1‑65 (MQ=18)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTg                                >  1:38002546/1‑67 (MQ=255)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTg                                >  1:8646796/1‑67 (MQ=255)
                                      ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTg                                >  1:19638394/1‑67 (MQ=255)
                                          tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCAcc                           <  1:17801427/68‑1 (MQ=255)
                                          tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCAcc                           <  1:34245405/68‑1 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                          >  1:23148130/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                          >  1:2554018/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGTGGCACCg                          >  1:31317662/1‑68 (MQ=255)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGCGGCACCg                          >  1:4551418/1‑68 (MQ=18)
                                           ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGCGGCACCg                          >  1:27368386/1‑68 (MQ=18)
                                           ccAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCAATCACCAACGGGGGCACCg                          >  1:16577660/1‑68 (MQ=21)
                                           ccAGGTTGACGGCGACCAGGCCGTGGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGAGGCACCg                          >  1:3304267/1‑68 (MQ=17)
                                             aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:2665806/1‑68 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:11940918/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:36323171/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                        >  1:23537507/1‑67 (MQ=255)
                                              ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGTGGCACCGGc                        >  1:8154848/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:36767420/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      >  1:20316076/1‑67 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                      <  1:1806964/67‑1 (MQ=255)
                                                ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGTGGCACCGGCAc                      >  1:17018304/1‑67 (MQ=255)
                                                  gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCa                   <  1:23956074/68‑1 (MQ=255)
                                                  gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCa                   <  1:14523019/68‑1 (MQ=255)
                                                  gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCa                   <  1:18342392/68‑1 (MQ=255)
                                                  gACGGCGACCAGGCGGGGGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCa                   <  1:37499365/68‑1 (MQ=255)
                                                    cGGGGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:30304422/68‑1 (MQ=18)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:761887/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:23749302/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:16663712/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:6091115/68‑1 (MQ=255)
                                                    cGGCGACCAGGCGGTTGTTAACAACCAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:7174506/68‑1 (MQ=255)
                                                    cGGAGACAAGGCGGTTGTTATTAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                 <  1:21554147/68‑1 (MQ=17)
                                                     ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                >  1:5244617/1‑68 (MQ=255)
                                                          ccAGGCGTTTGTTAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACg           <  1:1101646/68‑1 (MQ=18)
                                                            aGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt         <  1:30989036/68‑1 (MQ=255)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg        <  1:28808479/66‑1 (MQ=255)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg        <  1:33708033/66‑1 (MQ=255)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGtaaa      >  1:19891282/1‑65 (MQ=21)
                                                               cGGTTGTTAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg        <  1:17417943/66‑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 < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: