Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,590 A→C 5.1% 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.1% 2206.9 / 10.0 700N63T (AAC→ACC) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (485/169);  new base C (3/43);  total (488/212)
Fisher's exact test for biased strand distribution p-value = 1.21e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 717 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCA  >  CP009273/999525‑999657
                                                                 |                                                                   
 gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAAc                                                                    <  1:24813696/66‑1 (MQ=255)
   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGcc                                                                <  1:19525615/68‑1 (MQ=255)
   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGc                                                                 >  1:25047097/1‑67 (MQ=255)
     cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCa                                                              <  1:20952311/68‑1 (MQ=255)
         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                          <  1:18651624/68‑1 (MQ=255)
         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                          <  1:8753352/68‑1 (MQ=255)
          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                         >  1:36954242/1‑68 (MQ=255)
              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAATAGTAAAATCAACCGCCCAACGCGTTgg                                                     >  1:14056326/1‑68 (MQ=255)
              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                     >  1:23668952/1‑68 (MQ=255)
              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                     >  1:21775960/1‑68 (MQ=255)
               tttACTGACTAATACCACTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                     >  1:26344652/1‑67 (MQ=255)
                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                    >  1:10840673/1‑67 (MQ=255)
                 aaCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCACAAAGCGTTGTGcc                                                  >  1:12275794/2‑68 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:10000184/1‑68 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:36173320/1‑68 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:15094211/1‑68 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                 >  1:28252129/1‑68 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGAGTTGGGCCg                                                 >  1:35554447/1‑68 (MQ=255)
                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                               >  1:38214025/1‑68 (MQ=255)
                    tGACTAATACCAGTCAGGAATCCTGGCTTATTAACACTAAAATCAACCGCCCAACGCGTAGGGCCggg                                               >  1:35316564/1‑68 (MQ=255)
                    tGACTAATACCAGGCAGGAATCCTGGCTTATAAACAGTAAAACCAACCGCCCAACGCGTTGGGCCggg                                               <  1:306730/68‑1 (MQ=255)
                       cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                            >  1:15350921/1‑68 (MQ=255)
                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                        <  1:13884748/68‑1 (MQ=255)
                           tACCAGTCAGGAATCCTGGCTTATAAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAg                                        <  1:16025051/68‑1 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:7070354/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGc                                       >  1:28816240/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:25919242/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:29442562/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:11147838/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:25365760/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:30985679/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:31350065/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:38623180/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:32638945/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:32936006/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:37836107/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:34817661/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:35839226/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:10221443/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:2037084/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:13158119/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:23258004/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:15668161/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:16735808/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:23847637/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:20541830/1‑68 (MQ=255)
                            aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                       >  1:8078471/1‑68 (MQ=255)
                             ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                      >  1:9035172/1‑68 (MQ=255)
                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                      >  1:6368830/1‑67 (MQ=255)
                              cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGAGTTGGGCCGGGGGGGAAGCg                                      >  1:19110461/1‑67 (MQ=255)
                                 tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                   >  1:29035512/1‑67 (MQ=255)
                                 tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                   >  1:1994374/1‑67 (MQ=255)
                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAg                               >  1:6536362/1‑68 (MQ=255)
                                     gAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                              >  1:33900675/1‑68 (MQ=255)
                                      aaTCCTTGATAATTAACATAAAAATCACCCGCCCATCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:5525218/68‑1 (MQ=255)
                                      aaTCCTGTATTATTAACAGTAAAATCACCCCCCCCCCGCGGGGGGCCGGGGGGGAAGCGTCGACAGt                              <  1:11457794/67‑1 (MQ=255)
                                      aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:22102429/68‑1 (MQ=255)
                                      aaTCCTGGCTTATTAACAGTAAAATCAACCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:27575852/68‑1 (MQ=255)
                                       aTCGTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:11868099/67‑1 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:4735495/67‑1 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:29916105/1‑67 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:3983338/1‑67 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:6023757/1‑67 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:18420605/67‑1 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:8381123/1‑67 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             >  1:17410538/1‑67 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:27051375/67‑1 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:33503487/67‑1 (MQ=255)
                                       aTCCTGGCTTATAAACAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                             <  1:24293508/67‑1 (MQ=255)
                                         ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:5786924/1‑68 (MQ=255)
                                         ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                          >  1:29700045/1‑68 (MQ=255)
                                           tGGCTTATTAACATTAAAAACAACCCCCCACCGCGGGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGc                        <  1:14381436/68‑1 (MQ=255)
                                            ttcTTATTAACTGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:2863067/66‑1 (MQ=255)
                                            gtcTTATAAACATTAAAATCACCCGCCCCCCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:17849229/66‑1 (MQ=255)
                                            ggCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:33350675/68‑1 (MQ=255)
                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       <  1:26817228/68‑1 (MQ=255)
                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       >  1:9517784/1‑67 (MQ=255)
                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       >  1:14725306/1‑67 (MQ=255)
                                             gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                       >  1:25555630/1‑67 (MQ=255)
                                              ctaattaACAGTAAATTCACCGCCCCCCCGCGTGGGGCCGGGGGGGAAGAGGCGACAGTGCCAGCAcc                     <  1:31671100/65‑1 (MQ=255)
                                              cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     >  1:27786048/1‑68 (MQ=255)
                                               ttataaACAGTAAAACCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                     <  1:31987071/67‑1 (MQ=255)
                                                 attaACAGTAAAATCACCCGCCCACCGGGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  <  1:5020551/68‑1 (MQ=255)
                                                   taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:37523261/1‑66 (MQ=255)
                                                   taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:13418176/1‑66 (MQ=255)
                                                   taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                  >  1:21410379/1‑66 (MQ=255)
                                                    aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTTCCAGCACCAtt                  >  1:18046909/1‑65 (MQ=255)
                                                     aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:22753569/1‑68 (MQ=255)
                                                     aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:17071723/1‑68 (MQ=255)
                                                     aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg              >  1:36524538/1‑68 (MQ=255)
                                                          aaaaTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg             <  1:7485602/64‑1 (MQ=255)
                                                           aaaTCACCCGCCCACCCCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt        <  1:7823166/68‑1 (MQ=255)
                                                           aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt        <  1:32389933/68‑1 (MQ=255)
                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:18817594/1‑68 (MQ=255)
                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc      >  1:6751737/1‑68 (MQ=255)
                                                               cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg     >  1:16266261/1‑67 (MQ=255)
                                                               cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg     >  1:34516514/1‑67 (MQ=255)
                                                                aaCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc   <  1:4407691/68‑1 (MQ=255)
                                                                 aCCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa  <  1:12829913/68‑1 (MQ=255)
                                                                 aCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa  <  1:30118524/68‑1 (MQ=255)
                                                                 |                                                                   
AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCA  >  CP009273/999525‑999657

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: