Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,388,866	A→G	100%	G8G (GGA→GGG)	USA300HOU_2249 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,388,866	0	A	G	100.0%	80.1 / NA	25	G8G (GGA→GGG)	USA300HOU_2249	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/14);  total (11/14)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATGTTTCTTCTTAGATGTTTGTAATTCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGACAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTAGTGATTCTTAAGAGTGGGATAGAAATGAT  >  CP000730/2388732‑2389003                                                                                                                                       |                                                                                                                                          tATGTTTCTTCTTAGATGTTTGTAATTCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACtt                                                                                                                                     >  1:141101/1‑141 (MQ=255)             tAGATGTTTGTAATTCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCg                                                                                                                         <  1:29268/141‑1 (MQ=255)                      gTAATTCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACAttt                                                                                                                >  1:275176/1‑141 (MQ=255)                           tCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATAcc                                                                                                           >  1:39345/1‑141 (MQ=255)                                   ttatACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAg                                                                                                   <  2:39345/141‑1 (MQ=255)                                                 aaaTTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATGTATTACAACta                                                                                     >  2:79399/1‑140 (MQ=255)                                                  aaTTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCa                                                                                    >  2:213006/1‑141 (MQ=255)                                                  aaTTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCa                                                                                    >  1:163419/1‑141 (MQ=255)                                                    tttATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAAt                                                                                  <  2:141101/141‑1 (MQ=255)                                                       attatATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGc                                                                               >  2:5074/1‑141 (MQ=255)                                                        ttatATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCa                                                                              >  2:41455/1‑141 (MQ=255)                                                           tatatCCTACTTAGCGTGTGATCTATTCCAAATATTATAGTGGTGGCTATTTCTGGTTGAGATTTATATATTTTGGGCAACTTGGTGTTTATTCGACAACATTTATACCTAAAAAAGATAGATTACAACTCAATGTCATgg                                                                           <  2:167912/141‑1 (MQ=255)                                                           tatatCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATTGCATgg                                                                           >  1:145071/1‑141 (MQ=255)                                                                   aCTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGAc                                                                   <  2:24022/141‑1 (MQ=255)                                                                   aCTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGAc                                                                   <  2:163419/141‑1 (MQ=255)                                                                        gCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCaa                                                              <  2:79276/141‑1 (MQ=255)                                                                             ttATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTa                                                         <  1:243394/141‑1 (MQ=255)                                                                             gtaggatatACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTa                                                         <  2:97404/133‑1 (MQ=255)                                                                              tATCGATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTaa                                                        <  1:25730/141‑1 (MQ=255)                                                                                        aacttttgaaTTGCTGGCTACTTCGGGTGTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTa                                              <  1:134671/133‑1 (MQ=255)                                                                                                       ggCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTa                               <  1:75444/141‑1 (MQ=255)                                                                                                       ggCTATATCTTGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTa                               <  2:107751/141‑1 (MQ=255)                                                                                                                ttGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCAATAAAGATTTGTTATGTAATGATTCtt                      >  1:1473/1‑141 (MQ=255)                                                                                                                 tGTTTAGATTTATATATTTTGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTAGTGATTCTTa                     >  2:61203/1‑141 (MQ=255)                                                                                                                                    tGGGCAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTAGTGATTCTTAAGAGTGGGATAGAAATGat  <  1:224048/141‑1 (MQ=255)                                                                                                                                       |                                                                                                                                          TATGTTTCTTCTTAGATGTTTGTAATTCATATTATTATACAATTACTTAAATTTATTATATATCAAACTTAGCATTTTATCTATTACAAATATTAAATTGCTGGCTATATCTTGTTTAGATTTATATATTTTGGACAACTTGTTGTTTATTCGACAACATTTATACCTAAAAAAGATATATTACAACTCAATGGCATGGAGTGTGACATCAATGTTAAAGAACCACTAAAGATTTATTATGTAGTGATTCTTAAGAGTGGGATAGAAATGAT  >  CP000730/2388732‑2389003

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