Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 7.3% 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 7.3% 1870.7 / 19.5 718T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (476/177);  new base C (6/57);  total (482/236)
Fisher's exact test for biased strand distribution p-value = 1.28e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 721 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999533‑999665
                                                                 |                                                                   
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:32490899/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:8753352/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:31999047/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:29519745/68‑1 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:36954242/1‑68 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 <  1:28829168/68‑1 (MQ=255)
     aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTg                                                              >  1:34664523/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAATAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:14056326/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:27802015/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:15314969/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:28252129/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:38293590/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:4782426/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:24385862/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:37642618/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:14975376/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:36173320/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:10000184/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:6931785/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:18471571/68‑1 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCggg                                                       >  1:25500883/1‑68 (MQ=255)
            tGACTAATACCAGGCAGGAATCCTGGCTTATAAACAGTAAAACCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:306730/68‑1 (MQ=255)
               cTAATACCAGTCAGTAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:14374116/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:1765390/1‑68 (MQ=255)
                   tACCCGCCAGGAATCCTGGCTTATAAACAGTAACACCCACCGCCCACCGCGTGGGGCCGGGGGGGAAg                                                <  1:37949349/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCACCGCGTTGGGCCGGGGGGGAAg                                                <  1:25764874/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:2984530/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:13884748/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:17360373/68‑1 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:31350065/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:30985679/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:35839226/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:32638945/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:3533745/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:29297552/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:32936006/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:38889005/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:407822/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:27182002/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:26044375/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:509138/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:34817661/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:10221443/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:14630044/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18207015/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:20448138/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:21029504/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:21752327/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22020606/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:23188354/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:15668161/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:23847637/1‑68 (MQ=255)
                     ccATTCAGTTATCCTCGATTATTACCTGATAAATCACGGGCCCAACGCGTGGGGCCGGGGGGGAAGCg                                              <  1:27657773/68‑1 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:11162169/1‑68 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:9035172/1‑68 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:6368830/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           >  1:29035512/1‑67 (MQ=255)
                          cAGGAATCGTGGCTTATTAACATTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGa                                          <  1:1319571/67‑1 (MQ=255)
                          cAGGAATCCTGGCTTATTAACACTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGa                                          <  1:17940983/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:8058907/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:33537535/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:6097381/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:25095177/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:19227575/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:8381123/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:26926823/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:16374923/67‑1 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                  >  1:20304386/1‑68 (MQ=255)
                                    ttcTTATTAACTGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:2863067/66‑1 (MQ=255)
                                    gggttattaACATTAAAATCAACCGCCCAACGGGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:1992382/65‑1 (MQ=255)
                                    ggCTTATTAACAGTAAACTCCCCCCCCCACCGCTTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:11530506/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:26560129/68‑1 (MQ=255)
                                     gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               >  1:9517784/1‑67 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:1883536/1‑68 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:35573104/1‑68 (MQ=255)
                                       ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:11888645/67‑1 (MQ=255)
                                       ttataaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:12300268/67‑1 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:33986001/68‑1 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:2485932/68‑1 (MQ=255)
                                        tattaACACTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:33351607/68‑1 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:12865250/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:15963450/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:16189435/1‑66 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:36524538/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:22666299/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:6558580/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:18761465/1‑68 (MQ=255)
                                                  aaaaTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                     <  1:7485602/64‑1 (MQ=255)
                                                   aaaTCACCCGCCCACCCCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                <  1:7823166/68‑1 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:17020799/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:23546537/1‑68 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc            <  1:31699832/68‑1 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc            <  1:13378042/68‑1 (MQ=255)
                                                        aaCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           <  1:4407691/68‑1 (MQ=255)
                                                        aaCCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           <  1:5821465/68‑1 (MQ=255)
                                                         aCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:30118524/68‑1 (MQ=255)
                                                          ccccccAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc         <  1:7203943/65‑1 (MQ=255)
                                                            gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCATTACTGGCCGCGCCGCCACtt       >  1:5677955/1‑68 (MQ=255)
                                                              ccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt     <  1:10831867/68‑1 (MQ=255)
                                                               cAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt    <  1:20688929/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 < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: