Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,556,033	A→G	Y273H (TAC→CAC)	rtcA ←	RNA 3'‑terminal phosphate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,033	0	A	G	100.0%	153.7 / NA	43	Y273H (TAC→CAC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/20);  total (23/20)

CACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  NC_000913/3555902‑3556153                                                                                                                                    |                                                                                                                         caccacCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTc                                                                                                                   <  2:946824/139‑1 (MQ=255)                 gTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACt                                                                                                   >  2:567514/1‑139 (MQ=255)                     gcagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  2:338616/1‑123 (MQ=255)                     gcagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  1:338616/123‑1 (MQ=255)                      cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  2:459828/1‑122 (MQ=255)                      cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  2:584904/122‑1 (MQ=255)                      cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  1:584904/1‑122 (MQ=255)                      cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  1:459828/122‑1 (MQ=255)                        gATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                             <  1:678226/122‑1 (MQ=255)                        gATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                             >  2:678226/1‑122 (MQ=255)                                     aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  1:529880/107‑1 (MQ=255)                                     aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  2:529880/1‑107 (MQ=255)                                     aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCa                                                                               >  2:502125/1‑139 (MQ=255)                                       gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACt                                                                             >  2:471604/1‑139 (MQ=255)                                        ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               <  2:595322/104‑1 (MQ=255)                                        ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                               >  1:595322/1‑104 (MQ=255)                                           gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                         >  2:775233/1‑139 (MQ=255)                                             ccGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACGTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                       >  1:605433/1‑139 (MQ=255)                                              cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                       >  1:4885/1‑138 (MQ=255)                                              cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                       <  2:4885/138‑1 (MQ=255)                                                       ccccGCGCCCGCCAGCGCCATCGGTAGCTCCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt                                                                                    >  2:625418/1‑116 (MQ=255)                                                       ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt                                                                                    <  1:625418/116‑1 (MQ=255)                                                       ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcg                                                             >  2:549050/1‑139 (MQ=255)                                                        cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            >  2:851892/1‑139 (MQ=255)                                                            cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            <  2:249414/135‑1 (MQ=255)                                                            cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                            >  1:249414/1‑135 (MQ=255)                                                               ccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACaaa                                                     <  1:549050/139‑1 (MQ=255)                                                                       gcCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt                                             <  1:775233/139‑1 (MQ=255)                                                                                       gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         >  2:804918/1‑79 (MQ=255)                                                                                       gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                         <  1:804918/79‑1 (MQ=255)                                                                                        ttGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACtt                            >  2:791789/1‑139 (MQ=255)                                                                                         tGGTCAGCAAGATATTCCCCCACCGCCGCCGTACTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACttt                           >  1:599246/1‑139 (MQ=255)                                                                                           gTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCg                         >  1:918860/1‑139 (MQ=255)                                                                                                cAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTc                    >  1:345242/1‑139 (MQ=255)                                                                                                       ttCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgag             <  1:171495/139‑1 (MQ=255)                                                                                                        tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCTCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            <  1:6061/139‑1 (MQ=255)                                                                                                        tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            <  2:562685/139‑1 (MQ=255)                                                                                                        tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            <  1:256549/139‑1 (MQ=255)                                                                                                         cccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAc           >  2:577141/1‑139 (MQ=255)                                                                                                               ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTTACTTTCGACTTCAAGCGAGACGGTAtt     >  1:406237/1‑139 (MQ=255)                                                                                                                  ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:406237/139‑1 (MQ=255)                                                                                                                  ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:605433/139‑1 (MQ=255)                                                                                                                  ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  1:791789/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                         CACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  NC_000913/3555902‑3556153

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