Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,455,698 2 bp→AT coding (1216‑1217/1260 nt) USA300HOU_2315 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,455,6980CA93.9% 115.0 / 1.7 33coding (1216/1260 nt)USA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (4/27);  minor base . (0/2);  total (4/29)
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.65e-01
*CP0007302,455,6990AT93.9% 109.5 / ‑4.0 33H406L (CAT→CTT) USA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (4/27);  total (4/29)
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.96e-01

CAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAC  >  CP000730/2455558‑2455830
                                                                                                                                            ||                                                                                                                                   
cagcagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGaaaaa                                                                                                                           >  2:421083/1‑151 (MQ=255)
    agAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTg                                                                                                                       <  1:990367/151‑1 (MQ=255)
    agAATCAGTTAAACAAGCACCACAAAGTGACGAGTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTg                                                                                                                       <  1:673168/151‑1 (MQ=255)
      aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaa                                                                                                                     <  1:820552/151‑1 (MQ=255)
      aaTCAGTTAAACAAGCACCACAAAGTGACGAATTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCCACTTATAAAGAAGACAGACGTTAATCAGTTGTAATTCATGCA‑AAAAAATTATCGAAAAAACTgaa                                                                                                                     <  1:849368/151‑1 (MQ=255)
       aTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaag                                                                                                                    <  1:160313/151‑1 (MQ=255)
       aTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaag                                                                                                                    <  1:64913/151‑1 (MQ=255)
           gTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACa                                                                                                                <  1:1003148/151‑1 (MQ=255)
                 cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTc                                                                                                          <  1:31275/151‑1 (MQ=255)
                  aaGCACCACAAAGTGAAGAATAAACAAAAAATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAAAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAAAAAAAAAAAAAACTgaagattgttcttca                                                                                                         <  2:114793/151‑11 (MQ=255)
                               tGACGATTTAACAAATGATATAGATGAAGATGAAATAGAGAAATCATTAAAAGAACATGCGACTTATAGAGAAGACAGACGTTAATCAATTGCAATTGCTGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                            <  2:585515/151‑1 (MQ=255)
                               tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                            <  2:341975/151‑1 (MQ=255)
                                   gATTTAACAAAAGATTCAAATGAAGATAAAATAGAGCAATCATTAAAACAACCTGCGACTTACAAAGAAGACAGAAGTTCATCAGTTGTGATTGATG‑AAAAAAAAATAAAGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                        <  2:370307/151‑1 (MQ=255)
                                     tttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataatatc                                                                                      <  2:566494/151‑3 (MQ=255)
                                                  tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTa                                                                         <  1:945381/151‑1 (MQ=255)
                                                               aaaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaa                                                            <  1:691113/151‑1 (MQ=255)
                                                                aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                           <  1:427543/151‑1 (MQ=255)
                                                                aaTAGAGCAATCATTAAAAGAACCTGCGACTAATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                           <  1:712316/151‑1 (MQ=255)
                                                                 aTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg                                                          <  1:627038/151‑1 (MQ=255)
                                                                  tAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAga                                                         <  1:87828/151‑1 (MQ=255)
                                                                   agagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgat                                                        <  1:1017719/151‑1 (MQ=255)
                                                                      gCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatgat                                                     <  1:421083/151‑1 (MQ=255)
                                                                         aTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                                  <  1:54060/151‑1 (MQ=255)
                                                                                    aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                       <  2:634330/151‑1 (MQ=255)
                                                                                       cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                    <  1:532855/151‑1 (MQ=255)
                                                                                        tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                   <  2:986114/151‑1 (MQ=255)
                                                                                        tGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                   >  2:593968/1‑151 (MQ=255)
                                                                                         gCGACTTATAAAGAAAACAGACGTTGATCAGTTGTAATTGATGCA‑AAAAAATTAGCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                  <  2:572341/151‑1 (MQ=255)
                                                                                             cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                              <  1:375468/151‑1 (MQ=255)
                                                                                             cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                              <  1:745086/151‑1 (MQ=255)
                                                                                             cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                              <  2:166943/151‑1 (MQ=255)
                                                                                               tataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactat                            >  1:93391/1‑151 (MQ=255)
                                                                                                                         tGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAacac  >  2:273264/1‑151 (MQ=255)
                                                                                                                                            ||                                                                                                                                   
CAGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAC  >  CP000730/2455558‑2455830

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