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
*	CP000730	2,455,698	0	C	A	100.0%	43.9 / NA	13	H406I (CAT→ATT)	USA300HOU_2315	hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (2/11);  total (2/11)
*	CP000730	2,455,699	0	A	T	100.0%	41.9 / NA	13	H406I (CAT→ATT)	USA300HOU_2315	hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (2/11);  total (2/11)

AATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAA  >  CP000730/2455564‑2455827                                                                                                                                       ||                                                                                                                                 aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaa                                                                                                                  <  1:269862/151‑1 (MQ=255)              aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                     <  1:110165/151‑1 (MQ=255)                      aaaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTCCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAAc                                                                                             <  1:409884/151‑1 (MQ=255)                          tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAAAAAAAAAA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                         <  2:166478/151‑1 (MQ=255)                                            ttCAAATGAAGATGAAATAGAGCAATAATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCt                                                                       <  2:245765/151‑1 (MQ=255)                                                  tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACTTTAATCAGTTGTAAATGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                                                                 <  1:291220/151‑1 (MQ=255)                                                  tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                                                                 <  2:400333/151‑1 (MQ=255)                                                   gAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACct                                                                >  1:404491/1‑151 (MQ=255)                                                                    aTCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCt                                               <  1:174197/151‑1 (MQ=255)                                                                                  cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                 <  1:108857/151‑1 (MQ=255)                                                                                           ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactata                        <  1:313735/151‑1 (MQ=255)                                                                                                  agacagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATattatt                 <  2:404491/151‑1 (MQ=255)                                                                                                                 aGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTaa  >  2:233711/1‑151 (MQ=255)                                                                                                                                       ||                                                                                                                                 AATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAA  >  CP000730/2455564‑2455827

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