Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,527 C→A 5.6% pseudogene (879/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,5270CA 5.6% 2762.3 / 8.5 882pseudogene (879/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (519/277);  new base A (74/2);  total (602/280)
Fisher's exact test for biased strand distribution p-value = 4.00e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 898 total aligned reads displayed.

ACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCC  >  CP009273/1460460‑1460594
                                                                   |                                                                   
aCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATc                                                                     >  1:41475847/1‑68 (MQ=255)
aCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATc                                                                     >  1:28507118/1‑68 (MQ=255)
     gcgACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                  <  1:19297163/66‑1 (MQ=255)
       gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt                                                               <  1:22566601/67‑1 (MQ=255)
       gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTg                                                              <  1:33518059/68‑1 (MQ=255)
         cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGa                                                             <  1:13844402/67‑1 (MQ=255)
             ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAgg                                                        <  1:46410092/68‑1 (MQ=255)
               aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGc                                                       <  1:50961129/67‑1 (MQ=255)
               aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATAGACGGTCACAAGGcc                                                      >  1:24135906/1‑68 (MQ=11)
                 aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAt                                                    <  1:19669769/68‑1 (MQ=255)
                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTAGACGGGGACAAGGccctg                                                   >  1:40522203/1‑65 (MQ=11)
                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc                                                   >  1:7027746/1‑68 (MQ=255)
                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc                                                   >  1:12190517/1‑68 (MQ=255)
                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc                                                   >  1:11830134/1‑68 (MQ=255)
                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATc                                                   >  1:34289028/1‑68 (MQ=255)
                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCg                                                  >  1:22879628/1‑67 (MQ=255)
                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCg                                                  >  1:17279587/1‑67 (MQ=255)
                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCg                                                  >  1:26699293/1‑67 (MQ=255)
                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCg                                                  >  1:23231835/1‑67 (MQ=255)
                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGt                                                 <  1:38442266/68‑1 (MQ=255)
                     aGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGTGACAAGGCCATCGTc                                                >  1:10912549/1‑68 (MQ=255)
                       ggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCaa                                              >  1:35194881/1‑68 (MQ=255)
                        gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCaa                                              <  1:20844339/67‑1 (MQ=255)
                        gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAAt                                             >  1:10583976/1‑68 (MQ=255)
                          aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTAGACGGTGACAAGGCCGTCGTaaaaa                                            >  1:26205801/1‑65 (MQ=14)
                             aTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGa                                        >  1:24549442/1‑68 (MQ=255)
                              tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTGAATAACGaa                                       >  1:15459195/1‑68 (MQ=255)
                                aCCGTCACCGATCCGGAGTCCATGGGTATCCAGCTCGACGGTGACAAGGCCCTCGTGAATAACGAAg                                      >  1:16693349/1‑67 (MQ=17)
                                aCCGTCACCGATCCGGAGTCCATGGGTATCCAGCTCGACGGTGACAAGGCCATCGTCAATAACGAAg                                      >  1:47153898/1‑67 (MQ=255)
                                 ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAgg                                     <  1:9012630/67‑1 (MQ=255)
                                 ccGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGc                                    >  1:46094661/1‑68 (MQ=255)
                                  cGTCACCGATCCGGAGTCCATGGGTATCCAGATGGACGGTGACAAGGCCATCGTCAATAACGAAGGCg                                   >  1:31222519/1‑68 (MQ=255)
                                  cGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCg                                   >  1:25586426/1‑68 (MQ=255)
                                  cGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCg                                   >  1:28790720/1‑68 (MQ=255)
                                   gTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACTAGGCCATCGTCAATAACGAAGGCga                                  <  1:48237417/68‑1 (MQ=255)
                                     cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgaga                                <  1:39977710/68‑1 (MQ=255)
                                      aCCGATCCGGAGTCCATGGGTTTACAGATCTACGCTGAAAAGGCCAGCGTCAAAAACGAAGGCgagc                                >  1:30549804/1‑66 (MQ=17)
                                       ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGCCAAGGCCATCGTCAATAACGAAGGCGagag                               >  1:25784325/1‑67 (MQ=255)
                                        cGATCCGGAGTCTATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGc                              >  1:37256807/1‑67 (MQ=255)
                                        cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCa                             <  1:45306754/68‑1 (MQ=255)
                                         gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                            <  1:42438195/68‑1 (MQ=255)
                                           tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAcc                           <  1:33045117/67‑1 (MQ=255)
                                            ccGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCAt                         <  1:21797678/68‑1 (MQ=255)
                                              ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                       <  1:16294416/68‑1 (MQ=255)
                                              ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                       <  1:4090187/68‑1 (MQ=255)
                                              ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                       <  1:9463427/68‑1 (MQ=255)
                                              ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                       <  1:5662391/68‑1 (MQ=255)
                                              ggAGTCCATGGGTAACCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCa                       <  1:13419245/68‑1 (MQ=255)
                                                 gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                    >  1:48609582/1‑68 (MQ=255)
                                                  tCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAct                           <  1:7608248/60‑2 (MQ=255)
                                                  tCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCa                    >  1:6109494/1‑67 (MQ=255)
                                                  tCCATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                        >  1:3970245/1‑63 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:32658509/68‑1 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:42885253/68‑1 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:17758494/68‑1 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:12819566/68‑1 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:47734689/68‑1 (MQ=255)
                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACg                 <  1:37722783/68‑1 (MQ=255)
                                                      tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGt               >  1:49963492/1‑68 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAATggcg              >  1:35994921/1‑65 (MQ=11)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTg              >  1:14217292/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTg              >  1:1996965/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACgggg              >  1:1835612/1‑65 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:36937833/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:28181283/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:44929730/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:45826509/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:7957547/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:10011027/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:35656721/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTg              >  1:36334279/1‑67 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAgggg              >  1:41931715/1‑65 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAggag              >  1:50268549/1‑65 (MQ=255)
                                                        ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAATgggg              >  1:44022714/1‑65 (MQ=14)
                                                        ggTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                        >  1:5167550/1‑57 (MQ=11)
                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcacc         <  1:27445345/68‑1 (MQ=255)
                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcacc         <  1:8482482/68‑1 (MQ=255)
                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccg        >  1:24809590/1‑68 (MQ=255)
                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccg        >  1:29939987/1‑68 (MQ=255)
                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATAAAAAACGGTtgcaccg        >  1:38808839/1‑68 (MQ=18)
                                                              cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccg        >  1:27692082/1‑67 (MQ=255)
                                                              cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccg        >  1:42184842/1‑67 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAACGGGggcaccggc      >  1:14962486/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:29073188/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:29467772/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:6707536/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:49500037/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAAGGGggcaccggc      >  1:39270377/1‑68 (MQ=255)
                                                               aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAAGGCggcaccggc      >  1:347426/1‑68 (MQ=14)
                                                                gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:49038151/1‑67 (MQ=255)
                                                                gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:44430580/1‑67 (MQ=255)
                                                                gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:31578082/1‑67 (MQ=255)
                                                                gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      >  1:18931624/1‑67 (MQ=255)
                                                                gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAACGGAggcaccggc      >  1:20536738/1‑67 (MQ=255)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcac    >  1:25210767/1‑67 (MQ=255)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggca     <  1:11302210/66‑1 (MQ=255)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      <  1:43495885/65‑1 (MQ=255)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggc      <  1:3227269/65‑1 (MQ=255)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAAGGGAggcaccggcac    >  1:46787768/1‑67 (MQ=14)
                                                                  tCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAACGGAggcaccggcac    >  1:36623640/1‑67 (MQ=255)
                                                                   |                                                                   
ACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCC  >  CP009273/1460460‑1460594

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: