Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 6.0% 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.0% 1313.6 / 6.7 481T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (315/128);  new base C (2/33);  total (317/164)
Fisher's exact test for biased strand distribution p-value = 7.14e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 484 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTC  >  CP009273/999533‑999664
                                                                 |                                                                  
tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                   >  1:18101095/1‑67 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                 <  1:23942521/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                 <  1:4529283/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                 <  1:10019106/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                            >  1:15251807/1‑68 (MQ=255)
        ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                           >  1:14339129/1‑67 (MQ=255)
        ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGc                                                          <  1:12546392/68‑1 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:5005411/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:18763478/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:14177688/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:13872500/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:13568719/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:28408960/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:28012556/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:1150782/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:27664045/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:26523898/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:17677703/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:24696712/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCCg                                                        >  1:14006506/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                      >  1:14908223/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                   >  1:21924935/1‑68 (MQ=255)
                   tACCAGTCAGTAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                               <  1:6623039/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGGAAg                                               <  1:8682649/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATAAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAg                                               <  1:23315206/68‑1 (MQ=255)
                   tACCAGTCAGGAATCCGGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAg                                               <  1:21580429/68‑1 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:2935440/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:21513114/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:21141437/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:24406034/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:25126760/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:25348275/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:22532592/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:3391653/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:4586057/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:6117428/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:6869948/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:70891/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:9569134/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:9831230/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:16857357/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:11677755/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:12173660/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:17337945/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:19512814/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:12961507/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:13343157/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:15407547/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGACCAACGCGTTGGGCCGGGGGGGAAGc                                              >  1:20082611/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACACTAAAATCAACCGCCCAACGCGTTGGCCCGGGGGGGAAGc                                              >  1:7869194/1‑68 (MQ=255)
                     cccgTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                             <  1:14713475/65‑1 (MQ=255)
                     ccAGTCAGGAATCCTTGCTTATTAACACTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCg                                             <  1:17080419/68‑1 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                          >  1:10811060/1‑67 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAc                                        >  1:3852879/1‑68 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAc                                        >  1:18887640/1‑68 (MQ=255)
                              aaTCCTTGCTTATTAACAGTAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    <  1:11281519/68‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                     <  1:19809212/67‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAAACAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                     <  1:26022012/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAACTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    <  1:16073146/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    >  1:20461247/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    >  1:21806735/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    >  1:27166731/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    >  1:26741836/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                    >  1:16354559/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                    <  1:10314029/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAAACAACCCCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                    <  1:15797055/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGAAAAATCACCCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGc                                   <  1:261957/68‑1 (MQ=255)
                               aTCCTGGATTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTg                                    <  1:10611650/67‑1 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                 >  1:9325703/1‑68 (MQ=255)
                                  cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGGCGGGGGGGAAGCGTCGACAGTGCCAg                                <  1:11016775/68‑1 (MQ=255)
                                    ggCTTATTAACCGTAAATTCCACCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              <  1:8410477/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              <  1:20916856/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              <  1:9177985/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              <  1:6764622/68‑1 (MQ=255)
                                    ggCTTATTAACAGAAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              <  1:2753188/68‑1 (MQ=255)
                                     gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              >  1:4127664/1‑67 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                            >  1:1892152/1‑68 (MQ=255)
                                       ttatttTCAGTTAAATCAACCGGCAATCGGGGTGGGGCGGGGGGGAAGCGTCGACAGTGCGAGCAcc                            <  1:9036971/67‑1 (MQ=255)
                                       ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCa                           >  1:26347900/1‑68 (MQ=255)
                                       ttattaACAGAAAAATCACCCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                            <  1:8536169/67‑1 (MQ=255)
                                        tattaACAGTAAAAACCCCCGCCCCAGGCGGGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                          <  1:17945058/68‑1 (MQ=255)
                                        tattaACAGTAAAAACAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                          <  1:3108874/68‑1 (MQ=255)
                                        taataACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                          <  1:2314399/65‑1 (MQ=255)
                                         attaACAGTAAAATCACCCGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                         <  1:16505247/68‑1 (MQ=255)
                                         ataaaCAGTAAAATCAACCCCCCACCGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                         <  1:548018/65‑1 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                              >  1:28649839/1‑61 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                         >  1:28896028/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                         >  1:7477471/1‑66 (MQ=255)
                                               aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                    >  1:22127716/1‑67 (MQ=255)
                                                gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                  >  1:14491286/1‑68 (MQ=255)
                                                   aaaTCACCCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt               <  1:27809995/68‑1 (MQ=255)
                                                   aaaTCAACCGCCCCACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt               <  1:23666485/68‑1 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc             >  1:5862850/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc             >  1:20778643/1‑68 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGCCAGTGCCAGCACCATTACTGGCCGCTCCg            >  1:24805647/1‑67 (MQ=255)
                                                         aCCTCCCATCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa         <  1:22850744/68‑1 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc          >  1:15002640/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc          >  1:223026/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc          >  1:20317387/1‑67 (MQ=255)
                                                              ccACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt    <  1:24281980/68‑1 (MQ=255)
                                                                 |                                                                  
TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTC  >  CP009273/999533‑999664

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: