Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,125,492	C→T	Y15Y (TAC→TAT)	SAUSA300_RS10825 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,125,492	0	C	T	100.0%	80.2 / NA	29	Y15Y (TAC→TAT)	SAUSA300_RS10825	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/4);  total (25/4)

CTATTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAGAGAAACTGTAATATACGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCA  >  NC_007793/2125416‑2125585                                                                             |                                                                                              cTATTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAgggtagatcg                                                                        >  1:6640644/1‑91 (MQ=255)   aTTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAAc                                                                      >  1:6510785/1‑100 (MQ=255)   aTTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAAc                                                                      >  1:5299556/1‑100 (MQ=255)      tCAGCAGAGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAt                                                                   >  1:4117481/1‑100 (MQ=255)      gcagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAt                                                                   >  1:803873/2‑100 (MQ=255)       tagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:4569477/2‑100 (MQ=255)       cagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:1404005/1‑100 (MQ=255)       cagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:4781277/1‑100 (MQ=255)       cagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:3445294/1‑100 (MQ=255)       cagcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:4081601/1‑100 (MQ=255)       atgcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCAtt                                                                  >  1:1286978/3‑100 (MQ=255)        tgcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTc                                                                 >  1:100307/2‑100 (MQ=255)        cgcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTc                                                                 >  1:1425217/2‑100 (MQ=255)         gcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCa                                                                >  1:4564095/1‑100 (MQ=255)         gcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCa                                                                >  1:6126906/1‑100 (MQ=255)         gcagTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCa                                                                >  1:3159391/1‑100 (MQ=255)              gcttGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTAt                                                           >  1:6628691/3‑100 (MQ=255)                 tGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAGGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCGCTTATGAg                                                        >  1:6604513/1‑100 (MQ=255)                       gAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCt                                                  <  1:304800/100‑1 (MQ=255)                              gtaATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGAtt                                           >  1:86080/2‑100 (MQ=255)                               taATGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGAttt                                          <  1:1500195/100‑1 (MQ=255)                                 ctGAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTAt                                        >  1:3615638/2‑100 (MQ=255)                                   gAAAATAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCg                                      >  1:6543564/1‑100 (MQ=255)                                      aaTAACTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTca                                   >  1:1822105/1‑100 (MQ=255)                                         aaCTAATTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTCTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACAt                                >  1:5268472/1‑100 (MQ=255)                                             aaTTGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATg                            >  1:3305120/1‑100 (MQ=255)                                               ttGCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATgcg                          <  1:2665578/100‑1 (MQ=255)                                                 gCAAAATAAAAAAAGAAACTATAGTATATGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATgcgca                        <  1:4951513/100‑2 (MQ=255)                                                                       agtataTGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTca  >  1:4600034/3‑100 (MQ=255)                                                                             |                                                                                              CTATTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAGAGAAACTGTAATATACGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCA  >  NC_007793/2125416‑2125585

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