Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	479,553	G→A	A194T (GCA→ACA)	USA300HOU_0458 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	479,553	0	G	A	100.0%	85.8 / NA	26	A194T (GCA→ACA)	USA300HOU_0458	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

AGATGAACAAGGTGAGCCACTTAGCGAACAAGTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTCATTTTTCAATAGATCAGGAA  >  CP000730/479429‑479681                                                                                                                             |                                                                                                                                 aGATGAACAAGGTGAGCCACTTAGCGAACAAGTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAgc                                                                                                                  >  1:277543/1‑141 (MQ=255)         aaGGTGAGCCACTTAGCGAACAAGTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATg                                                                                                          <  2:264099/141‑1 (MQ=255)          aGGTGAGCCACTTAGCGAACAAGTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGa                                                                                                         <  2:187862/141‑1 (MQ=255)                       gttcgctaaGAGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAAACTCTATCGAGTTGGACAATGCCGAAGCGTGCGCAAAGTTCTTACCATACATGGTTGCATTTATCGCAACATGACCATAGATTTatt                                                                                            >  2:220355/8‑139 (MQ=255)                                gTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAg                                                                                   >  1:56313/1‑141 (MQ=255)                                  gAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAgcg                                                                                 <  2:243508/141‑1 (MQ=255)                                    aTCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGcc                                                                               >  2:162465/1‑141 (MQ=255)                                         cAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTg                                                                          >  2:199601/1‑141 (MQ=255)                                          aaaTGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATATTTTTACTAAAGCACAGCGCCAAGTGa                                                                         >  2:250207/1‑141 (MQ=255)                                              gATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTaa                                                                     >  2:128018/1‑141 (MQ=255)                                                      ggACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTAc                                                             <  2:55373/141‑1 (MQ=255)                                                         cGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAGGCACAGGGGCAAGTGATTAAAGGGTTACCCa                                                          >  2:155573/1‑141 (MQ=255)                                                           aaaCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAACGTCTTTAACATACCTGGTCGCATTCAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAAt                                                        <  1:59003/141‑1 (MQ=255)                                                            aaCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATg                                                       <  2:277543/141‑1 (MQ=255)                                                             aCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGa                                                      <  1:17008/141‑1 (MQ=255)                                                             aCTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGa                                                      >  1:247057/1‑141 (MQ=255)                                                                ttATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGGGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTa                                                                                           <  2:298770/101‑1 (MQ=255)                                                                ttATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTa                                                                                           >  1:298770/1‑101 (MQ=255)                                                                   tatCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGa                                                >  2:278334/1‑141 (MQ=255)                                                                          cAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACg                                         <  1:205086/141‑1 (MQ=255)                                                                            aTTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGAt                                       <  1:117622/141‑1 (MQ=255)                                                                             tttCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATa                                      >  1:34013/1‑141 (MQ=255)                                                                                             cAATGCCGAAGCGTGAGCAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTc                      <  1:56856/141‑1 (MQ=255)                                                                                                              cAAAGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTCATTTTTCAATAGATCAg     >  2:258272/1‑141 (MQ=255)                                                                                                               aaaGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTCATTTTTCAATAGATCAgg    <  2:137739/141‑1 (MQ=255)                                                                                                                 aGTTTTTACCATACATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTCATTTTTCAATAGATCAGGaa  <  1:126603/141‑1 (MQ=255)                                                                                                                             |                                                                                                                                 AGATGAACAAGGTGAGCCACTTAGCGAACAAGTGAATCGTCAAATGATTAAATGGACGAAACTTTATATCGATCAATTTCTATCGAGTTGGACAATGCCGAAGCGTGAGCAAAGTTTTTACCATGCATGGTTGCATTTAGCGCAACATGACCATAGTTTTACTAAAGCACAGCGCCAAGTGATTAAAGGCTTACCCAATGATCCTGAAATGACGATAGAGTCAGTATTAACTCATTTTTCAATAGATCAGGAA  >  CP000730/479429‑479681

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