Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 8.4% 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,5980AC 8.4% 1390.0 / 21.6 551T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (347/150);  new base C (4/50);  total (351/200)
Fisher's exact test for biased strand distribution p-value = 9.39e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 557 total aligned reads displayed.

ATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999532‑999665
                                                                  |                                                                   
atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                    >  1:11960974/1‑68 (MQ=255)
 tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                    >  1:14516645/1‑67 (MQ=255)
 tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                    <  1:395503/67‑1 (MQ=255)
  atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:16186108/68‑1 (MQ=255)
  atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:19825063/68‑1 (MQ=255)
  atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:2278613/68‑1 (MQ=255)
  atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:27400673/68‑1 (MQ=255)
  atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:25372706/68‑1 (MQ=255)
   tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:16954172/1‑68 (MQ=255)
   tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:5188443/1‑68 (MQ=255)
       ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             <  1:11807204/68‑1 (MQ=255)
       ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:3592887/1‑68 (MQ=255)
         ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGc                                                           >  1:8769497/1‑68 (MQ=255)
           aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcgg                                                         >  1:1939805/1‑66 (MQ=255)
           aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:29232925/1‑68 (MQ=255)
           aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:18279698/1‑68 (MQ=255)
           aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:9002705/1‑68 (MQ=255)
           aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:20198429/1‑68 (MQ=255)
             tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:29508613/1‑68 (MQ=255)
             tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:24941298/68‑1 (MQ=255)
             tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:26935397/68‑1 (MQ=255)
             tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggc                                                       >  1:11203746/1‑67 (MQ=255)
                    tACCAGTCAGGAATCCGGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:9741228/68‑1 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGc                                               >  1:27256812/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:15835995/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:5373351/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:12066391/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7999565/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:462933/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22969323/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:12947720/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:24409424/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18871405/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:2472259/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:6066630/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:9321854/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7148557/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGCCCGGGGGGGACGc                                               >  1:10032411/1‑68 (MQ=255)
                     aCCAGTCAGGAATCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTCGGCCGGGGGGGaacc                                               >  1:12773928/1‑66 (MQ=255)
                      ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:15998922/1‑68 (MQ=255)
                      ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:5556783/1‑68 (MQ=255)
                      ccACTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGGGTTGGGCCGGGGGGGAAGCg                                              <  1:10802880/68‑1 (MQ=255)
                       cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:25643048/1‑67 (MQ=255)
                       cAGTCAGGAATCCTGGCTTAATAACAGTAAAATCAACCGCCAAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:28569154/1‑67 (MQ=255)
                         gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:5795185/1‑68 (MQ=255)
                          tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:6621279/1‑67 (MQ=255)
                          tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:15284494/1‑67 (MQ=255)
                           cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           <  1:17998940/66‑1 (MQ=255)
                             ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAg                                       >  1:19108648/1‑68 (MQ=255)
                             ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAg                                       >  1:17505066/1‑68 (MQ=255)
                               aaTCCTGTCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:21187611/68‑1 (MQ=255)
                               aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:21441689/68‑1 (MQ=255)
                               aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGTCAGt                                      >  1:18611586/1‑67 (MQ=255)
                               aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAg                                       <  1:31210875/66‑1 (MQ=255)
                               aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      >  1:30148398/1‑67 (MQ=255)
                               aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      >  1:14086778/1‑67 (MQ=255)
                               aaTCCTGGCTTATAAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:14675703/68‑1 (MQ=255)
                               aaTCCTGGCTTATAAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGt                                      <  1:24222095/67‑1 (MQ=255)
                                aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:30612783/1‑67 (MQ=255)
                                 tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGcc                                   >  1:29917148/1‑68 (MQ=255)
                                     ggCTTATTAACATTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:31973452/68‑1 (MQ=255)
                                     ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGc                                >  1:28095446/1‑67 (MQ=255)
                                     ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:28951054/68‑1 (MQ=255)
                                     ggCTTATTAACAGTAAAAGCAACCGCCCCCCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:24101136/68‑1 (MQ=255)
                                     ggCTTATTAACAGTAAAACCCCCCGCCCACCGGGGGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:10309191/68‑1 (MQ=255)
                                     ggCTAATAACCATTAAAATCACCCCCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:4326869/68‑1 (MQ=255)
                                      gCTTATTAACAGTAAAAACAACCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:7063062/67‑1 (MQ=255)
                                       cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:31735844/1‑68 (MQ=255)
                                       cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:4820264/1‑68 (MQ=255)
                                        ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:7533732/67‑1 (MQ=255)
                                         tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:2373820/68‑1 (MQ=255)
                                         tataaACATAAAAACCACCCCCCCCCCGGGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:13285235/68‑1 (MQ=255)
                                         tataaACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:11095333/68‑1 (MQ=255)
                                          attaACATTAAAATCAACAGACCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:15507745/68‑1 (MQ=255)
                                            taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:30973122/1‑66 (MQ=255)
                                            taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:31330299/1‑66 (MQ=255)
                                            taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:5054131/1‑66 (MQ=255)
                                            taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACAAtt                          >  1:6345064/1‑66 (MQ=255)
                                             aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:20691589/1‑65 (MQ=255)
                                              aCAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:8449683/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:7727040/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:16680981/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:17235651/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:12325432/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:27765706/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:1140757/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:10716224/1‑68 (MQ=255)
                                              aCAGTAAAATCAACCGCCCAACGCGTTGGGCCCGGGGGGAAGCGTCGACACTGCCAGCACCATTAcag                      >  1:23311350/1‑66 (MQ=255)
                                                aGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                     >  1:3316060/1‑67 (MQ=255)
                                                aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                     >  1:21781143/1‑67 (MQ=255)
                                                    aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                <  1:22502725/68‑1 (MQ=255)
                                                      aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:16422103/1‑68 (MQ=255)
                                                      aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:24542752/1‑68 (MQ=255)
                                                        cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                         >  1:15036574/1‑55 (MQ=255)
                                                        cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg             >  1:26597569/1‑67 (MQ=255)
                                                        cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg             >  1:20963315/1‑67 (MQ=255)
                                                          aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:23493631/1‑67 (MQ=255)
                                                           ccGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:28224500/1‑66 (MQ=255)
                                                             gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       >  1:14767563/1‑68 (MQ=255)
                                                             gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:13698487/68‑1 (MQ=255)
                                                                  |                                                                   
ATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999532‑999665

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