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%	123.0 / NA	35	Y273H (TAC→CAC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  NC_000913/3555904‑3556164                                                                                                                                  |                                                                                                                                    ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc                                                                                                                            <  2:418308/139‑1 (MQ=255)  caccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAcc                                                                                                                           <  1:702138/139‑1 (MQ=255)          tatTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCAc                                                                                                                   <  2:540697/139‑1 (MQ=255)                    cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc                                                                                                         >  1:559486/1‑139 (MQ=255)                              gAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTACCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt                                                                                               >  1:680898/1‑139 (MQ=255)                                   aTGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTCCACCAACTGTGCCGCAACCACCTCGGCa                                                                                          >  2:639016/1‑139 (MQ=255)                                    tGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCa                                                                                                                             >  2:945567/1‑103 (MQ=255)                                    tGGGCGACCATAAATTCCCCCGCGCCCTCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTCTCa                                                                                                                             <  1:945567/103‑1 (MQ=255)                                     gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACt                                                                                        >  1:825767/1‑139 (MQ=255)                                     gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACt                                                                                        >  2:1082154/1‑139 (MQ=255)                                      ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTg                                                                                       >  2:182834/1‑139 (MQ=255)                                         gACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:680898/139‑1 (MQ=255)                                           ccGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcgcg                                                                                  <  1:116420/139‑1 (MQ=255)                                            cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          >  1:598587/1‑98 (MQ=255)                                            cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          <  2:598587/98‑1 (MQ=255)                                                    cccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCAcc                                                                         <  2:825767/139‑1 (MQ=255)                                                         gcgcCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcgacga                                                                    >  2:597305/1‑139 (MQ=255)                                                               gccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          <  1:1347209/79‑1 (MQ=255)                                                               gccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                          >  2:1347209/1‑79 (MQ=255)                                                                         tCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGt                                                    <  1:597305/139‑1 (MQ=255)                                                                                cACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    <  2:1291285/84‑1 (MQ=255)                                                                                cACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                                    >  1:1291285/1‑84 (MQ=255)                                                                                     gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    >  1:763009/1‑95 (MQ=255)                                                                                     gTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg                                                                                    <  2:763009/95‑1 (MQ=255)                                                                                             gCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACtt                                >  2:786004/1‑139 (MQ=255)                                                                                                 gATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAg                            <  1:544970/139‑1 (MQ=255)                                                                                                 gATATTCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAg                            <  1:182834/139‑1 (MQ=255)                                                                                                      tCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga                       <  2:77154/139‑1 (MQ=255)                                                                                                         cccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                                       >  1:841486/1‑88 (MQ=255)                                                                                                         cccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga                                                                       <  2:841486/88‑1 (MQ=255)                                                                                                            aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAg                            >  2:405225/1‑128 (MQ=255)                                                                                                            aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAg                            <  1:405225/128‑1 (MQ=255)                                                                                                               gccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc                                                                                     <  1:631100/68‑1 (MQ=255)                                                                                                               gccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc                                                                                     >  2:631100/1‑68 (MQ=255)                                                                                                                           gCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGt  <  2:1351070/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    CCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGT  >  NC_000913/3555904‑3556164

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