Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,590 A→C 5.6% 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.6% 2218.0 / 14.0 703N63T (AAC→ACC) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (480/174);  new base C (3/46);  total (483/220)
Fisher's exact test for biased strand distribution p-value = 1.31e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 729 total aligned reads displayed.

GAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCA  >  CP009273/999526‑999657
                                                                |                                                                   
gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAAcc                                                                   >  1:42739114/1‑67 (MQ=255)
  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGC                                                                 <  1:19358774‑M1/67‑1 (MQ=255)
     gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCA                                                              >  1:23248337‑M1/1‑67 (MQ=255)
     gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCA                                                              >  1:15013381‑M1/1‑67 (MQ=255)
       tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaac                                                            <  1:10592222‑M2/67‑8 (MQ=255)
        atCTATTTTACTGTCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCTCCCACCGC                                                          <  1:15157754‑M1/68‑1 (MQ=255)
        atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaacgc                                                          >  1:6855023‑M2/1‑59 (MQ=255)
        atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                          <  1:43262294‑M1/68‑1 (MQ=255)
         tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaacgcg                                                         >  1:29971854‑M2/1‑58 (MQ=255)
         tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaacgcg                                                         <  1:36166903‑M2/68‑11 (MQ=255)
          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGT                                                        >  1:17089667‑M1/1‑68 (MQ=255)
           tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                          <  1:21614871‑M1/65‑1 (MQ=255)
            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgccccacgcgttg                                                      <  1:24283572‑M2/68‑14 (MQ=255)
             ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaacgcgttgg                                                     <  1:39312885‑M2/68‑15 (MQ=255)
             ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGG                                                     >  1:32741496‑M1/1‑68 (MQ=255)
             ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGG                                                     >  1:7520929‑M1/1‑68 (MQ=255)
              tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCgcccaacgcgttggg                                                    >  1:13433927‑M2/1‑53 (MQ=255)
              tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGG                                                     >  1:7284657‑M1/1‑67 (MQ=255)
                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                  >  1:27660718/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCg                                                 <  1:37630335/68‑1 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:27061325/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:28452792/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:32040005/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:6007150/1‑68 (MQ=255)
                 aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCAAACGCGTTGGGCCg                                                 >  1:1553/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:3508555/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:38961375/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:39875560/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:28707666/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:37233095/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:42466032/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:9300707/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:19687571/68‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:35379579/1‑68 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               <  1:8291511/68‑1 (MQ=255)
                   tGACTAAAACCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCAACGCGTTGGGCCggg                                               <  1:12436761/68‑1 (MQ=255)
                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                            <  1:22263056/68‑1 (MQ=255)
                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                            >  1:37925042/1‑68 (MQ=255)
                          tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAg                                        <  1:11907125/68‑1 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCCACGCGTTGGGCCGGGGGGGAAGc                                       >  1:10105968/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAGCGCGTTGGGCCGGGGGGGAAGc                                       >  1:17850204/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:34675063/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:33393551/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:4141516/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:28122344/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:38659034/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:40721668/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:2461935/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:41974959/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:4368218/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:811170/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:8860808/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:16097497/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:13799987/1‑68 (MQ=255)
                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:14579680/1‑68 (MQ=255)
                           aCCACTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGACCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:18393591/1‑68 (MQ=255)
                            ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                      >  1:32412455/1‑68 (MQ=255)
                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                   >  1:40455953/1‑67 (MQ=255)
                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                   >  1:3827780/1‑67 (MQ=255)
                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAc                                 <  1:1154980/68‑1 (MQ=255)
                                 cAGGAATCCTGGCTTATTAACAGAAAAAACACCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGa                                  <  1:25153858/67‑1 (MQ=255)
                                     caTCCTGCCTTATTAACAGTAAATTCAACCCCCCAACCCGGGGGGCCGGGGGGGAAGCGTCGACAGt                              <  1:11085723/66‑1 (MQ=255)
                                     aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                              <  1:21230912/67‑1 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAACTCACCCGCCCAACGCGTTGGGCggttgattttactgttaataag                             <  1:20768803‑M1/67‑23 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:7685279/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:5679963/1‑67 (MQ=255)
                                      aTCCTGGCTTATTAACACTAAAACCACCCCCCCACCGCGGGGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:35625066/67‑1 (MQ=255)
                                      aTCCTGGATTATTAACATTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:5244311/67‑1 (MQ=255)
                                           gtcTTATAAACATAAAAAACAACCGCCCACCGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:25507797/66‑1 (MQ=255)
                                           ggattattaACACAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:3423129/65‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:14117724/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:17568445/68‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:22879174/68‑1 (MQ=255)
                                           ggCTTATTAACAGGAAAATCACCCGCCCAACGCGTTGGGCCGGTGGGGAAGCGTCGACAGTGCCAGCa                       <  1:978251/68‑1 (MQ=255)
                                           ggCTTATAAACATTAAAATCCCCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:26812313/68‑1 (MQ=255)
                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       >  1:26354330/1‑67 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCAAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     >  1:27321123/1‑68 (MQ=255)
                                              ttattaACATTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:36749972/67‑1 (MQ=255)
                                              ttattaACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:22605831/67‑1 (MQ=255)
                                              ttattaACACTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:4875628/67‑1 (MQ=255)
                                              ttataaACAGAAAAATCAACCCCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCa                    <  1:30730851/68‑1 (MQ=255)
                                               tattaAAAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                   <  1:19062987/68‑1 (MQ=255)
                                               tataaACAGAAAAATCACCCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                   <  1:32968805/68‑1 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:12674631/1‑66 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:25915903/1‑66 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:29667118/1‑66 (MQ=255)
                                                  taACAGTAAAAACAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGccaccaccaat                  >  1:33765826/1‑64 (MQ=255)
                                                   aaCAGTAAAATCAACCGCCCAGCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:20445608/1‑65 (MQ=255)
                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:23913519/1‑65 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCggttgattttactggtattaagccaggattccgg                 >  1:12821801‑M1/1‑31 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:30262736/1‑68 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:20007556/1‑68 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             >  1:655935/1‑67 (MQ=255)
                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:29057413/1‑68 (MQ=255)
                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:33943993/1‑68 (MQ=255)
                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGAGCCGTCGAAAGTGCCAGCACCATTACGGGCCCCTCCGc    >  1:31319102/1‑67 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   >  1:6200422/1‑67 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   >  1:36437326/1‑67 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   >  1:15925340/1‑67 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa  <  1:29341756/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: