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

CACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  NC_000913/3555902‑3556153                                                                                                                                     |                                                                                                                          caccacCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTc                                                                                                                    <  2:790982/139‑1 (MQ=255)     accGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                <  1:89441/139‑1 (MQ=255)       cgcgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                              <  1:1865897/139‑1 (MQ=255)              ttGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCa                                                                                                       >  2:1567029/1‑139 (MQ=255)                ggTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAAc                                                                                                     <  1:723405/139‑1 (MQ=255)                      cagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc                                                                                               >  2:1858700/1‑139 (MQ=255)                           ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                          <  2:178586/139‑1 (MQ=255)                           ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                          >  1:1436456/1‑139 (MQ=255)                           ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac                                                                                          >  2:1332274/1‑139 (MQ=255)                                  gggATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCg                                                                                   >  1:342450/1‑139 (MQ=255)                                       gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCTCCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACt                                                                              >  1:1330881/1‑139 (MQ=255)                                       gagCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCGGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACt                                                                              >  1:447265/3‑139 (MQ=255)                                        ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                >  2:1109941/1‑104 (MQ=255)                                        ggCGACCGTAAATTCCCCCGCGACCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                <  1:1109941/104‑1 (MQ=255)                                          cGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc                                                                           <  1:1598924/139‑1 (MQ=255)                                            aCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                >  1:984686/1‑100 (MQ=255)                                            aCCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct                                                                                                                <  2:984686/100‑1 (MQ=255)                                                        cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCCACCGCCGCCGTGCATGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcg                                                             >  2:1169032/1‑139 (MQ=255)                                                         ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcgac                                                           <  2:1330881/139‑1 (MQ=255)                                                          cgcgCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                            >  2:634239/1‑106 (MQ=255)                                                          cgcgCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCa                                                                                            <  1:634239/106‑1 (MQ=255)                                                          cgcgCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcgacg                                                          >  2:1347862/1‑139 (MQ=255)                                                           gcgcCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTACTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcgacga                                                         <  1:1347862/139‑1 (MQ=255)                                                            cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCAcc                                                              >  2:1240744/1‑133 (MQ=255)                                                            cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCAcc                                                              <  1:1240744/133‑1 (MQ=255)                                                                    aGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGc                                                >  2:522722/1‑139 (MQ=255)                                                                       gcCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGtt                                             <  1:239921/139‑1 (MQ=255)                                                                       gcCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGtt                                             <  1:693876/139‑1 (MQ=255)                                                                       gcCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGtt                                             <  1:1858700/139‑1 (MQ=255)                                                                         cATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCg                                           <  2:1610119/139‑1 (MQ=255)                                                                                       gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcgacga                                                         <  2:25968/111‑1 (MQ=255)                                                                                       gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcgacga                                                         >  1:25968/1‑111 (MQ=255)                                                                                       gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcga                                                            >  2:1712504/1‑108 (MQ=255)                                                                                       gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACcga                                                            <  1:1712504/108‑1 (MQ=255)                                                                                       gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACAACC‑CCGCACCGGCGCCCCTTTCCCCCGCCGCCAAAAAA‑CGTTCGGTGAATATTTCACtt                            >  1:1228746/1‑139 (MQ=255)                                                                                             cAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAc                       >  1:26409/1‑139 (MQ=255)                                                                                             cAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTCTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAc                       <  1:1332274/139‑1 (MQ=255)                                                                                               gCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACtt                     >  1:1094185/1‑139 (MQ=255)                                                                                                  aGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                              <  1:808486/48‑1 (MQ=255)                                                                                                  aGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACctct                                                                                                              >  2:808486/1‑48 (MQ=255)                                                                                                  aGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCaa                  >  2:759870/1‑139 (MQ=255)                                                                                                  aGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCaa                  >  2:59241/1‑139 (MQ=255)                                                                                                        t‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            <  1:1567029/139‑1 (MQ=255)                                                                                                        t‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga            <  1:522722/139‑1 (MQ=255)                                                                                                           ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCAc                                                                                                         <  2:1067430/45‑1 (MQ=255)                                                                                                           ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCAc                                                                                                         >  1:1067430/1‑45 (MQ=255)                                                                                                              cACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg‑cg                                                                       >  2:838638/1‑75 (MQ=255)                                                                                                              cACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAcg‑cg                                                                       <  1:838638/75‑1 (MQ=255)                                                                                                              cACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTa       >  1:683220/1‑139 (MQ=255)                                                                                                               aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCttt                                                                   <  2:1098524/78‑1 (MQ=255)                                                                                                               aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACGGACG‑CGCttt                                                                   >  1:1098524/1‑78 (MQ=255)                                                                                                                  gccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAc   >  1:218176/1‑139 (MQ=255)                                                                                                                   ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:218176/139‑1 (MQ=255)                                                                                                                   ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  >  1:292839/1‑139 (MQ=255)                                                                                                                   ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  2:1436456/139‑1 (MQ=255)                                                                                                                   ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc  <  1:1169032/139‑1 (MQ=255)                                                                                                                    cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc                                                                           >  1:602620/1‑66 (MQ=255)                                                                                                                    cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGAc                                                                           <  2:602620/66‑1 (MQ=255)                                                                                                                    cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACtt                            <  2:992463/112‑1 (MQ=255)                                                                                                                    cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACtt                            >  1:992463/1‑112 (MQ=255)                                                                                                                              gCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTg                                                                             <  2:875349/54‑1 (MQ=255)                                                                                                                              gCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTg                                                                             >  1:875349/1‑54 (MQ=255)                                                                                                                                     |                                                                                                                          CACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACG‑CGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC  >  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: ‑