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% 1380.1 / 9.3 506T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (336/132);  new base C (3/35);  total (339/167)
Fisher's exact test for biased strand distribution p-value = 3.29e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 508 total aligned reads displayed.

TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAT  >  CP009273/999533‑999662
                                                                 |                                                                
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTACCAGTAAAATCAACCGCCCAACGC                                                               <  1:15889642‑M1/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                               <  1:25522612‑M1/68‑1 (MQ=255)
 atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGC                                                               <  1:6421481‑M1/68‑1 (MQ=255)
  tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCG                                                              <  1:11778925‑M1/68‑1 (MQ=255)
   cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGT                                                             <  1:25293179‑M1/68‑1 (MQ=255)
   cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGT                                                             >  1:14996503‑M1/1‑68 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCACCGCGTTGG                                                          <  1:3286303‑M1/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGG                                                          <  1:3744929‑M1/68‑1 (MQ=255)
      ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGG                                                          <  1:1523568‑M1/68‑1 (MQ=255)
       tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGG                                                         <  1:7850467‑M1/68‑1 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                      >  1:28506373/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                      >  1:14903511/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                      >  1:17696090/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                      >  1:27509303/1‑68 (MQ=255)
          aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                      >  1:2221166/1‑68 (MQ=255)
           ctcaCTAATCCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCgg                                                     <  1:5191368/65‑1 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCCCCCACCGCGTTGGGCCggg                                                    <  1:28907725/68‑1 (MQ=255)
            tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                    >  1:29725817/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                 >  1:1379358/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                 >  1:6802828/1‑68 (MQ=255)
               cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                 >  1:1921128/1‑68 (MQ=255)
                 aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGa                                               >  1:20595715/1‑68 (MQ=255)
                   tACCAGTCAGGAATCCTGGCTTATAAACAGAAAAATCAACCGCCCACCGCGGTGGGCCGGGGGGGAAg                                             <  1:25648194/68‑1 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTGTTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:20020746/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTGACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:3588483/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:26072790/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:937322/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:22531519/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:2340176/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:2431789/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:24914990/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:8749984/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:8258296/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:25842141/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:26312919/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:5877215/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:27925713/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:5681701/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:28749726/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:4963563/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:29278730/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:10321223/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:4741164/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:4609796/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:21682684/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:10044616/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:10797167/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:13561786/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:15475762/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:16124154/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:16607491/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:18875505/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGc                                            >  1:18971109/1‑68 (MQ=255)
                    aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGACCAACGCGTTGGGCCGGGGCGGAAGc                                            >  1:5486471/1‑68 (MQ=255)
                    aCCAGTCAGCAATCCTGGCTTATTAACAGTAAAAACAACCGCCCAACGCGTTGGGCCGGGGGGGAAga                                            >  1:28375137/1‑67 (MQ=255)
                     ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                           >  1:3050408/1‑68 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                           >  1:19324749/1‑67 (MQ=255)
                      cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCg                                           >  1:23753785/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                        >  1:20216146/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                        >  1:19645010/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                        >  1:9137630/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                        >  1:10549045/1‑67 (MQ=255)
                         tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                        >  1:11023303/1‑67 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCggttgattttactgtta                                       <  1:12249585‑M1/67‑18 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGa                                       <  1:24789522/67‑1 (MQ=255)
                          cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCggttgattttactgt                                         <  1:22915460‑M1/65‑16 (MQ=255)
                              aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTg                                  <  1:21968876/68‑1 (MQ=255)
                              aaTCCTGGCTTATAAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGt                                   <  1:14229410/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  >  1:11471830/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  <  1:14051699/67‑1 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  >  1:442615/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  >  1:2733725/1‑67 (MQ=255)
                               aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  >  1:17405210/1‑67 (MQ=255)
                               aTCCTGGATTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                  <  1:21057329/67‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCACCCGCCCCACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                            <  1:18215517/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                            <  1:11291869/68‑1 (MQ=255)
                                    ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                            <  1:8036891/68‑1 (MQ=255)
                                     gCTTATAAACATTAAAATCAACCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                            <  1:27702042/67‑1 (MQ=255)
                                     ccTTATTAACAGTAAAAACAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCa                            <  1:896284/66‑1 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                          >  1:14508336/1‑68 (MQ=255)
                                      cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                          >  1:17301546/1‑68 (MQ=255)
                                      cTTATTAACAGTAAAATCAAACGCCAAACGCGTTGGGCCGGGGGGGAAGCGTCGAAAGTTCCAGAAca                          >  1:24005201/1‑67 (MQ=255)
                                           taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                       >  1:4096608/1‑66 (MQ=255)
                                             aCAGTAAAATCAACCGCCCACCGGGTTGGGCggttgattttacggttaataagccaggattcctg                      >  1:21335249‑M1/1‑31 (MQ=255)
                                                     ataaaCCGCCCAACGCGTTGGGCCGGGGGGGAAGCCTCGACAGTGCCAGCACCATTACTGGCCGCAcc           >  1:7183082/4‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc           >  1:16420100/1‑68 (MQ=255)
                                                     aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc           <  1:13057050/68‑1 (MQ=255)
                                                       cacccGCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc         <  1:24464482/65‑1 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg          >  1:12404372/1‑67 (MQ=255)
                                                       cAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg          >  1:7450153/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattaggcag       >  1:26731988‑M1/1‑19 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc        >  1:18517151/1‑67 (MQ=255)
                                                         aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc        >  1:20796712/1‑67 (MQ=255)
                                                           cGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt     >  1:7343516/1‑68 (MQ=255)
                                                            gCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt    <  1:23259326/68‑1 (MQ=255)
                                                            gCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt    >  1:12673780/1‑68 (MQ=255)
                                                             cccAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattagtcagtaaa   >  1:18029134‑M1/1‑15 (MQ=255)
                                                             cccAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattagtcagtaaa   >  1:16836773‑M1/1‑15 (MQ=255)
                                                              ccAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattagtcagtaaaa  >  1:6100813‑M1/1‑14 (MQ=255)
                                                              ccAACGCGTTGGGCggttgattttactgttaataagccaggattcctgactggtattagtcagtaaaa  >  1:26980563‑M1/1‑14 (MQ=255)
                                                                 |                                                                
TATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAT  >  CP009273/999533‑999662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: