Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,930,186 A→G 12.1% *34* (TAA→TAG)  USA300HOU_1796 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,930,1860AG12.1% 95.8 / 5.4 43*34* (TAA→TAG) USA300HOU_1796hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/13);  new base G (3/2);  total (27/16)
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.02e-01

TTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGACCATTTGAT  >  CP000730/1930046‑1930326
                                                                                                                                            |                                                                                                                                            
ttGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATAAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTaa                                                                                                                                              >  1:259313/1‑141 (MQ=255)
       aaTTATACTGCTAATTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCtt                                                                                                                                       <  1:270462/141‑1 (MQ=255)
          tataCTGCTAATTGTTCTGTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTATCTTCCttttt                                                                                                                                    <  2:82193/141‑1 (MQ=255)
           ataCTGCTAATTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTATATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTa                                                                                                                                   <  1:388343/141‑1 (MQ=255)
                   aaTTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTgg                                                                                                                           <  2:124205/141‑1 (MQ=255)
                     tttttCTTTTGGTAAGATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAg                                                                                                                         <  1:21401/141‑1 (MQ=255)
                          cTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTaa                                                                                                                    <  2:349916/141‑1 (MQ=255)
                             tttGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGt                                                                                                                 >  2:192820/1‑141 (MQ=255)
                                    aaTATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAAc                                                                                                          >  1:176827/1‑141 (MQ=255)
                                    aaTATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAAc                                                                                                          >  1:236643/1‑141 (MQ=255)
                                     atatGCAAGGTTTTAAAGAGAAACATCAAGAACTAATAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATGGCTTAACTTCCTTTTTACTGTGTGGAGCGTAGAGTTTTGAACa                                                                                                         >  2:77738/1‑141 (MQ=255)
                                     atatGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACa                                                                                                         >  2:192061/1‑141 (MQ=255)
                                      tatGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACat                                                                                                        >  2:215356/1‑141 (MQ=255)
                                          cAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGCCAAATTGGCCTGATGAGTTCAATATCGGAAGTTAGGCAACTAAACATTGCTTAGCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACataata                                                                                                    >  1:304743/1‑141 (MQ=255)
                                             ggTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACataataata                                                                                                 >  2:184959/1‑141 (MQ=255)
                                                     agagAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTg                                                                                         <  2:69809/141‑1 (MQ=255)
                                                       agaAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTgcg                                                                                       <  2:174115/141‑1 (MQ=255)
                                                              tCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGa                                                                                >  1:25628/1‑141 (MQ=255)
                                                              tCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGa                                                                                <  2:6112/141‑1 (MQ=255)
                                                               cAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGAt                                                                               >  2:325072/1‑141 (MQ=255)
                                                                    aCTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTCGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTaa                                                                          <  1:350011/141‑1 (MQ=255)
                                                                     cTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAAc                                                                         >  2:55900/1‑141 (MQ=255)
                                                                          aaaGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCa                                                                    >  2:453439/1‑141 (MQ=255)
                                                                          aaaGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCa                                                                    >  2:416602/1‑141 (MQ=255)
                                                                             ggCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATaa                                                                 >  1:86731/1‑141 (MQ=255)
                                                                               ctctATGCCAAATTGGCCTGATGAGTTCAATATCGGAAGTTAGGCAACTAAACATTGCTTAGCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTTCGCAAATGATTATAACTTcctta                                                                  <  2:418833/138‑4 (MQ=25)
                                                                               ctctATGCCAAATTGGCCTGATGAGTTCAATATCGGAAGTTAGGCAACTAAACATTGCTTAGCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTTCGCAAATGATTATAACTTCCTTa                                                                  >  1:418833/1‑138 (MQ=255)
                                                                                             ggACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGtt                                                 <  1:97604/141‑1 (MQ=255)
                                                                                                            atatCCGAAGTTAAGCAACTAAACGTTGCTTAACTTCCTTTTTACTTTTTGTAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGTTGTACGTTTAATTAattttattt                                  >  2:356464/1‑141 (MQ=255)
                                                                                                                ccGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGtt                              <  2:89982/141‑1 (MQ=255)
                                                                                                                ccGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGtt                              >  1:55930/1‑141 (MQ=255)
                                                                                                                   aaGTTAGTCAGCTTATCATTGCTTAGCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTTCGCAAATTTTTATAACTTCCTTCACCAAAACATGTACGTTTAATTAATTTTATTTTgttgtt                           <  1:280680/141‑1 (MQ=255)
                                                                                                                         aGCAACTAAACATTGCTTAACTTCCTTTTTTTTTTTTTTTTCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTTTTATTTATAcc                     >  2:254148/1‑141 (MQ=255)
                                                                                                                         aGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTTTTTTTTTTTTTTTTGAACATAATAATATTCGATTTCTCAAATGATTTTAACTTCCATAACCAAAAGATTTACGTTTAATTATTTTTATTTTATTATTTATAcc                     >  2:410240/1‑141 (MQ=255)
                                                                                                                         aGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATAcc                     >  1:336022/1‑141 (MQ=255)
                                                                                                                          gCAACTAAACATTGTATAACTTTCTTTTTTTTTTTTTTTTTTTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAATCAAATGATGTACGTTTAATTAATTTTATTATGttagttattcat                    >  2:300871/1‑137 (MQ=255)
                                                                                                                           caatgtttagttGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCtt                   >  2:350737/11‑141 (MQ=255)
                                                                                                                              tgtttagttGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGGACATAATAATAATCGATTGTGCACATGATTGTAACTTCCATAACCAAAAGATGTACGTTAAATTAATTTTATTTTGTTTTTTATAtcttctt                >  2:141169/8‑140 (MQ=255)
                                                                                                                               tAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTaa               <  1:56038/141‑1 (MQ=255)
                                                                                                                                 ttagttGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAg             >  1:41695/5‑141 (MQ=255)
                                                                                                                                   cATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGa           >  2:123049/1‑141 (MQ=255)
                                                                                                                                   agttGCTTAACTTCCTTTTTACTTTTTGGTGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTTTAACTTCCATAGGa           >  2:253849/3‑141 (MQ=255)
                                                                                                                                      tGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGACCa        <  1:324363/141‑1 (MQ=255)
                                                                                                                                       gCTTAGCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTTCGCAAATGATTATAACTTCCTTAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTGTTTATCCCTTCTAAAGGACCAt       >  1:201900/1‑141 (MQ=255)
                                                                                                                                           aaCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGACCATTTGa   <  1:369738/141‑1 (MQ=255)
                                                                                                                                            aCTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGACCATTTGAt  <  1:76862/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATATGCAAGGTTTTAAAGAGAAACATCAAGAACTAAAAAAGGCTCTATGTCAAATTGGACTGATGCGTTCAATATCCGAAGTTAAGCAACTAAACATTGCTTAACTTCCTTTTTACTTTTTGGAGCGTAAAGTTTTGAACATAATAATATTCGATTGCGCAAATGATTGTAACTTCCATAACCAAAAGATGTACGTTTAATTAATTTTATTTTGTTATTTATACCTTCTAAAGGACCATTTGAT  >  CP000730/1930046‑1930326

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: ‑