Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,125,476	G→A	R10K (AGA→AAA)	SAUSA300_RS10825 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,125,476	0	G	A	100.0%	77.0 / NA	28	R10K (AGA→AAA)	SAUSA300_RS10825	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/4);  total (24/4)

CTATTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAGAGAAACTGTAATATACGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCG  >  NC_007793/2125416‑2125563                                                             |                                                                                        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)                                                             |                                                                                        CTATTTCAGCAGTGTTTGAAAGGAAGTTTATAATGAAAATAACTAATTGCAAAATAAAAAGAGAAACTGTAATATACGAAGTTTTAACTAGTGGTAATCAACCATTCACTTATGAGTTACCTAAAGATTTATCGTCACATAATGCGCG  >  NC_007793/2125416‑2125563

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