Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 6.6% 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 6.6% 1869.4 / 16.7 703T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (477/168);  new base C (2/53);  total (481/222)
Fisher's exact test for biased strand distribution p-value = 2.01e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 706 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999533‑999665
                                                                 |                                                                    
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                                   <  1:43262294‑M1/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                                   <  1:28454690‑M1/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAAGAGTAAAATCAACCGCCCAACGCGTTGc                                                              >  1:21234164‑M1/1‑67 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:41769498/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:12487924/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:12683730/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:28189839/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:33899094/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:33076809/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:32140628/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:36349322/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:42466032/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:19687571/68‑1 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:2619789/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:9300707/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:4278907/68‑1 (MQ=255)
             gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggg                                                      >  1:25440753/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:15211650/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:25672995/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:32966780/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:7845019/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:37925042/1‑68 (MQ=255)
                   tcccACTCAGGAATCCTGGCTTATAAACAGAAAAATCCACCGCCCACCGCGTTGGGCCGGGGGGGAAg                                                <  1:15982938/66‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGTCTTATAAACACAAAAAACAACCCCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:40660792/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:7594523/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:29378467/68‑1 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCCACGCGTTGGGCCGGGGGGGAAGc                                               >  1:10105968/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:26534135/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:40217333/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:10635970/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:3852511/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:37220548/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:2411158/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:1503364/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:23354784/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:6142143/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:17754399/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:6722483/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7217027/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:13566151/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:8231311/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18034590/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAAACAACAGCCCAACGCGTGGGGCCGGGGGGGAAGc                                               >  1:1568505/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGGAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:41259897/1‑68 (MQ=255)
                     ccAGTCAGGAATCCTGTCTTATTAACACAAAAATCCACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              <  1:5915647/68‑1 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGCGCCGGGGGGGAAGCg                                              >  1:14644199/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:3422828/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:3827780/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:19348899/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:17443893/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAACTCACCCGCCCAACGCGTTGGGCggttgattttactgttaataag                                      <  1:20768803‑M1/67‑23 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:27680009/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:33736580/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGGCGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:28734300/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACACTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:22261554/67‑1 (MQ=255)
                               aTCCTGGCTTATAAACAGTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:35823718/67‑1 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGcc                                   >  1:29660431/1‑68 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                  >  1:18439155/1‑68 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                  >  1:27390646/1‑68 (MQ=255)
                                 ccTGGATTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTgaaa                                  >  1:30496618/1‑65 (MQ=255)
                                   tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGc                                >  1:7483354/1‑68 (MQ=255)
                                    ggattattaACAGAAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:24935040/65‑1 (MQ=255)
                                    ggCTTATTACCAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:11683776/68‑1 (MQ=255)
                                    ggCTTATTAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:14247825/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:31632950/68‑1 (MQ=255)
                                    ggCTTATTAACAGGAAAATCACCCGCCCAACGCGTTGGGCCGGTGGGGAAGCGTCGACAGTGCCAGCa                               <  1:978251/68‑1 (MQ=255)
                                     gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               >  1:35469764/1‑67 (MQ=255)
                                      cTTATTAACAGTAAAATCCCCCGCCCACCGCGGGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:27143812/68‑1 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:13940363/1‑68 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:29067217/1‑68 (MQ=255)
                                       ttattaACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:22605831/67‑1 (MQ=255)
                                       ttattaACACTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:4875628/67‑1 (MQ=255)
                                       ttattaAAAGTAAAATCAACCCCCCACCGGGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:26900681/67‑1 (MQ=255)
                                        tattaACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:30974610/68‑1 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:39216805/68‑1 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:2149130/68‑1 (MQ=255)
                                        tattaAAAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:19062987/68‑1 (MQ=255)
                                         attaACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:37618908/68‑1 (MQ=255)
                                         ataaaCAGTAAAATCAACCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:5257641/65‑1 (MQ=255)
                                         ataaaCAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:11982030/65‑1 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:16604515/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:32681297/1‑66 (MQ=255)
                                            aaCAGTAAAATCAACCGCCCAGCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:20445608/1‑65 (MQ=255)
                                            aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:24493007/1‑65 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCggttgattttactggtattaagccaggattccgg                          >  1:12821801‑M1/1‑31 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:20007556/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:11171249/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:38938838/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:35030459/1‑68 (MQ=255)
                                               aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                     >  1:8800432/1‑67 (MQ=255)
                                                   aaaTCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                <  1:30018912/68‑1 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:39741242/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:35179061/1‑68 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:22942835/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:4443966/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:36748358/1‑67 (MQ=255)
                                                            gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       >  1:4207980/1‑68 (MQ=255)
                                                            ccccAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:641579/67‑1 (MQ=255)
                                                              ccAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattagtcagtaaa       <  1:566487‑M1/67‑54 (MQ=255)
                                                              ccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt     <  1:34487290/68‑1 (MQ=255)
                                                                 |                                                                    
TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999533‑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: