Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,653 G→T 11.3% pseudogene (1005/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,6530GT11.3% 843.8 / 7.6 317pseudogene (1005/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (129/126);  new base T (5/37);  total (136/181)
Fisher's exact test for biased strand distribution p-value = 1.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

CGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460587‑1460720
                                                                  |                                                                   
cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:29737631/68‑1 (MQ=40)
         gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                           >  1:18706840/1‑68 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACGACCGTTGACGGCAAAGATTCCACCggg                                                         >  1:21231026/1‑67 (MQ=11)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:10744639/68‑1 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:389279/68‑1 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:38528112/68‑1 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:36953722/68‑1 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:30397761/68‑1 (MQ=35)
           ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:32622144/68‑1 (MQ=35)
           ttAATGGTGATGAAGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:42111878/67‑1 (MQ=255)
            tAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         >  1:28518364/1‑67 (MQ=35)
               tGGTGATGACGCCACGGCAAACAACAACGGCAAATCCACCGTTGACGGCAAGGAt                                                                  >  1:39660070/1‑55 (MQ=9)
               tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                     >  1:43128355/1‑68 (MQ=31)
                    atgaCGCCACGGCAAAAAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACGGAAATc                                                >  1:33772862/1‑68 (MQ=11)
                     tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                >  1:5731588/1‑67 (MQ=31)
                     tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               >  1:15385164/1‑68 (MQ=33)
                     tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCa                                               >  1:8463470/1‑68 (MQ=255)
                       aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGATTGAAAAGATTCCACCGGTACCCAAATcgcc                                             >  1:28102317/1‑65 (MQ=11)
                       aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                             >  1:5889700/1‑68 (MQ=37)
                        cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:13149228/68‑1 (MQ=39)
                        cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:3791078/68‑1 (MQ=39)
                        cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:37166976/68‑1 (MQ=39)
                          ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGATATCCACGGt                                          >  1:21588970/1‑68 (MQ=255)
                           cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGt                                          <  1:16864270/67‑1 (MQ=255)
                              ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:42054438/62‑1 (MQ=39)
                              ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:33497330/68‑1 (MQ=255)
                              cgCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      <  1:17351617/67‑1 (MQ=255)
                               gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:11939804/1‑67 (MQ=255)
                                   acaacaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGCACGGAAATCAACGGTAATAACGGa                                 <  1:6258919/68‑1 (MQ=11)
                                   acaacaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                 <  1:43623743/68‑1 (MQ=255)
                                      acaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:31771101/68‑1 (MQ=255)
                                      acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:23014658/68‑1 (MQ=255)
                                      acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:1926731/68‑1 (MQ=255)
                                      acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:6975334/68‑1 (MQ=255)
                                      acaatGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:26036722/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAATGATTCCACAGGTACGGAAATAAAAGAAAATAACGGaaa                               >  1:2249855/1‑67 (MQ=14)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACTGTAATAACGGaaaa                              >  1:35728533/1‑68 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:43322186/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:8849009/1‑68 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:11222558/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:10121942/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:5454279/1‑68 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:24378012/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:23750766/1‑68 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:39985644/68‑1 (MQ=255)
                                      acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:19871444/1‑68 (MQ=255)
                                        tatggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:36649682/65‑1 (MQ=255)
                                        aatggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:41727617/65‑1 (MQ=255)
                                        aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCTACGGTAATAACGGAAAAGt                            >  1:36128548/1‑68 (MQ=255)
                                        aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:4967325/68‑1 (MQ=255)
                                        aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:35120002/68‑1 (MQ=255)
                                        aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGAAATAACGGAAAAGt                            >  1:7970198/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:29266434/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:38395438/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:23316772/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATGAACGGGAATAACGGAAAAGTGAt                         >  1:17419639/1‑68 (MQ=14)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGAAATAACGGAAAAGTTAt                         >  1:30749570/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAATGTTAt                         >  1:10865602/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:27109932/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:20922772/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:40453293/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:1806517/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGAAAAAGTTAt                         <  1:19325497/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGGAATAACGGAAAAGTTAt                         >  1:2212192/1‑68 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:12020062/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:37369672/68‑1 (MQ=255)
                                           ggCAAAACCACCGTTGACGGAAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:19642728/68‑1 (MQ=255)
                                              aaaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCa                      <  1:36439393/68‑1 (MQ=255)
                                                aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:21822952/68‑1 (MQ=255)
                                                aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:32172689/68‑1 (MQ=255)
                                                 accaccGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                   <  1:15926097/68‑1 (MQ=255)
                                                 accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGTAATAACGGAAAAGTGATTCAgg                    >  1:29857036/1‑67 (MQ=255)
                                                 accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:39004937/1‑67 (MQ=255)
                                                 accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:12574971/1‑67 (MQ=255)
                                                 accaccGTTGACGGCAAGGATTCCACCGGGACGGAAATCAACGGTAATAACGGAAAAGTGATTCggg                    >  1:24181680/1‑67 (MQ=14)
                                                   caccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg                 <  1:26255496/68‑1 (MQ=255)
                                                   caccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg                 <  1:7391852/68‑1 (MQ=255)
                                                     ccGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:2821722/68‑1 (MQ=255)
                                                     ccGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:31195887/68‑1 (MQ=255)
                                                     ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:33962780/68‑1 (MQ=255)
                                                       gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             <  1:42467155/68‑1 (MQ=255)
                                                       gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:34397296/1‑68 (MQ=255)
                                                        ttGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:40949744/1‑67 (MQ=255)
                                                         tGACGGGAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:33101042/67‑1 (MQ=255)
                                                         tGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:31371562/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:41542624/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:41243380/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:24653973/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:11432044/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:1473186/67‑1 (MQ=255)
                                                         tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:6628932/67‑1 (MQ=255)
                                                         tGACGGCAACGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:20651753/67‑1 (MQ=14)
                                                          gACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt          <  1:22392751/68‑1 (MQ=255)
                                                            cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGTGATCTgg        <  1:9042269/68‑1 (MQ=255)
                                                            cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:14074189/1‑68 (MQ=255)
                                                            cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        >  1:16352781/1‑68 (MQ=255)
                                                            cGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        <  1:33639978/68‑1 (MQ=255)
                                                             ggaaaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa       <  1:25047099/65‑1 (MQ=255)
                                                              gCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAt      >  1:27466439/1‑68 (MQ=255)
                                                              gCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAt      >  1:28966557/1‑68 (MQ=255)
                                                                  |                                                                   
CGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460587‑1460720

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: