Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,787,510	T→C	intergenic (+56/‑184)	USA300HOU_2639 → / → USA300HOU_2640	hypothetical protein/BglG family transcriptional antiterminator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,787,510	0	T	C	100.0%	77.4 / NA	23	intergenic (+56/‑184)	USA300HOU_2639/USA300HOU_2640	hypothetical protein/BglG family transcriptional antiterminator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/14);  total (9/14)

AATAAATCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGTGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACGTATGTATAAGCTCTCAT  >  CP000730/2787389‑2787626                                                                                                                          |                                                                                                                     aataaATCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCt                                                                                            <  2:56116/148‑1 (MQ=255)  ataaATCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCtt                                                                                           >  1:238162/1‑148 (MQ=255)     aaTCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTAct                                                                                        <  1:54011/148‑1 (MQ=255)       tCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTActgc                                                                                      >  2:74478/1‑148 (MQ=255)            taataCCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTc                                                                                 <  1:198410/148‑1 (MQ=255)            taataCCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTc                                                                                 <  1:323319/148‑1 (MQ=255)                aCCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTa                                                                             >  1:311093/1‑148 (MQ=255)                  cAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTAtt                                                                           >  1:250936/1‑148 (MQ=255)                            tCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTc                                                                 >  2:66004/1‑148 (MQ=255)                                               cacaATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACCGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAATATCata                                              <  2:354712/148‑1 (MQ=255)                                                          aaataTTCATCAATTGCCCGACACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGc                                   <  2:129759/148‑1 (MQ=255)                                                           aaaaTTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCt                                  >  2:46067/1‑148 (MQ=255)                                                           aaaaTTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCt                                  >  2:208620/1‑148 (MQ=255)                                                           aaaaTTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCt                                  <  2:270676/148‑1 (MQ=255)                                                            aaaTTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCtt                                 <  2:53186/148‑1 (MQ=255)                                                            aaaTTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCtt                                 <  1:339360/148‑1 (MQ=255)                                                                         gCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGAc                    <  1:172036/148‑1 (MQ=255)                                                                            tGATACGTCCACAAGTACTATTTGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACgta                 <  1:22997/148‑1 (MQ=255)                                                                            tGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACgta                 <  2:196605/148‑1 (MQ=255)                                                                              aTACGTCCACAAGTACCAGATGCATAAATTCAGAACAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACgtatg               <  2:175080/148‑1 (MQ=255)                                                                                  gTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACGTATGtata           >  1:184676/1‑148 (MQ=255)                                                                                  gTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACGTATGtata           >  1:273617/1‑148 (MQ=255)                                                                                           tACCAGATGCATAAATTCAGAAAAACAGCAGCGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACGTATGTATAAGCTCTCAt  <  1:319508/148‑1 (MQ=255)                                                                                                                          |                                                                                                                     AATAAATCAAATAATACCAATACGCCTTCTAAACAAACACCTCAAGCACAATCACCTAAAAATTAATCATTTGCCTGATACGTCCACAAGTACCAGATGCATAAATTCAGAAAAACAGCAGTGAGATATTTTTCTATTTGAAATATCTTACTGCTGTTCTTTATTTGCATAATTCGCACTCTGTTAAAATCATATAAAAGTTCGCTTTAATGTCATTGACGTATGTATAAGCTCTCAT  >  CP000730/2787389‑2787626

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