Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,266,814	G→T	R251R (CGA→AGA)	USA300HOU_2140 ←	SAP domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,266,814	0	G	T	100.0%	59.2 / NA	17	R251R (CGA→AGA)	USA300HOU_2140	SAP domain protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/8);  total (9/8)

CTACTTTCTTTTTTTACATTAGACACTACATAGTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCGCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAATTTTTTTAAGAAATCTCTCAAACCAAAGTA  >  CP000730/2266672‑2266952                                                                                                                                               |                                                                                                                                           cTACTTTCTTTTTTTACATTAGACACTACATAGTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACa                                                                                                                                    <  1:289600/151‑1 (MQ=255)                      gACACTACATAGTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGa                                                                                                               >  2:62115/1‑151 (MQ=255)                       acacTACATAGTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGAt                                                                                                              <  2:237374/151‑1 (MQ=255)                                 gTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTc                                                                                                    <  1:42389/151‑1 (MQ=255)                                   ttATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAg                                                                                                  >  2:48781/1‑151 (MQ=255)                                       gAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCg                                                                                              >  1:10036/1‑151 (MQ=255)                                       gAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCg                                                                                              >  1:27772/1‑151 (MQ=255)                                       gAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCg                                                                                              >  1:303172/1‑151 (MQ=255)                                                              ccAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCtat                                                                       <  1:142905/151‑1 (MQ=255)                                                              acAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGCGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCtat                                                                       <  2:207354/150‑1 (MQ=255)                                                                        ataATCTTGTGTGTCCTTTTTAAATAATGCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTCCAAATCCATATTCCCTTTa                                                             >  2:263016/1‑151 (MQ=255)                                                                                            tAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGGTGATAAATCGATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAAtt                                         >  1:112052/1‑151 (MQ=255)                                                                                              aataatTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTAc                                       >  1:124791/1‑151 (MQ=255)                                                                                                         cATAATATTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAAtttt                            >  1:18002/1‑151 (MQ=255)                                                                                                          ataataTTATTGGATAAATTTTCATTATATATTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAAttttt                           <  2:95370/151‑1 (MQ=255)                                                                                                                                   tatataTTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAATTTTTTTAAGAAATCTCTCAAACCAAAGTa  <  1:117703/151‑1 (MQ=255)                                                                                                                                   tatataTTAATCTCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAATTTTTTTAAGAAATCTCTCAAACCAAAGTa  <  1:1910/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                           CTACTTTCTTTTTTTACATTAGACACTACATAGTTTATGAACTTTTCAATCAATTGTTCTTCCAAATCACTATAATCTTGTGTGTCCTTTTTAAATAATTCAAACATAATATTATTGGATAAATTTTCATTATATATTAATCGCTCGTACATACCAGAAGCGCTAGAGTTGATAAATCTATTCAGATCGTTTTGTATTTCGTTCAAATCTATATTCCCTTTACTATCGTAATATGAATAATTACTCATTAATTTTTTTAAGAAATCTCTCAAACCAAAGTA  >  CP000730/2266672‑2266952

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