Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	469,786	A→G	100%	G118G (GGA→GGG)	USA300HOU_0445 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	469,786	0	A	G	100.0%	80.7 / NA	26	G118G (GGA→GGG)	USA300HOU_0445	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/5);  total (21/5)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACT  >  CP000730/469652‑469922                                                                                                                                       |                                                                                                                                         gaTGGAGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCatta                                                                                                                                     <  1:343965/140‑1 (MQ=255)  aTGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGAAATGAAGGTATGGTCCTATACTTAGAAAGAAATACACGGACTGCAAAGGGGCattat                                                                                                                                    >  2:182297/1‑140 (MQ=255)    ggtgAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTAttt                                                                                                                                  <  2:233568/137‑1 (MQ=255)       gAATTTAAAAAGGGTGACAAAAAGTGGTGGACGCTATATACAGATTTAGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGCCATATACTTAGATAGAAATACACGAACTGTAAAGGGGGATtttttttt                                                                                                                               >  2:164152/1‑138 (MQ=255)                 aGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTc                                                                                                                     <  1:269663/140‑1 (MQ=255)                  gggTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCt                                                                                                                    <  2:105900/140‑1 (MQ=255)                      gATAAAGGACTGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAGATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTataa                                                                                                                >  2:112434/1‑140 (MQ=255)                                  tgGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAAtt                                                                                                    >  1:72401/1‑140 (MQ=255)                                             aaCAGATTGCTCCAGAAGTAATAAGCAGGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCAAGATAGaa                                                                                         >  2:298859/2‑140 (MQ=255)                                                     tCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATtat                                                                                 >  2:238104/1‑140 (MQ=255)                                                      cGCCAAAAGTACTAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTGAGATAGACATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAAGGCCAATTCCCAGATAGaaaaacctcta                                                                                >  2:155520/1‑134 (MQ=255)                                                      cGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTTTGTAATAAGGGCAAATTCGCAGATTGAAAAAATtata                                                                                >  1:194134/1‑140 (MQ=255)                                                              gTAATAATCAAGGTGGATTGGGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTGGTTAAGACATTCTATAATAAGGGCAAATTCCCAGCTAGAAAAAATTATAAAGTtgaa                                                                        >  1:224686/1‑140 (MQ=38)                                                                 ataaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatg                                                                     >  1:132002/1‑140 (MQ=255)                                                                           tGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAataata                                                           >  1:343397/1‑140 (MQ=255)                                                                                          aGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATATAATTAg                                            >  2:154748/1‑140 (MQ=255)                                                                                                           tagatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACa                           >  2:16387/1‑140 (MQ=255)                                                                                                           tagatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATGATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAGGATACa                           >  2:176772/1‑140 (MQ=255)                                                                                                             gatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACaaa                         <  1:182297/140‑1 (MQ=255)                                                                                                               tagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGCGACTCATTTATTGGGGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATc                       >  1:259290/1‑140 (MQ=25)                                                                                                                      ccacGGACTGCAAAGGGGCATTTTTTTGTTAAGACATTCTATAATGAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTCGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCAaaaaca                >  2:255683/2‑138 (MQ=21)                                                                                                                       cacGGACTGCAAAGGGGCATTATTTTGATAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATAAGATAAAGTAGAAGATACAAAACTAAAAAag               >  2:213151/1‑140 (MQ=255)                                                                                                                         cGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAACagag             >  1:256286/1‑140 (MQ=255)                                                                                                                         cGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGAAGAAGATACAAATCTAAAAAagag             >  1:105885/1‑140 (MQ=255)                                                                                                                          ggACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGACATGCAATTTAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAagagc            >  1:168927/1‑139 (MQ=21)                                                                                                                          ggACTGCAAAGGGGCATTATTTTCTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAaaga            >  2:255041/1‑140 (MQ=255)                                                                                                                                    ggggCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACt  >  2:21425/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         GATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACT  >  CP000730/469652‑469922

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