Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,593 C→G 8.7% 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 8.7% 1088.8 / 46.9 360pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (189/132);  new base G (15/13);  total (215/145)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 421 total aligned reads displayed.

TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCC  >  CP009273/1460526‑1460659
                                                                   |                                                                  
tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                    <  1:10377594/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCACTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:24546793/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAg                                                                 <  1:39426310/67‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                <  1:4469570/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:16753162/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:28088214/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:19034051/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:20227920/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:22122042/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:24979290/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:2616125/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:26719896/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:29513211/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:30527973/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:32041367/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:33023619/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:38497783/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:13990764/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:11371896/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:15309040/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:10627364/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:14596785/68‑1 (MQ=255)
    cGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                <  1:23298297/68‑1 (MQ=255)
     ggTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                               >  1:29813104/1‑68 (MQ=255)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACTCAGATaaacg                                                          >  1:37061106/1‑66 (MQ=14)
          cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCa                                                                  >  1:18753093/1‑60 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                         <  1:8579473/68‑1 (MQ=255)
           aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATgg                                                         <  1:33438749/68‑1 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:20505503/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:23718773/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:1003272/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                       >  1:28837970/1‑68 (MQ=255)
             ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGtg                                                       >  1:8906120/1‑68 (MQ=255)
              gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                      >  1:25853891/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATGAACGGtgaa                                                     >  1:37510129/1‑67 (MQ=14)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:14795097/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                     >  1:27800685/1‑68 (MQ=255)
               ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                      <  1:34055661/67‑1 (MQ=255)
               ccAGCGTCAATAACGAAGGCGAGAGCACCATCCCCAACGGTGGCACCGTCACCCAGATTAATGGtga                                                      <  1:23193119/67‑1 (MQ=18)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                  >  1:34277943/1‑67 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAc                                                  >  1:34764626/1‑67 (MQ=255)
                   cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                 <  1:15667730/68‑1 (MQ=255)
                    gTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACGCAGATTAAGGGTTAGGACg                                                 >  1:21904035/1‑67 (MQ=14)
                     tCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                <  1:30043717/51‑1 (MQ=11)
                       aaTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAc                                             <  1:5597233/68‑1 (MQ=255)
                        acaaCGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            <  1:21629660/66‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGTTGACGCCACg                                            <  1:35576953/68‑1 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:1359039/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                            >  1:36540067/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTCATGACGCCACg                                            >  1:26493583/1‑68 (MQ=255)
                        aTAACGAAGGCGAGAGCACCAT‑‑CCAACGGTGGCACCTGCACCCAGATTAATGGTGATTACGCCACGGc                                          <  1:12091391/68‑1 (MQ=17)
                         tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGAAGACGCCACg                                            >  1:32436228/1‑67 (MQ=11)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:32888517/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:13111814/68‑1 (MQ=255)
                          aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                          <  1:18368435/68‑1 (MQ=255)
                           aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                         >  1:2840175/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCCCGGCaa                                        >  1:17452427/1‑68 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        <  1:35105001/68‑1 (MQ=255)
                            cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                        <  1:34445704/68‑1 (MQ=255)
                             gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaaa                                       <  1:37381701/68‑1 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     <  1:29107177/68‑1 (MQ=255)
                               aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     <  1:9198398/68‑1 (MQ=255)
                                ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                     >  1:12436994/1‑67 (MQ=255)
                                    agagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                >  1:36084290/1‑68 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                             >  1:32596350/1‑68 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                             >  1:25613062/1‑68 (MQ=255)
                                       gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                             >  1:30832714/1‑68 (MQ=255)
                                        cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACACCGGCa                            >  1:771532/1‑68 (MQ=255)
                                          ccATCACCAACGTTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                          <  1:1190544/68‑1 (MQ=255)
                                            aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                         >  1:4019447/1‑67 (MQ=40)
                                              cACCAACGGTGTCTCCGGCACCCAGATTAATGGTGATGAAGCCACGGCAAACAACAACGGCAAAacca                      <  1:3492767/68‑1 (MQ=14)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacag                   >  1:11537194/1‑66 (MQ=255)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                   >  1:7003644/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                   >  1:9347008/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                   >  1:31366809/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCg                   >  1:27408347/1‑68 (MQ=40)
                                                 cAACGGTGGCACCGGCACCCAGATTAATGGTGAAGACGCCACCGCGAAAAACACGGGAAAAacc                       >  1:5839911/1‑64 (MQ=14)
                                                   aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                 >  1:16029859/1‑68 (MQ=40)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                 >  1:21147781/1‑67 (MQ=40)
                                                    cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:17849735/1‑68 (MQ=40)
                                                     ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                >  1:6885540/1‑67 (MQ=40)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCGAAAAACAGAGGAAAAACCACTGTTGaa              >  1:19385970/1‑67 (MQ=14)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc              >  1:24051440/1‑68 (MQ=40)
                                                      gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc              >  1:8228125/1‑68 (MQ=40)
                                                      gTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAAAAACACAGGAAAAACCACGGTTGag              >  1:5105663/1‑67 (MQ=11)
                                                       tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAGCAACGGCAAAACCACCGTTGACg             >  1:5178565/1‑68 (MQ=255)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             <  1:3799433/67‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             <  1:25412679/67‑1 (MQ=40)
                                                        ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg             >  1:16436858/1‑67 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCGAAAAACACAGGCAAAACCACAGTTGAAggaa          >  1:35816315/1‑65 (MQ=14)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa         >  1:12743359/1‑68 (MQ=40)
                                                           accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa         >  1:6092984/1‑68 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:25521560/68‑1 (MQ=40)
                                                             cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg       <  1:19744512/68‑1 (MQ=40)
                                                              ggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa      <  1:28543355/68‑1 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt    >  1:6574586/1‑67 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt    <  1:10907400/67‑1 (MQ=40)
                                                                 accCAGATTAATGGTGATGACGCCACGGCAAAAAACACAGGAAAAACCACCGTTGACGGAAAGGAtt    >  1:21441353/1‑67 (MQ=14)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAaacaacaa                                 >  1:22778445/1‑36 (MQ=255)
                                                                   ccAGATTAATGGTGATGACGCCACGGCAAAAAACAACGGAAAAACCACGGTTGAAGGAAAGGATTcc  >  1:39052201/1‑67 (MQ=11)
                                                                   |                                                                  
TCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCC  >  CP009273/1460526‑1460659

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: