Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,590 A→C 5.3% N63T (AAC→ACC)  ycbF → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273999,5900AC 5.3% 1468.3 / 5.3 484N63T (AAC→ACC) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (314/137);  new base C (2/31);  total (316/168)
Fisher's exact test for biased strand distribution p-value = 2.09e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 501 total aligned reads displayed.

GAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCA  >  CP009273/999526‑999657
                                                                |                                                                   
                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                  >  1:14489497/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:9658641/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:19589495/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:16758366/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:15677873/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCggg                                               <  1:18653098/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:321856/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:15867741/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:13238315/1‑68 (MQ=255)
                     aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggg                                             >  1:1063885/1‑68 (MQ=255)
                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggtggg                                            <  1:25523015/68‑1 (MQ=255)
                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCcgggcg                                            >  1:24829073/1‑66 (MQ=255)
                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGa                                          >  1:14067054/1‑68 (MQ=255)
                          tACCAGTCCGGAATCCTGTCTTAATAACACTAAAATCACCCGCCCAACGCGTGGGGCCGGGGGGGAAg                                        <  1:18212453/68‑1 (MQ=255)
                          tACCAGTCAGGATTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                        <  1:15373216/68‑1 (MQ=255)
                          tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                        <  1:6225557/68‑1 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:24393384/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:23250085/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:27874926/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:2373353/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:29543149/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:20605980/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:30641562/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:3142589/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:10875511/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:23435246/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:33571129/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:33915380/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:3651591/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:8043655/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:7639308/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:5894101/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGATGGGAAGc                                       >  1:14916251/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCACGGGGGGAAGc                                       >  1:20965572/1‑68 (MQ=255)
                            ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                      >  1:28996965/1‑68 (MQ=255)
                            ccAGTCAGGAATCCTGGCTTAATAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCg                                      <  1:17485855/68‑1 (MQ=255)
                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                   >  1:28231713/1‑67 (MQ=255)
                                   ggtaTCCTGGCTTTTTAACTGTAATATCAACCTCCCTCCGGGTTGGGCCGGGGGGGAAGCGTCGACAg                               <  1:5023778/65‑1 (MQ=255)
                                     aaTCCTGGCTTATTAACATTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:9017580/68‑1 (MQ=255)
                                     aaTCCTGGCTTATTAACATTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:15494914/68‑1 (MQ=255)
                                     aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                              >  1:3811366/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGc                            <  1:4761270/68‑1 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:30415660/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:34168512/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:9460367/67‑1 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:26811172/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:240967/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:1643541/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:11387950/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:33474323/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:2269335/67‑1 (MQ=255)
                                      aTCCTGGCTTATTAACACTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:3575476/67‑1 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:17671212/1‑68 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:33015966/1‑68 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:32575837/1‑68 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:22990298/1‑68 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:20150023/1‑68 (MQ=255)
                                           gtcTTATTAACACAAAAATCAACCCCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:2863635/66‑1 (MQ=255)
                                           gtattattaACACTAAAAACAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:1242299/65‑1 (MQ=255)
                                           ggCTTATTAACAGTAAACCCAACCGCCCACCGGGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:25992628/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:2198366/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:25359120/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:19088922/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:25214389/68‑1 (MQ=255)
                                           ggCTTATTAACAGAAAAATCACCCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:13404615/68‑1 (MQ=255)
                                           ggCTTATAAACACTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:3957482/68‑1 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     >  1:3189705/1‑68 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     >  1:19983882/1‑68 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     >  1:21347870/1‑68 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGACAGCAcc                     >  1:14605119/1‑68 (MQ=255)
                                              ttattaACATTAAAATCAACCGCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:21740881/67‑1 (MQ=255)
                                              ttattaACATTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:1293469/67‑1 (MQ=255)
                                              ttattaACAGTAAATTCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:20302491/67‑1 (MQ=255)
                                              ttattaACAGTAAAATCAACCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:20767876/67‑1 (MQ=255)
                                              ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:13824058/67‑1 (MQ=255)
                                               tattaACCCTAAAATCAACCCCCCCCCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                   <  1:16163755/68‑1 (MQ=255)
                                               tataaACATTAAAACCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                   <  1:7142974/68‑1 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:4851844/1‑66 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:11964742/1‑66 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:6799390/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)
                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:10967343/1‑65 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:27159686/1‑68 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:8405244/1‑68 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             >  1:3163498/1‑67 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             >  1:31166780/1‑67 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             >  1:29831051/1‑67 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             >  1:14250729/1‑67 (MQ=255)
                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc           >  1:30901674/1‑68 (MQ=255)
                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:19226938/1‑68 (MQ=255)
                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:3015875/1‑68 (MQ=255)
                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:27441609/1‑68 (MQ=255)
                                                              cAACCGCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc    <  1:8605280/68‑1 (MQ=255)
                                                              cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc    <  1:23155600/68‑1 (MQ=255)
                                                               aaCCCCCCCCCGCGGGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   <  1:26416921/68‑1 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   >  1:22441588/1‑67 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   >  1:4419465/1‑67 (MQ=255)
                                                                aCCCCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa  <  1:25803297/68‑1 (MQ=255)
                                                                |                                                                   
GAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCA  >  CP009273/999526‑999657

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: