Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 108,148 A→G 11.8% K1037R (AAA→AGA)  USA300HOU_0106 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730108,1480AG11.8% 86.9 / 4.7 34K1037R (AAA→AGA) USA300HOU_0106hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/21);  new base G (2/2);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG  >  CP000730/108009‑108288
                                                                                                                                           |                                                                                                                                            
tATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGaaa                                                                                                                                             <  2:154309/141‑1 (MQ=255)
  ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc                                                                                                                                           <  1:56091/141‑1 (MQ=255)
  ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc                                                                                                                                           <  1:443731/141‑1 (MQ=255)
  ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTGGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATATAGATGAAAcc                                                                                                                                           <  1:60695/141‑1 (MQ=255)
      gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt                                                                                                                                       >  1:432536/1‑141 (MQ=255)
      gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt                                                                                                                                       >  1:293801/1‑141 (MQ=255)
                    taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt                                                                                                                         >  1:71913/1‑141 (MQ=255)
                    taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt                                                                                                                         >  1:302309/1‑141 (MQ=255)
                        ccccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:359129/140‑1 (MQ=255)
                        aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:84766/141‑1 (MQ=255)
                        aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa                                                                                                                     <  2:269728/141‑1 (MQ=255)
                          ccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaaga                                                                                                                   >  2:189597/1‑141 (MQ=255)
                             aGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaagaaag                                                                                                                >  1:227057/1‑141 (MQ=255)
                                tGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAgaaagaaa                                                                                                             <  1:143909/141‑1 (MQ=255)
                                 ggTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAAt                                                                                                            <  1:442946/141‑1 (MQ=255)
                                  gTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATg                                                                                                           <  1:216517/141‑1 (MQ=255)
                                  gTGCTATGAAAAGGTCATGAGAAAAACCCAAATTGTTAAAAGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAACGAATACAGCTGAAACCATCTTATGAGACGATTTTTAAAGCAAGAAATg                                                                                                           <  2:298691/141‑1 (MQ=255)
                                               gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc                                                                                              <  1:460406/141‑1 (MQ=255)
                                               gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc                                                                                              <  2:372951/141‑1 (MQ=255)
                                                            ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGtatata                                                                                 <  2:432536/141‑1 (MQ=255)
                                                            ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAAAAAAAAAAAAAGtatata                                                                                 <  2:293801/141‑1 (MQ=255)
                                                                aCTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAgg                                                                             >  1:221888/1‑141 (MQ=255)
                                                                 cTTGTTAAACGTTGCAGTGACAAGAGGAAAGAGAGAGGTGTGTGGGATGGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGa                                                                            >  2:20480/1‑141 (MQ=255)
                                                                            ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACAtt                                                                 <  2:364321/141‑1 (MQ=255)
                                                                            ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAACAAAAAAAAAAGTATATAGAGGAAGTTATACAtt                                                                 <  2:418026/141‑1 (MQ=255)
                                                                             tGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACAttt                                                                <  1:98540/141‑1 (MQ=255)
                                                                                   gACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAg                                                          <  2:221888/141‑1 (MQ=255)
                                                                                    aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg                                                         <  2:317835/141‑1 (MQ=255)
                                                                                    aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg                                                         <  2:348154/141‑1 (MQ=255)
                                                                                                                 tGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTaac                            <  1:298101/141‑1 (MQ=255)
                                                                                                                     gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt                        <  2:302309/141‑1 (MQ=255)
                                                                                                                     gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt                        <  2:71913/141‑1 (MQ=255)
                                                                                                                           cAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGatat                  >  1:161637/1‑141 (MQ=255)
                                                                                                                             aaGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTa                >  2:178362/1‑141 (MQ=255)
                                                                                                                                   aCAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCa          >  1:378502/1‑141 (MQ=255)
                                                                                                                                           aaCCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTgg  <  1:177476/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG  >  CP000730/108009‑108288

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