Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 8.9% 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.9% 1467.9 / 22.9 563T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (357/147);  new base C (5/53);  total (362/201)
Fisher's exact test for biased strand distribution p-value = 9.07e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 573 total aligned reads displayed.

GATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999531‑999665
                                                                   |                                                                   
  tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                    >  1:6178503/1‑67 (MQ=255)
   atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAAGAGTAAAATCAACCGCCCAAcgc                                                                  <  1:18415229/68‑1 (MQ=255)
   atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:14948200/68‑1 (MQ=255)
    tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 <  1:32620761/68‑1 (MQ=255)
        ttttACTGCCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             <  1:33286502/68‑1 (MQ=255)
        ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:16460259/1‑68 (MQ=255)
        ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:27858271/1‑68 (MQ=255)
         tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                            <  1:24046588/68‑1 (MQ=255)
          ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                            >  1:4870958/1‑67 (MQ=255)
           tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                          >  1:29122527/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:32469379/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:30529943/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:23597075/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:8692019/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:14411964/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:25850762/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:16172002/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:7016761/1‑68 (MQ=255)
            aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:12143566/1‑68 (MQ=255)
              tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCggg                                                       <  1:7698897/68‑1 (MQ=255)
              tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:17935688/68‑1 (MQ=255)
              tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:11319156/68‑1 (MQ=255)
              tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:12674634/1‑68 (MQ=255)
                 cTAATACCAGTCAGGAATCCTGGCTTATTTACAGTAAAATCAACCGCCCAACGCTTTGTGCCgggggg                                                    >  1:32100477/1‑68 (MQ=255)
                     tACCTGTCAGTAATCCTTTCTTATTAACAGTAAAAACAACCGCCCCACGCGTTGGGCCGGGGGGGAAg                                                <  1:12904379/68‑1 (MQ=255)
                     tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:24175484/68‑1 (MQ=255)
                     tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:31239280/68‑1 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:26864158/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:27024521/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:34591691/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:417646/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:24490434/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:10132387/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:29308340/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:4351131/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:33690564/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22598486/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22069459/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:21430451/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:17053127/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18599680/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:5567383/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:20106526/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18917883/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7210952/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:9019365/1‑68 (MQ=255)
                      aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTCCGCCGGGGGGGAAGc                                               >  1:18251970/1‑68 (MQ=255)
                        cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:31785149/1‑67 (MQ=255)
                        cAGTCAGGAACCCTGGCTTATAAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGt                                             <  1:6654919/68‑1 (MQ=255)
                           tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:10982447/1‑67 (MQ=255)
                           tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:27551949/1‑67 (MQ=255)
                           tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:20756014/1‑67 (MQ=255)
                           gcAGGAATCCTCGCTTATTAACAGAAAAATCAACCGACCAACGCGTGGGGCCGGGGGGGAAGCGGCGa                                          <  1:22431076/67‑1 (MQ=255)
                            cAGTAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCg                                           <  1:922557/66‑1 (MQ=255)
                                aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      >  1:5443222/1‑67 (MQ=255)
                                 ttCCTTGATTATAAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:23278440/66‑1 (MQ=255)
                                 aTCCTGGCTTATTACCTGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:15505325/67‑1 (MQ=255)
                                 aTCCTGGCTTATTAACATTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:33999773/67‑1 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:24746471/1‑67 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:14042399/1‑67 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:30835151/67‑1 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:6510303/1‑67 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:9557697/67‑1 (MQ=255)
                                 aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:19021664/67‑1 (MQ=255)
                                 aTCCTGGATTATAAACATAAAAATCAACCCCCCAAGGGTGTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:13324176/67‑1 (MQ=255)
                                      ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:29907086/68‑1 (MQ=255)
                                      ggCTTATTAACAGTAAAAACAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:26536136/68‑1 (MQ=255)
                                      ggCTTATAAACATTAAAATCAACCGCCCCACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:2882015/68‑1 (MQ=255)
                                      ggCTTATAAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:30144873/68‑1 (MQ=255)
                                        cTTATTAANAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:5842207/1‑68 (MQ=255)
                                        cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:3797249/1‑68 (MQ=255)
                                        cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:752139/1‑68 (MQ=255)
                                        cTTATTAACACTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:10458354/66‑1 (MQ=255)
                                         ttattaACAGTAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:19681388/67‑1 (MQ=255)
                                         ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCa                            <  1:17371830/68‑1 (MQ=255)
                                         ttataaACAGTAAAATCCACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:8216879/67‑1 (MQ=255)
                                          tattaACAGTAAAATCACCCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:19202913/68‑1 (MQ=255)
                                          tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:16925556/68‑1 (MQ=255)
                                          tattaACAGTAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:29654693/68‑1 (MQ=255)
                                          tataaACATTAAAATCAACCGCCCAACGGTTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:20305427/68‑1 (MQ=255)
                                           attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:15136175/68‑1 (MQ=255)
                                             taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:9786605/1‑66 (MQ=255)
                                              aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:16715133/1‑65 (MQ=255)
                                               aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACGGTGCCAGCCCCATTACTg                      >  1:20964253/1‑68 (MQ=255)
                                               aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:15988236/1‑68 (MQ=255)
                                                   tAAAATCAACCGCCCAACGCGTTGGGCAGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   >  1:25120885/1‑67 (MQ=255)
                                                    aaaaTCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   <  1:12787863/66‑1 (MQ=255)
                                                       aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:2386809/1‑68 (MQ=255)
                                                       aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:19508867/1‑68 (MQ=255)
                                                       aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:3538358/1‑68 (MQ=255)
                                                       aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:32987745/1‑68 (MQ=255)
                                                           aCCGCCCACCGCGCTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:11700354/68‑1 (MQ=255)
                                                           aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:7962365/1‑67 (MQ=255)
                                                           aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCACCACCATTACTGTCCGCTCCGcc           >  1:27285002/1‑67 (MQ=255)
                                                            ccGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc         <  1:28278898/68‑1 (MQ=255)
                                                             cGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt        >  1:21687620/1‑68 (MQ=255)
                                                              gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCAATACAGGCCCCTCCGCCAcat       >  1:14109148/1‑66 (MQ=255)
                                                              ccccACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:26694956/67‑1 (MQ=255)
                                                                ccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt     >  1:5076317/1‑68 (MQ=255)
                                                                   aCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct  <  1:13731269/68‑1 (MQ=255)
                                                                   |                                                                   
GATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999531‑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: