Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	396,540	G→A	100%	A208T (GCA→ACA)	USA300HOU_0364 →	possible lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	396,540	0	G	A	100.0%	69.5 / NA	24	A208T (GCA→ACA)	USA300HOU_0364	possible lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/8);  total (16/8)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  CP000730/396421‑396677                                                                                                                        |                                                                                                                                          aaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAt                                                                                                                      >  1:156835/1‑141 (MQ=255)                  tCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAATTATTTAAACCTATTTTagag                                                                                                     >  2:228672/1‑141 (MQ=255)                          ccTATTAAATGAAGTTCCACCCTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGCCTTATTTAAACCTATTTTAGAGAAAAAAGa                                                                                             <  1:22941/141‑1 (MQ=255)                                aaaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaa                                                                                       <  1:213306/141‑1 (MQ=255)                                    gAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAtaa                                                                                   >  1:242676/1‑139 (MQ=255)                                      aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAGCCTATTTTAGTGAAAAAGGATAAAAAATTAAg                                                                                 >  1:486/1‑141 (MQ=255)                                             aCTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAATAAAGATAAAAAATTAAGTgatgat                                                                          >  2:23298/1‑141 (MQ=255)                                                 cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATAAAGTGATGATATcc                                                                      >  1:95924/1‑141 (MQ=255)                                                    aaaTCAAAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa                                                                   <  2:65973/141‑1 (MQ=255)                                                      aTCACAGGTGAAGAAGAAATTTATTCACATAAAGATTTATATGATATTAAAGCGCACGTTGAAGGCACACAAAAAATGTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                                                                 >  1:332992/1‑141 (MQ=255)                                                          catgTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGAGAACGTTAAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt                                                             <  2:196869/138‑1 (MQ=255)                                                          caGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt                                                             <  2:130223/141‑1 (MQ=255)                                                                aagaagaaATTTACTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATa                                                       <  2:74970/141‑1 (MQ=255)                                                                 agaagaaATTTGTTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaa                                                      >  2:302827/1‑141 (MQ=255)                                                                 agaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaa                                                      >  2:291623/1‑141 (MQ=255)                                                                                         tttATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATCAAGTGAATCAATTATTGGCTAAatt                              >  1:114773/1‑140 (MQ=255)                                                                                         tttATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAata                              >  1:394841/1‑141 (MQ=255)                                                                                               tGATTTTTAAGCGAACGTGGAAGGCACACGGAAAGTTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGa                        >  2:331483/1‑141 (MQ=255)                                                                                                     tAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATaacaa                  >  1:335736/1‑141 (MQ=255)                                                                                                       aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg                >  1:359440/1‑141 (MQ=255)                                                                                                       aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg                >  1:370403/1‑141 (MQ=255)                                                                                                        aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg               <  1:317511/141‑1 (MQ=255)                                                                                                            aaCGTTGAAGGCACACAAAAAATTTATGAGTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAATGAATCAATTATTGGTGAAATATAAAGATAACAAcggcggg           >  2:45423/1‑140 (MQ=255)                                                                                                                     ggCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCa  <  1:215492/141‑1 (MQ=255)                                                                                                                        |                                                                                                                                          AATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  CP000730/396421‑396677

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