Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,653 G→T 13.5% 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,6530GT13.5% 732.0 / 10.2 298pseudogene (1005/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (133/99);  new base T (3/42);  total (137/161)
Fisher's exact test for biased strand distribution p-value = 4.77e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

CCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460586‑1460720
                                                                   |                                                                   
ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                     >  1:33966982/1‑68 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:1064691/68‑1 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:34866303/68‑1 (MQ=40)
 cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                    <  1:6807081/68‑1 (MQ=40)
  ggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                   <  1:21170936/68‑1 (MQ=40)
     accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                 <  1:19465986/67‑1 (MQ=40)
       ccAGATTAATGGTGATGACGCCACGGCCAACAACAGGGGCCAACCCACAGTTGACGGCAAGGATTcc                                                               >  1:2341826/1‑67 (MQ=11)
       ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                               >  1:5687515/1‑67 (MQ=37)
          gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                           >  1:12575940/1‑68 (MQ=35)
          gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                           >  1:5533988/1‑68 (MQ=35)
          gATTAATGGTGATGACGCCACGGCAAAAAACAACGGAAAAACCACCGTTGAAGGAAAGGATTCCACCg                                                           >  1:36592050/1‑68 (MQ=11)
            ttACTGGTGATGACGCCACTGCAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                          <  1:38368233/67‑1 (MQ=11)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:32119421/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         >  1:17013330/1‑68 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:14224835/68‑1 (MQ=35)
            ttAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                         <  1:8736093/68‑1 (MQ=35)
               aTGGTGATGACGCCACGGCAAACAACAAGGGAAAAAACAAAGTTGACGGAAATGAATCCACCGGTAc                                                       >  1:34604902/1‑67 (MQ=11)
               aTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                              >  1:4507513/1‑60 (MQ=33)
               aTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                      >  1:17176809/1‑68 (MQ=33)
               aTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                      >  1:31181661/1‑68 (MQ=33)
                      tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                               <  1:3739301/68‑1 (MQ=33)
                        aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                             >  1:32480831/1‑68 (MQ=37)
                        aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                             >  1:33539930/1‑68 (MQ=37)
                        aCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                             >  1:26472207/1‑68 (MQ=37)
                         tgCCACTGAAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            <  1:27422016/67‑1 (MQ=11)
                         cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATGAACg                                            >  1:6492899/1‑68 (MQ=18)
                         cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACg                                            >  1:20568881/1‑68 (MQ=39)
                         cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAACg                                            >  1:6973842/1‑68 (MQ=18)
                          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACgg                                           <  1:34237610/68‑1 (MQ=40)
                            cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGta                                         <  1:19788879/68‑1 (MQ=255)
                                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGgaata                                      >  1:1888522/1‑67 (MQ=255)
                                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                      >  1:11224078/1‑67 (MQ=255)
                                  aaaCAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:1843423/68‑1 (MQ=25)
                                  aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTTCGGAAATCGACGGGAATAACg                                   >  1:8337173/1‑68 (MQ=14)
                                  aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                   <  1:35851970/68‑1 (MQ=255)
                                    acaacaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                 <  1:16488114/68‑1 (MQ=255)
                                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                 >  1:14741016/1‑68 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:8559777/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:28224350/68‑1 (MQ=255)
                                       acaatGGCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:15709739/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCCCCGGGAATAACGGcaat                              >  1:12193647/1‑67 (MQ=11)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:10889589/1‑68 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:15509962/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:32696296/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:26941345/68‑1 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              >  1:9047046/1‑68 (MQ=255)
                                       acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAACCGGTAATAACGGaaaa                              >  1:29611308/1‑68 (MQ=255)
                                       acaacGGCAAAACCACCGGTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:31838885/68‑1 (MQ=255)
                                       acaaaGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                              <  1:17709583/68‑1 (MQ=255)
                                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAg                             <  1:36162587/68‑1 (MQ=255)
                                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            >  1:13006375/1‑68 (MQ=255)
                                         aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                            <  1:10394027/68‑1 (MQ=255)
                                          atggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                           <  1:4037775/66‑1 (MQ=255)
                                          acGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGGAATAACGGAAAAGtt                           >  1:25674155/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:24654604/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCTACGGTAATAACGGAAAAGTTAt                         >  1:27849801/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:20912342/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         >  1:33745527/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         >  1:25439448/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:1491506/68‑1 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAAAGGTAATAACGGAAAAGTTAt                         >  1:4961363/1‑68 (MQ=255)
                                            ggCAAAACCACCGTTGACGGCAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                         <  1:10252645/68‑1 (MQ=255)
                                              cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                       >  1:32196975/1‑68 (MQ=255)
                                               aaaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                       <  1:3819059/67‑1 (MQ=255)
                                                 aaCCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:6225240/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:25017237/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:10025869/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:38840250/68‑1 (MQ=255)
                                                 aaCCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:32976117/68‑1 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAGGGTAATAACGGAAAAGTTATTCAgg                    >  1:22636977/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:6733676/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:20690586/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    >  1:21650571/1‑67 (MQ=255)
                                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGGAATAACGGAAAAGTTATTCAgg                    >  1:26329581/1‑67 (MQ=255)
                                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                  <  1:30383954/68‑1 (MQ=255)
                                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                  >  1:17390768/1‑68 (MQ=255)
                                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                  >  1:35588321/1‑68 (MQ=255)
                                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                  >  1:27051608/1‑68 (MQ=255)
                                                    caccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGAAATAACGGAAAAGTTATTAAGCACg                 >  1:27934180/1‑68 (MQ=17)
                                                     accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTATAACGGAAAAGTTATTCAGGACGGc               <  1:36278292/68‑1 (MQ=255)
                                                      ccGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:7871046/68‑1 (MQ=255)
                                                      ccGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:17903271/68‑1 (MQ=255)
                                                      ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:1192691/68‑1 (MQ=255)
                                                      ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:29121315/68‑1 (MQ=255)
                                                      ccGTTGACGGCAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc               <  1:23034031/68‑1 (MQ=255)
                                                        gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             <  1:22275026/68‑1 (MQ=255)
                                                         ttGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                    <  1:13727029/60‑1 (MQ=255)
                                                          tGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:13356208/67‑1 (MQ=255)
                                                          tGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:15157804/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:19124368/1‑66 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:5895094/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:7204477/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:1808571/67‑1 (MQ=255)
                                                          tGACGGCAAGGATTCCACCGGTACGGAAATAAACGGGAATAACGGAAAAGTTATTCAGGACGGCGa             >  1:20410170/1‑66 (MQ=17)
                                                          tGACGGCAACGATTCCACCGGTACGGAAAGCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGAt            <  1:14392040/67‑1 (MQ=255)
                                                           gACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt          <  1:23665898/68‑1 (MQ=255)
                                                             cGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        <  1:11704327/68‑1 (MQ=255)
                                                             cGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg        <  1:9266306/68‑1 (MQ=255)
                                                               gCAAGGATTCCACCGGTACGGAAATCCACGGTAATAACGGAAAAGTAAATAAGGACGGCGAACTGGAt      >  1:17548039/1‑68 (MQ=14)
                                                               gCAAGGATTACACCGGTACGGAAATCAACGGAAATAACGGCAATGTTATTCAGGACGGAGATATTGAt      >  1:4683019/1‑68 (MQ=14)
                                                                   |                                                                   
CCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCA  >  CP009273/1460586‑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: