Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,390,754	G→T	A252A (GCC→GCA)	rsgA ←	ribosome small subunit‑dependent GTPase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,390,754	0	G	T	100.0%	67.7 / NA	20	A252A (GCC→GCA)	rsgA	ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/5);  total (15/5)

AACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTT  >  NC_000913/4390692‑4390818                                                               |                                                                 aaCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGt                                                             <  1:2880692/68‑1 (MQ=255)    tcccGCACTCGTGGGGAGTCAATCACATCCCCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCgt                                                          <  1:720463/65‑1 (MQ=255)        gCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCgtggtg                                                      >  1:2318066/1‑68 (MQ=255)          aCTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGtgtg                                                    >  1:1100020/1‑68 (MQ=255)            tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCt                                                  >  1:3842524/1‑68 (MQ=255)            tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCt                                                  >  1:2418889/1‑68 (MQ=255)             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATGCAGCCAAGCTGCGGTCGCGGGGTGCTg                                                 >  1:3393223/1‑68 (MQ=255)                  ggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGAccct                                            >  1:2956151/1‑66 (MQ=255)                  ggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGGGGTGTGCTGACCGa                                            >  1:1972660/1‑68 (MQ=255)                   gAGTCAATCACATCCCCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAg                                           <  1:1787003/68‑1 (MQ=255)                   gAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAg                                           <  1:1739591/68‑1 (MQ=255)                    aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGc                                          >  1:1625393/1‑68 (MQ=255)                    aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGc                                          >  1:2971625/1‑68 (MQ=255)                    aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGCGTGCTGACCGAGc                                          >  1:3291615/1‑68 (MQ=255)                    aGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGCGGTGTGCTGACCGAGc                                          >  1:911883/1‑68 (MQ=255)                         atcacatcacCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGa                                     >  1:2740626/1‑68 (MQ=255)                            acatcacCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGtt                                  >  1:542727/1‑68 (MQ=255)                              atcccCGCCGTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGt                                <  1:2947183/68‑1 (MQ=255)                                   cgcctTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGCCCGAGCCCCGAGTTGTCagag                           <  1:3104823/68‑1 (MQ=255)                                    gccgTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAgaga                          >  1:2533403/1‑68 (MQ=255)                                    gccgTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAgaga                          >  1:1595551/1‑68 (MQ=255)                                    gccgTGCGGGAAGTGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAaaga                          >  1:2955311/1‑68 (MQ=255)                                        ggCGGGAAGGGATACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                      <  1:1384145/67‑1 (MQ=255)                                                  gATGCAGCCGAGCTGCGGGCGTGGTGTGCGGGCCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAgg            <  1:2604158/68‑1 (MQ=255)                                                   aTACAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGa           <  1:1609641/68‑1 (MQ=255)                                                     aCAGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtc         >  1:4022221/1‑68 (MQ=255)                                                       aGCCGAGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGGTGTCAGAGATATCGTTGGCCAGGCtctc       >  1:2165590/1‑68 (MQ=255)                                                        gCCGAGCTGCGGTCGTGGTGTGCTGGCCGACCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATctct      >  1:466284/1‑68 (MQ=255)                                                            aGCTGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCttttt  <  1:3492284/68‑1 (MQ=255)                                                               |                                                                 AACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTT  >  NC_000913/4390692‑4390818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑