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
*CP0007302,679,4550A.96.9% 140.6 / ‑1.5 32intergenic (+327/‑62)USA300HOU_2535/cbpDhypothetical 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)
            aTCACTATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTTTTTTTTTACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGtt                                                                                                                     >  2:386193/1‑148 (MQ=255)
                 tATGACAATTAATTGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTATTTTTTTTTTCTTTATATTATAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAATAATTTGTTATCAt                                                                                                                >  2:179893/1‑148 (MQ=255)
                       aattaattGTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTaa                                                                                                          >  1:128654/1‑148 (MQ=255)
                               gTTAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAAt                                                                                                  >  2:71006/1‑148 (MQ=255)
                                 tAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATaa                                                                                                <  2:62268/148‑1 (MQ=255)
                                 tAAATACATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATaa                                                                                                >  1:375689/1‑148 (MQ=255)
                                       cATGAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTCCTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGaaaa                                                                                          >  2:394421/1‑148 (MQ=255)
                                          gAATTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGaaaaaaa                                                                                       <  1:273362/148‑1 (MQ=255)
                                            aTTTCTACATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAAtt                                                                                     >  1:326527/1‑148 (MQ=255)
                                                  aCATTTTATGAAAATATCCATTTTTATTACAATTAAACAATTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACATTAAAAAAAAAAAAAAATTAACTTg                                                                               <  2:188415/148‑1 (MQ=255)
                                                  aCATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAGATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTg                                                                               <  1:349806/148‑1 (MQ=255)
                                                  aCATTTTATGAAAAAATCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTg                                                                               <  2:45523/148‑1 (MQ=255)
                                                                aaTCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:130017/148‑1 (MQ=255)
                                                                aaTCCATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:1961/148‑1 (MQ=255)
                                                                aaTCCATTTTTATTACAAGTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  1:335509/148‑1 (MQ=255)
                                                                aaTCCATTTTTAATACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCTTAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATgg                                                                 <  2:58198/148‑1 (MQ=255)
                                                                   ccATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAat                                                              <  2:95117/148‑1 (MQ=255)
                                                                   ccATTTTTATTACAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTATAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAat                                                              <  2:274905/148‑1 (MQ=255)
                                                                         ttattaCAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATaaaa                                                        <  1:13181/148‑1 (MQ=255)
                                                                              aCAATTCAACACTTTATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATa                                                   <  1:324668/148‑1 (MQ=255)
                                                                                         cTTTATATTACAAATTCATTACAGATACTACAAATATTGATTGCTTACAATGTATTCTAAAAAAAATTTTTATCATTATTAAAATAAAAAAAAAAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGttt                                        <  2:198157/148‑1 (MQ=255)
                                                                                           ttATATGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTttat                                      >  1:397706/1‑148 (MQ=255)
                                                                                                tGACAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAAtt                                 <  1:286766/148‑1 (MQ=255)
                                                                                                   cAACTTCATTACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCt                              >  1:319856/1‑148 (MQ=255)
                                                                                                            tACAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTa                     <  1:261903/148‑1 (MQ=255)
                                                                                                             aCAGTTACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTAc                    <  1:168446/148‑1 (MQ=255)
                                                                                                                 ttACTTTTATTGTTGATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACGTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTCGCAGTACtttt                >  2:48860/1‑148 (MQ=255)
                                                                                                                           tgttgATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTTATTTTCTCAGCAGAt      >  1:127797/1‑148 (MQ=255)
                                                                                                                            gttgATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGAtt     >  1:5270/1‑148 (MQ=255)
                                                                                                                              tgATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTc   <  2:170073/148‑1 (MQ=255)
                                                                                                                              tgATTGCTTACATTAATTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTc   <  2:218950/148‑1 (MQ=255)
                                                                                                                               gATTGCTTACATTGTTTTCTAAAAAAATTTGTTATCATAATTAACGTTGAATAAAGAAAAAAATTAACTTGGGAGATAAAAATGGAATATAAAAAGATACTAATTCGTTTATTAATTGCTTTTGCAGTACTTTTCTCAGCAGATTTCa  <  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: