Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 10.5% T66P (ACG→CCG)  ycbF → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273999,5980AC10.5% 1220.2 / 29.0 514T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (324/129);  new base C (3/58);  total (327/187)
Fisher's exact test for biased strand distribution p-value = 6.40e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

TACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999542‑999665
                                                        |                                                                   
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCg                                                         >  1:2261301/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:13322426/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:1718191/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:18055019/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:20017258/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:21276326/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:29664537/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:3080942/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:32104671/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:33015344/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:7594437/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:10885296/1‑68 (MQ=255)
 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:9658641/1‑68 (MQ=255)
   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:8029085/1‑68 (MQ=255)
   tGACTAAAACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:10982930/68‑1 (MQ=255)
      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:13004889/1‑68 (MQ=255)
          tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:25237544/68‑1 (MQ=255)
          tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAAACCCCCGCCCACCGCGTGGGGCCGGGGGGGAAg                                                <  1:2584953/68‑1 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGCAGc                                               >  1:24029502/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:5894101/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22441241/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:11552069/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:1206348/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7065800/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:23462215/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:2373353/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:6235623/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:33324862/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:26210842/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:27959915/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:3739727/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:30362477/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:30994390/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:3648781/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:15430456/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:20753248/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18310325/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:19703287/1‑68 (MQ=255)
           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:16576572/1‑68 (MQ=255)
                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:14313792/1‑67 (MQ=255)
                 cAGGAATCCTGTCTTATTAACAGAAAAACCAACCGCCCACCGCGGGGGGCCGGGGGGGAAGCGTCGa                                          <  1:27810702/67‑1 (MQ=255)
                   ggtaTCCTGGCTTTTTAACTGTAATATCAACCTCCCTCCGGGTTGGGCCGGGGGGGAAGCGTCGACAg                                       <  1:5023778/65‑1 (MQ=255)
                     aaTCCTGGCTTATTAACATTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:9017580/68‑1 (MQ=255)
                     aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAg                                       <  1:19408910/66‑1 (MQ=255)
                      aTCCTTGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:17603273/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:3954800/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:16360618/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:9460367/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:1553565/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:20858590/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:4566162/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:29156280/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:25627852/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:14738693/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:2023522/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:2036941/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:13845591/67‑1 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:21708013/1‑67 (MQ=255)
                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAAAGTg                                     >  1:33729785/1‑67 (MQ=255)
                          tGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:25365185/68‑1 (MQ=255)
                           tgCTTATTACCAGTAAAATCAACCGCCCACCGCGTAGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:19179069/67‑1 (MQ=255)
                           ggtttattaACATTAAAATCACCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:14521760/65‑1 (MQ=255)
                           ggattattaACATAAAACTCAACCGCCCAACGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:13435325/65‑1 (MQ=255)
                           ggCTTATTAACATTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:30454107/68‑1 (MQ=255)
                           ggCTTATTAACATTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:24360511/68‑1 (MQ=255)
                           ggCTTATTAACAGTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:2198366/68‑1 (MQ=255)
                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:28710301/68‑1 (MQ=255)
                           ggCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:32531174/68‑1 (MQ=255)
                           ggCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:3553696/68‑1 (MQ=255)
                           gacTTATTAACCGAAAAATAAACCGCCCAACTCTTTGTGCCGGTGGGGAAGCTTCGACAGTGCCAGCa                               <  1:22981497/66‑1 (MQ=255)
                             cTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGCGAATCGTCGACATTGCCAGCAcc                             >  1:26521292/1‑68 (MQ=255)
                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:31577368/1‑68 (MQ=255)
                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:16030912/1‑68 (MQ=255)
                              ttattaACAGTAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:14169581/67‑1 (MQ=255)
                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           >  1:6776302/1‑68 (MQ=255)
                                attaACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:15751319/68‑1 (MQ=255)
                                attaAAAGTAAAATCAACCGCCCCACGGGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:15043078/68‑1 (MQ=255)
                                ataaaCAGTAAACACACCCCCCCAACGCGCGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:31294158/65‑1 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:30120872/1‑66 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:2328563/1‑66 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:4851844/1‑66 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:6032709/1‑66 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:5521845/1‑66 (MQ=255)
                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:10076449/1‑66 (MQ=255)
                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:8405244/1‑68 (MQ=255)
                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:24809198/1‑68 (MQ=255)
                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:27239013/1‑68 (MQ=255)
                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                     >  1:33558228/1‑67 (MQ=255)
                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   >  1:1246486/1‑68 (MQ=255)
                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAATAGTCGACATTCCCAGCACAATTACTGGCAGCTcc              >  1:31746122/1‑68 (MQ=255)
                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCATTACTGGCCGCTcc              >  1:10673332/1‑68 (MQ=255)
                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCATTACGGGCCGCTcc              >  1:26964560/1‑68 (MQ=255)
                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:4359100/1‑68 (MQ=255)
                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:18667764/1‑68 (MQ=255)
                                              cAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc            <  1:3412399/68‑1 (MQ=255)
                                              cAACCGCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc            <  1:8605280/68‑1 (MQ=255)
                                                cccGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:23036622/67‑1 (MQ=255)
                                                aCCGGCCAAGGCTTTGGGGCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcct          <  1:25195691/68‑2 (MQ=255)
                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:28305426/1‑67 (MQ=255)
                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:23179473/1‑67 (MQ=255)
                                                        |                                                                   
TACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999542‑999665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: