Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 999,598 A→C 7.1% 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.1% 2742.2 / 38.6 1037T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (646/298);  new base C (6/81);  total (654/383)
Fisher's exact test for biased strand distribution p-value = 8.40e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1046 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCT  >  CP009273/999533‑999665
                                                                 |                                                                   
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:912693/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:56924302/68‑1 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                 <  1:46217229/68‑1 (MQ=255)
    tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGtt                                                               >  1:29240291/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:30149632/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:36538651/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                             >  1:4599073/1‑68 (MQ=255)
       tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                            <  1:13303086/68‑1 (MQ=255)
       tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                            <  1:43462133/68‑1 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:20944701/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:31978658/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:25711410/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:44201281/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:21419845/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:45586488/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:16067903/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:48207685/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                         >  1:13702742/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCggg                                                       <  1:24354277/68‑1 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                       >  1:955160/1‑68 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGCCCggg                                                       >  1:1754892/1‑68 (MQ=255)
              aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgc                                                                  <  1:53858240/55‑1 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCggggg                                                     <  1:26944390/67‑1 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:51978119/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                    >  1:5870925/1‑68 (MQ=255)
                   tACCTGTCAGGAATCCGGCCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAg                                                <  1:28433775/68‑1 (MQ=255)
                   aaCCAGTCATCAATCCTGGCTTATAAACACTAAAATCACCCGCCCAACGGGTTGGGCCGGGGGGGAAg                                                <  1:56550022/67‑1 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTCACAGTAAAATCAAACGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:28118026/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGTGGGAAGc                                               >  1:29868713/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:7154948/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:35938104/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:32945554/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:8234051/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:32642499/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:34401955/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:34852075/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:40159946/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:4526335/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:5237418/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:10802809/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:47265157/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:16552565/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:12307998/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:12974280/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:19779751/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                               >  1:20182812/1‑68 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:1002222/1‑68 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:54667209/1‑68 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                              >  1:42599586/1‑67 (MQ=255)
                          cAGGAATCCTGGCTTATAAACATTAAACTCAACCGCCCATGGCGTTGGGCCGGGGGGGAAGCGTCGa                                          <  1:35409745/67‑1 (MQ=255)
                             gcaTCCTGTCTTATTAACATTAAAAACAACCCCCCAACCCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                      <  1:5729659/66‑1 (MQ=255)
                              aaTCCTGGCTTATTAACAGAAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:46695193/68‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:1538736/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGCCAGTg                                     >  1:42942597/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:5324968/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:44675349/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:22145652/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:23961315/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:14201431/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                     >  1:7698213/1‑67 (MQ=255)
                               aTCCTGGATAAATAACAGTAAAATCACCCGCCCACCGCGGGGGGCCGGGGGGGAAGCGTCGACAGTg                                     <  1:33683254/67‑1 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                  >  1:14693060/1‑68 (MQ=255)
                                 ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCa                                  >  1:37998197/1‑68 (MQ=255)
                                    ggCTTATTAACATTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:41131522/68‑1 (MQ=255)
                                    ggCTTATTAACATTAAAATCAACCCCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:28920080/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:2341492/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:41793100/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:55647016/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:39463552/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:1856838/68‑1 (MQ=255)
                                    ggCTTATTAACAGAAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:51164625/68‑1 (MQ=255)
                                    gccTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                               <  1:37568538/66‑1 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:22914843/1‑68 (MQ=255)
                                      cTTAGTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                             >  1:49348351/1‑68 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:2513164/68‑1 (MQ=255)
                                        tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:30784884/68‑1 (MQ=255)
                                        tattaACAGAAAAATCACCCCCCCCACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                           <  1:40614014/68‑1 (MQ=255)
                                         attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:57692644/68‑1 (MQ=255)
                                         attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          <  1:47765259/68‑1 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCAtt                          >  1:49790101/1‑66 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                          >  1:6388438/1‑66 (MQ=255)
                                            aaCAGTAAAATCAACCGCCCAACGCGTTGGTCCGGGGGGGAAGAGTCGACAGTGCCATCACCAtt                          >  1:38635921/1‑65 (MQ=255)
                                            aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGACGCGTCGACAGTGCCAGCACCAtt                          >  1:24815839/1‑65 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:26262184/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:49030934/1‑68 (MQ=255)
                                             aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                      >  1:37044234/1‑68 (MQ=255)
                                                gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                   >  1:39072719/1‑68 (MQ=255)
                                                   aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATCACTGGCCGCt                <  1:11205021/68‑1 (MQ=255)
                                                     aTCAACCGCCCAACTCGTTGGGCAGCGGGGGAAGCGTCTACAGTTACAGCACCATTACAGGCCGCTca              >  1:55172583/1‑67 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:19200696/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:52715465/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc              >  1:43913949/1‑68 (MQ=255)
                                                         tccGCCTAACTGGTTGGGACGTGGGGGAAGCGTGGAGAGTGCCAGCACCATTACTGGCCTCTCCGCCa          <  1:56251664/67‑1 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:50255731/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc           >  1:31525380/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa          <  1:14918142/68‑1 (MQ=255)
                                                            gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt       <  1:11656046/68‑1 (MQ=255)
                                                              ccACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt     <  1:27470822/68‑1 (MQ=255)
                                                              ccAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt     <  1:42222478/68‑1 (MQ=255)
                                                                 aCGCGTTGGGCCGGGGGGGAAGCGGCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct  <  1:18163068/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: