Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	76,969	G→A	intergenic (+1027/‑107)	USA300HOU_0072 → / → USA300HOU_0075	universal stress protein/lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	76,969	0	G	A	100.0%	102.3 / NA	30	intergenic (+1027/‑107)	USA300HOU_0072/USA300HOU_0075	universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/15);  total (15/15)

CTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAGATAAT  >  CP000730/76830‑77111                                                                                                                                            |                                                                                                                                               ctctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAAttt                                                                                                                                       <  2:249788/149‑1 (MQ=255)  tctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTta                                                                                                                                      <  2:210803/149‑1 (MQ=255)       aGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAg                                                                                                                                 <  1:220905/149‑1 (MQ=255)          tCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAt                                                                                                                              <  1:188445/149‑1 (MQ=255)          tCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGATAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAt                                                                                                                              <  2:223786/149‑1 (MQ=255)                gATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGc                                                                                                                        >  2:263094/1‑149 (MQ=255)                 aTTTTCGCCATGTCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCa                                                                                                                       >  1:88841/1‑149 (MQ=255)                    ttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAAc                                                                                                                    >  1:248992/1‑149 (MQ=255)                         cATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAAAATTGa                                                                                                               >  2:187793/1‑149 (MQ=255)                             ccACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTg                                                                                                           <  1:10809/149‑1 (MQ=255)                                   aaTTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATA‑TTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaa                                                                                                    >  1:1879/1‑149 (MQ=255)                                    aTTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaa                                                                                                    >  2:190368/1‑149 (MQ=255)                                               gCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTa                                                                                         <  2:243907/149‑1 (MQ=255)                                               gCTAGCAAGTTAAGCGAACACGGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAATGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAtgggga                                                                                         <  2:103925/149‑3 (MQ=255)                                                 tAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAtt                                                                                       >  2:48730/1‑149 (MQ=255)                                                                   cTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCa                                                                     <  2:82821/149‑1 (MQ=255)                                                                                 tAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGc                                                       >  2:302064/1‑149 (MQ=255)                                                                                   gTGGTTAGTTATATATTTTTACTTTGCAACAGAACCTATAAAAAACATAGGCTTAAAAGGTTAAATTATATAGTATTACTGCAAACATAAAAAAAAAAAAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                                     <  2:138309/149‑1 (MQ=255)                                                                                     ggTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGAt                                                   >  1:100616/1‑149 (MQ=255)                                                                                          gCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAAAAAAAAAAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgag                                              <  2:248992/149‑1 (MQ=255)                                                                                            tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                            >  1:20019/1‑149 (MQ=255)                                                                                            tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                            >  1:131524/1‑149 (MQ=255)                                                                                            tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                            >  2:337785/1‑149 (MQ=255)                                                                                            tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                            >  1:273297/1‑149 (MQ=255)                                                                                                          ttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                              <  1:333743/149‑1 (MQ=255)                                                                                                          ttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                              >  1:145334/1‑149 (MQ=255)                                                                                                                         taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGa               >  2:277237/1‑149 (MQ=255)                                                                                                                                cATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCaaaa        <  1:333442/149‑1 (MQ=255)                                                                                                                                  tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAGa      <  2:247092/149‑1 (MQ=255)                                                                                                                                      cTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAGataat  <  2:88841/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                               CTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAAGATAAT  >  CP000730/76830‑77111

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