Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	747,922	C→T	intergenic (+3/‑100)	SAUSA300_RS03610 → / → SAUSA300_RS03615	cobalamin biosynthesis protein CobW/aldo/keto reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	747,922	0	C	T	100.0%	81.7 / NA	29	intergenic (+3/‑100)	SAUSA300_RS03610/SAUSA300_RS03615	cobalamin biosynthesis protein CobW/aldo/keto reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/1);  total (28/1)

TGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATCAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACTTTTATTGTGAATATTTTGTGACATTGAGGGATGGAGTG  >  NC_007793/747826‑748013                                                                                                 |                                                                                            tGGCATGCAATTACCCTTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGa                                                                                          >  1:2269324/1‑100 (MQ=255) tGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGa                                                                                          >  1:5185898/1‑100 (MQ=255) tGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGa                                                                                          >  1:4800453/1‑100 (MQ=255) tGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGa                                                                                          >  1:47779/1‑100 (MQ=255) tGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGa                                                                                          >  1:2450872/1‑100 (MQ=255)  ggCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGaa                                                                                         >  1:173990/1‑100 (MQ=255)  ggCATGCAATTACCATTAACGATTGTTATTATTGGGGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGaa                                                                                         >  1:2797981/1‑100 (MQ=255)  ggCATGCAATTACCATTAACGATTGTTATTATTGGAGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGaa                                                                                         >  1:4459970/1‑100 (MQ=255)    cATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGaaaa                                                                                       >  1:6217472/1‑100 (MQ=255)    cATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCAAATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGaaaa                                                                                       >  1:5308256/1‑100 (MQ=255)       gCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAAtgt                                                                                    >  1:2222147/1‑100 (MQ=255)       gCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAAtgt                                                                                    >  1:2966375/1‑100 (MQ=255)          aTTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGtgtt                                                                                 >  1:5439285/1‑100 (MQ=255)           ttACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTGCGTAAATTAGAAAATGtgttg                                                                                >  1:6159516/1‑100 (MQ=255)           ttACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGtgttg                                                                                >  1:6315602/1‑100 (MQ=255)           gaaCCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGtgttg                                                                                >  1:2723258/3‑100 (MQ=255)             aCCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTgttgtt                                                                              >  1:1714752/1‑100 (MQ=255)               cATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTc                                                                            >  1:6410234/1‑100 (MQ=255)                                attGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGc                                                           >  1:2870663/1‑100 (MQ=255)                                   ggTGAAACTTTAGATACAAATCACATACGTAATCAATTGGCTATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGcaca                                                        >  1:3150158/1‑100 (MQ=255)                                   ggTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGcaca                                                        >  1:1120774/1‑100 (MQ=255)                                     tGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACAtt                                                      >  1:2264647/1‑100 (MQ=255)                                       aaaCTTTAGATACAAATCACATACGTAATCAATTGGAAATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTaa                                                    >  1:5771622/1‑100 (MQ=255)                                                     cctggcaTACGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACttt                                      >  1:1369921/6‑100 (MQ=255)                                                     aggcccaTACGTAATCAACTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACttt                                      >  1:4535045/6‑100 (MQ=255)                                                             aCGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACTTTTATTGTGa                              >  1:3780521/1‑100 (MQ=255)                                                             aCGTAATCAATTGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACTTTTATTGTGa                              >  1:1858498/1‑100 (MQ=255)                                                                       ttGGATATGCTACAATTTACGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATCCTTTTATTGTGAATATTTtgtg                    <  1:4793008/100‑1 (MQ=255)                                                                                         aCGTAAATTAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACTTTTATTGTGAATATTTTGTGACATTGAGGACTGGAgat  >  1:4547419/1‑98 (MQ=255)                                                                                                 |                                                                                            TGGCATGCAATTACCATTAACGATTGTTATTATTGGTGAAACTTTAGATACAAATCACATACGTAATCAATTGGATATGCTACAATTTACGTAAATCAGAAAATGTGTTGTTTCAATTTAACGGAAATAGCACATTAATTTCTTAATACTTTTATTGTGAATATTTTGTGACATTGAGGGATGGAGTG  >  NC_007793/747826‑748013

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