Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,017,678	T→C	F391S (TTC→TCC)	USA300HOU_0971 →	AGCS family alanine or glycine:sodium (Na+) or proton (H+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,017,678	0	T	C	100.0%	100.7 / NA	28	F391S (TTC→TCC)	USA300HOU_0971	AGCS family alanine or glycine:sodium (Na+) or proton (H+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/13);  total (15/13)

AAATGCAACAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTTCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGATTATTTTGTTCGCTACA  >  CP000730/1017538‑1017811                                                                                                                                             |                                                                                                                                      aaaTGCAACAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCa                                                                                                                                >  1:112717/1‑148 (MQ=255) aaaTGCAACAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCa                                                                                                                                >  1:78962/1‑148 (MQ=255)       aaCAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGTTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTtact                                                                                                                          <  1:116950/148‑1 (MQ=255)       aaCAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTtact                                                                                                                          <  1:371404/148‑1 (MQ=255)                  attattCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAGGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTc                                                                                                               >  2:368441/1‑148 (MQ=255)                    tattCAGTTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGt                                                                                                             <  2:82216/148‑1 (MQ=255)                    tattCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGt                                                                                                             <  2:43606/148‑1 (MQ=255)                       tCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGtact                                                                                                          <  2:326974/148‑1 (MQ=255)                          ggTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGtactact                                                                                                       >  2:194168/1‑148 (MQ=255)                                               gCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTg                                                                                  <  1:338774/148‑1 (MQ=255)                                                   gCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTa                                                                              <  1:235347/148‑1 (MQ=255)                                                            aaGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAAcgcg                                                                     >  2:106475/1‑148 (MQ=255)                                                                 ttCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATc                                                                <  1:194168/148‑1 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAATTGTTGTTTTTTTAACGCGTAATCaa                                                              >  1:109237/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  1:126622/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  2:76004/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  2:52174/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  2:338805/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  2:355837/1‑148 (MQ=255)                                                                   ccATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCaa                                                              >  2:230222/1‑148 (MQ=255)                                                                      tGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAaat                                                           >  2:403774/1‑148 (MQ=255)                                                                       ggCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAaata                                                          <  1:181242/148‑1 (MQ=255)                                                                                           tCCTTCTNTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACTGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTATCATCGAtattt                                      <  2:312034/148‑1 (MQ=255)                                                                                             cTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGAtatttat                                    >  2:196631/1‑148 (MQ=255)                                                                                                              ggTTCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGa                   <  1:84679/148‑1 (MQ=255)                                                                                                                 tCGTACTTTATTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGatta                <  2:57123/148‑1 (MQ=255)                                                                                                                 tCGTACTTTAGTGCATTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGatta                <  2:50240/148‑1 (MQ=255)                                                                                                                               aTTTGCTTTATTCTCCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGATTATTTTGTTCGCTACa  >  2:75169/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                      AAATGCAACAGGTAAAGATTATTCAGGTACTGCGATGTATGCACAAGCCGGCATTGATAAAGCGTTCCATGGCAGTGGTTATCAATTTGATCCTACTTTCTCTGGCGTAGGTTCGTACTTTATTGCATTTGCTTTATTCTTCTTTGCATTTACTACAATTTTGTCGTACTACTACATTACAGAAACAAATGTTGCTTATTTAACGCGTAATCAAAATAATCAAGTTTCATCGATATTTATTAATATTGCTCGTGTGATTATTTTGTTCGCTACA  >  CP000730/1017538‑1017811

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