Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,629,478	G→A	H21Y (CAT→TAT)	rbbA ←	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,629,478	0	G	A	100.0%	95.3 / NA	31	H21Y (CAT→TAT)	rbbA	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

GACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCT  >  NC_000913/3629384‑3629601                                                                                               |                                                                                                                            gACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCa                                                                                 >  2:503129/1‑139 (MQ=255)             ccGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                     <  1:146638/139‑1 (MQ=255)                tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                  >  1:923466/1‑139 (MQ=255)                       cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                           <  1:503129/139‑1 (MQ=255)                          tCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                              >  2:635607/1‑117 (MQ=255)                          tCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                              <  1:635607/117‑1 (MQ=255)                                     aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                             <  2:923466/139‑1 (MQ=255)                                      tACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                            >  1:598818/1‑139 (MQ=255)                                             ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCAcg                                                                                                                   <  1:398482/61‑1 (MQ=255)                                             ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCAcg                                                                                                                   >  2:398482/1‑61 (MQ=255)                                             ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                     >  2:322734/1‑139 (MQ=255)                                             ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGGTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATCCGTCATTGCGGCAACCTCACCACGAGgtcgtcgg                                     >  1:858894/1‑137 (MQ=255)                                                cGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGcc                                  >  1:402708/1‑139 (MQ=255)                                                        cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                        <  1:357061/125‑1 (MQ=255)                                                        cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                        >  2:357061/1‑125 (MQ=255)                                                         gagaGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGtt                         <  1:322734/139‑1 (MQ=255)                                                              tGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttca                    <  2:598818/139‑1 (MQ=255)                                                                atatTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcatt                  >  1:206111/1‑139 (MQ=255)                                                                      ttCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg            <  1:887908/139‑1 (MQ=255)                                                                          gcgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                             >  1:753676/1‑102 (MQ=255)                                                                          gcgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                             <  2:753676/102‑1 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                    >  1:512556/1‑62 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                    <  2:512556/62‑1 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCGAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGtt                         <  1:651594/121‑1 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCGAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGtt                         >  2:651594/1‑121 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                           >  1:838097/1‑87 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                           <  2:838097/87‑1 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc              >  1:798860/1‑132 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc              <  2:798860/132‑1 (MQ=255)                                                                           cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc       >  1:754029/1‑139 (MQ=255)                                                                                aggTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCt  >  2:580431/2‑139 (MQ=255)                                                                                               |                                                                                                                            GACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCT  >  NC_000913/3629384‑3629601

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