Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 5.7% 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 5.7% 2386.2 / 19.8 888T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (566/255);  new base C (6/56);  total (574/314)
Fisher's exact test for biased strand distribution p-value = 8.89e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 894 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999533‑999665
                                                                 |                                                                   
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:28058705/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:39463997/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:10010505/68‑1 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:42611098/1‑68 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:38740584/1‑68 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 >  1:4870418/1‑68 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 <  1:21779280/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             <  1:11913528/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:32599366/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:43011866/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:45944856/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:33115892/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:8442088/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:17470384/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:31150971/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:51190825/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:37857167/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:24468857/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:11092991/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:238840/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:46882544/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:50239718/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       <  1:38258064/68‑1 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggg                                                     <  1:26320713/67‑1 (MQ=255)
                   aaCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                <  1:33600331/67‑1 (MQ=255)
                    aCCCGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:49705661/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAATTCAACCGCCCAACGAGTTGGGCCGGGGGGGAAGc                                               >  1:34580468/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:49231844/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:22568340/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                                >  1:41280303/1‑67 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:43457478/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:32076980/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:41907877/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:26811012/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:27804926/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:27049263/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:14583352/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:1274725/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:18055536/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:19632008/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:21225855/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7743756/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7268515/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAGCGTTGGGCCGGGGGGGAAGCg                                              >  1:39943481/1‑68 (MQ=255)
                    aCCAGTCAGGAATACTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:11386894/1‑68 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              <  1:16055878/68‑1 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:36094383/1‑67 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:5436318/1‑67 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                           <  1:181864/66‑1 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAc                                         >  1:44469388/1‑68 (MQ=255)
                              taTCCTTGATTATTAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      <  1:48142393/66‑1 (MQ=255)
                              aaTCGTGGCTTATTAACAGTAAAATCAACCGCCCATCGCGGTGGGCCGGGGGGGAAGCGGCGACAGTg                                     <  1:18921746/68‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      >  1:15268016/1‑67 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:13648988/68‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:44981408/68‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:36724287/68‑1 (MQ=255)
                              aaTCCTGGCTTATAAAACGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:50535337/68‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:40513412/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:9200885/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:41649039/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:42182570/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:47268503/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:7001537/67‑1 (MQ=255)
                               aTCCTGGCTTATAAACAGAAAAAACAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGc                                    <  1:35737300/68‑1 (MQ=255)
                                tcgtGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGcc                                   <  1:27293195/65‑1 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGCCAGTGCCa                                  >  1:34235170/1‑68 (MQ=255)
                                   tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGCCAGTGCCAGc                                >  1:16979623/1‑68 (MQ=255)
                                    tgCTAATTAACTGTAAAATCTACTGCCATACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:5740506/67‑1 (MQ=255)
                                    ggCTTATTAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:30447003/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:16643290/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:46435709/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:14405902/68‑1 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:25381384/1‑68 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:24816594/1‑68 (MQ=255)
                                       ttattaACACAAAAAACAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             <  1:10468781/67‑1 (MQ=255)
                                         attaACATTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:28619268/68‑1 (MQ=255)
                                         attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:8780588/68‑1 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:35443212/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:23001234/1‑66 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:30708245/1‑68 (MQ=255)
                                             aCAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      <  1:31658369/68‑1 (MQ=255)
                                                gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   >  1:5950780/1‑68 (MQ=255)
                                                gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   >  1:20006700/1‑68 (MQ=255)
                                                   aaaTCACCCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                <  1:37010961/68‑1 (MQ=255)
                                                   aaaTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                <  1:33882108/68‑1 (MQ=255)
                                                     aTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:9446398/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:47749734/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:36355735/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:29171739/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:13110465/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACAATTACTGGCCGCTcc              >  1:34839138/1‑68 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg             >  1:12331778/1‑67 (MQ=255)
                                                         cccGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:2327381/67‑1 (MQ=255)
                                                         aCCGCCCACCGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:20472913/68‑1 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGATGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:37315140/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:6588649/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:25587616/1‑67 (MQ=255)
                                                           cccccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt        <  1:22224316/66‑1 (MQ=255)
                                                            gCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:29670168/68‑1 (MQ=255)
                                                            ccccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:44044981/67‑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: