Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	677,191	C→T	P32L (CCG→CTG)	USA300HOU_0620 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	677,191	0	C	T	100.0%	123.7 / NA	38	P32L (CCG→CTG)	USA300HOU_0620	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

TTAATTTTTATGATTTGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACCGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCTTT  >  CP000730/677046‑677329                                                                                                                                                  |                                                                                                                                           ttAATTTTTATGATTTGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGc                                                                                                                                          <  1:146688/148‑1 (MQ=255) ttAATTTTTATGATTTGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGc                                                                                                                                          <  1:3736/148‑1 (MQ=255)    aTTTTTATGATTTGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCAtt                                                                                                                                       >  2:117841/1‑148 (MQ=255)              tttGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATtgatga                                                                                                                             >  2:62324/1‑148 (MQ=255)                tGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGata                                                                                                                           <  1:65867/148‑1 (MQ=255)                         gAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAgaga                                                                                                                  <  2:195212/148‑1 (MQ=255)                         gAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAgaga                                                                                                                  <  2:49716/148‑1 (MQ=255)                            aTTAGACCAACACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGATGATTACTTTATGTATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGGGAAAg                                                                                                               <  2:106648/148‑1 (MQ=255)                             ttAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGt                                                                                                              <  1:109565/148‑1 (MQ=255)                              tAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTa                                                                                                             <  1:4937/148‑1 (MQ=255)                                  ccAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTa                                                                                                         >  2:114053/1‑148 (MQ=255)                                     atacaataAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTg                                                                                                      <  2:42611/148‑1 (MQ=255)                                      tacaataAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGa                                                                                                     >  2:100977/1‑148 (MQ=255)                                      tacaataAATAAAGATGTTTCATTTTAAAGGAGTACTACCAAAGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGa                                                                                                     >  2:299475/1‑148 (MQ=255)                                         aataaataAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGatat                                                                                                  <  2:160916/148‑1 (MQ=255)                                                 aaGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAAtt                                                                                          <  2:175459/148‑1 (MQ=255)                                                   gATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCt                                                                                        >  1:261563/1‑148 (MQ=255)                                                   gATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCt                                                                                        >  1:98300/1‑148 (MQ=255)                                                          cATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAaaaga                                                                                 <  1:209643/148‑1 (MQ=255)                                                          cATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAaaaga                                                                                 >  1:7970/1‑148 (MQ=255)                                                                  aGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCaa                                                                         >  1:266543/1‑148 (MQ=255)                                                                               aTGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAat                                                            <  2:266133/148‑1 (MQ=255)                                                                               aTGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAat                                                            <  2:206029/148‑1 (MQ=255)                                                                                 gCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTTCTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAatat                                                          <  2:264310/148‑1 (MQ=255)                                                                                             tAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTATGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAAtct                                              >  2:192697/1‑148 (MQ=255)                                                                                                agattaaCTAAAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCa                                           <  2:166425/142‑1 (MQ=255)                                                                                                    ttACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGcac                                       >  1:248148/1‑148 (MQ=255)                                                                                                                   gaagaTTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGa                        <  1:126753/148‑1 (MQ=255)                                                                                                                    aagaTTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGAt                       >  1:246142/1‑148 (MQ=255)                                                                                                                    aagaTTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGAt                       >  1:240234/1‑148 (MQ=255)                                                                                                                     agaTTACTTTAATCATAAAGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGAtt                      >  1:287469/1‑148 (MQ=255)                                                                                                                                     aaaGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAAtttgtttg      >  1:77896/1‑148 (MQ=255)                                                                                                                                     aaaGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAAtttgtttg      >  2:278045/1‑148 (MQ=255)                                                                                                                                     aaaGGTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAAtttgtttg      >  1:57785/1‑148 (MQ=255)                                                                                                                                        ggTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCtt   >  2:113680/1‑148 (MQ=255)                                                                                                                                        ggTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACCGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCtt   >  2:264121/1‑148 (MQ=255)                                                                                                                                         gTTTATCACTGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCttt  <  1:107859/148‑1 (MQ=255)                                                                                                                                         gTTTATCACTACATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCttt  <  2:128251/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           TTAATTTTTATGATTTGCGAAGGGGAAATTAAACCAATACAATAAATAAAGATGTTTCATTTTAAAGGAGTACTACCAATGCGTATTAACGATAAAATTTTACTAGAAAATATAGAAGATTACTTTAATCATAAAGGTTTATCACCGCATTTGATTGATGATATTAAAGAGAAAGTAATTACTGATATAAAAAATTCTGAAAAGAAAGATCAAGATTATATTGAATATAAAAGAAAATCTCCAGCACAAATCATATTAATGATTCAAAGGAATTTGTTTGCTTT  >  CP000730/677046‑677329

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