| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description | 
| RA | CP000730 | 2,679,462 | (A)8→7 | intergenic (+334/‑55) | USA300HOU_2535 → / → cbpD | hypothetical protein/possible secretory antigen | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,679,455 | 0 | A | . | 96.9% | 140.6 / ‑1.5 | 32 | intergenic (+327/‑62) | USA300HOU_2535/cbpD | hypothetical protein/possible secretory antigen | 
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (12/19); total (12/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 | |||||||||||
TGACATCCTTGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCA  >  CP000730/2679308‑2679583                                                                                                                                                   |                                                                                                                                tGACATCCTTGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTa                                                                                                                                  >  2:335978/1‑148 (MQ=255)  aCATCCTTGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTaaa                                                                                                                                >  2:326692/1‑148 (MQ=255)  aCATCCTTGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTaaa                                                                                                                                >  2:326705/1‑148 (MQ=255)        ttGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTaaaaaaatt                                                                                                                          >  1:59593/1‑146 (MQ=255)         tGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTaaaaaaattt                                                                                                                         <  2:271971/148‑4 (MQ=255)            aTCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTTTTTTTTTACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGtt                                                                                                                     >  2:386193/1‑148 (MQ=255)                 tATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTATTTTTTTTTTCTTTATATTATAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAATAATTTGTTATCAt                                                                                                                >  2:179893/1‑148 (MQ=255)                       aattaattGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTaa                                                                                                          >  1:128654/1‑148 (MQ=255)                               gTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAAt                                                                                                  >  2:71006/1‑148 (MQ=255)                                 tAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATaa                                                                                                <  2:62268/148‑1 (MQ=255)                                 tAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATaa                                                                                                >  1:375689/1‑148 (MQ=255)                                       cATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTCCTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGaaaa                                                                                          >  2:394421/1‑148 (MQ=255)                                          gAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGaaaaaaa                                                                                       <  1:273362/148‑1 (MQ=255)                                            aTTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAAtt                                                                                     >  1:326527/1‑148 (MQ=255)                                                  aCATTTTATGAAAATATCCATTTTTATTACAATTAAACAATTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACATTAAAAAAAAAAAAAAATTAACTTg                                                                               <  2:188415/148‑1 (MQ=255)                                                  aCATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAGATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTg                                                                               <  1:349806/148‑1 (MQ=255)                                                  aCATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTg                                                                               <  2:45523/148‑1 (MQ=255)                                                                aaTCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:130017/148‑1 (MQ=255)                                                                aaTCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:1961/148‑1 (MQ=255)                                                                aaTCCATTTTTATTACAAGTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  1:335509/148‑1 (MQ=255)                                                                aaTCCATTTTTAATACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCTTAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:58198/148‑1 (MQ=255)                                                                   ccATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAat                                                              <  2:95117/148‑1 (MQ=255)                                                                   ccATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTAT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAat                                                              <  2:274905/148‑1 (MQ=255)                                                                         ttattaCAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATaaaa                                                        <  1:13181/148‑1 (MQ=255)                                                                              aCAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATa                                                   <  1:324668/148‑1 (MQ=255)                                                                                         cTTTATATTACAAATTCATTACAGATACTACAAATATTGATTGCTTACAATGTATTCTAAAAAAAATTT‑TTATCATTATTAAAATAAAAAAAAAAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGttt                                        <  2:198157/148‑1 (MQ=255)                                                                                           ttATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTttat                                      >  1:397706/1‑148 (MQ=255)                                                                                                tGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAAtt                                 <  1:286766/148‑1 (MQ=255)                                                                                                   cAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCt                              >  1:319856/1‑148 (MQ=255)                                                                                                            tACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTa                     <  1:261903/148‑1 (MQ=255)                                                                                                             aCAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTAc                    <  1:168446/148‑1 (MQ=255)                                                                                                                 ttACTTTTATTGTTGATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACGTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTCGCAGTACtttt                >  2:48860/1‑148 (MQ=255)                                                                                                                           tgttgATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTTATTTTCTCAGCAGAt      >  1:127797/1‑148 (MQ=255)                                                                                                                            gttgATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGAtt     >  1:5270/1‑148 (MQ=255)                                                                                                                              tgATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTc   <  2:170073/148‑1 (MQ=255)                                                                                                                              tgATTGCTTACATTAATTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTc   <  2:218950/148‑1 (MQ=255)                                                                                                                               gATTGCTTACATTGTTTTCT‑AAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCa  <  1:397394/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                TGACATCCTTGTATCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCA  >  CP000730/2679308‑2679583 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |