Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,011	C→T	*195* (TAG→TAA)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,011	0	C	T	100.0%	193.3 / NA	60	*195* (TAG→TAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/31);  total (29/31)

GCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAA  >  NC_000913/1145884‑1146145                                                                                                                                |                                                                                                                                       gCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAg                                                                                                                             >  1:1215642/1‑139 (MQ=255) gCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAAGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAg                                                                                                                             >  1:1209041/1‑139 (MQ=255)  cGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGc                                                                                                                            <  1:1057386/139‑1 (MQ=255)                   cAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCg                                                                                                           <  1:1023721/139‑1 (MQ=255)                        atatAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      <  2:716166/139‑1 (MQ=255)                        atatAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      <  1:31218/139‑1 (MQ=255)                          ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                                    >  2:612922/1‑139 (MQ=255)                          ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                                    >  1:136913/1‑139 (MQ=255)                                 tGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGac                                                                                             >  1:536626/1‑139 (MQ=255)                                 ggATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGac                                                                                             <  1:312302/138‑1 (MQ=255)                                  gATAATGGGCTGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAGAATTATCCCATCAGCGGGTTTTTCCCGTCCCGACGTAACATCTGa                                                                                              >  2:727380/1‑137 (MQ=255)                                  gATAATGGGCTGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGa                                                                                              <  1:727380/137‑1 (MQ=255)                                  gATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACAt                                                                                                  >  2:895633/1‑133 (MQ=255)                                  gATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACAt                                                                                                  <  1:895633/133‑1 (MQ=255)                                   aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                             <  2:345973/121‑1 (MQ=255)                                   aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                             >  1:345973/1‑121 (MQ=255)                                    tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGAcaca                                                                                          >  2:1036964/1‑139 (MQ=255)                                        gggCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAgcgc                                                                                      >  1:1028361/1‑139 (MQ=255)                                          gCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttccc                                                                                                            <  1:1200643/115‑1 (MQ=255)                                          gCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttccc                                                                                                            >  2:1200643/1‑115 (MQ=255)                                          gCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGa                                                                                    >  1:324929/1‑139 (MQ=255)                                           cGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGAt                                                                                   <  2:726245/139‑1 (MQ=255)                                             gTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTa                                                                                 >  1:497348/1‑139 (MQ=255)                                                ggCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCg                                                                              >  1:12324/1‑139 (MQ=255)                                                        aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaacca                                                                      >  2:366842/1‑139 (MQ=255)                                                           gCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaaccaacc                                                                   <  1:258176/139‑1 (MQ=255)                                                            cGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                                  <  2:324929/139‑1 (MQ=255)                                                             gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                                 <  2:497348/139‑1 (MQ=255)                                                             gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                                 <  2:536626/139‑1 (MQ=255)                                                             gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                                 <  2:12324/139‑1 (MQ=255)                                                             gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                                 <  1:1036964/139‑1 (MQ=255)                                                             gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                                 <  1:366842/139‑1 (MQ=255)                                                              tttACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                                  <  2:317729/137‑1 (MQ=255)                                                              tttACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                                  >  1:317729/1‑137 (MQ=255)                                                                tACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCgt                                                              >  2:551420/1‑139 (MQ=255)                                                                     cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttccc                                                                                                            >  2:1001808/1‑88 (MQ=255)                                                                     cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttccc                                                                                                            <  1:1001808/88‑1 (MQ=255)                                                                     cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAg                                                         >  2:293036/1‑139 (MQ=255)                                                                      gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGc                                                                                                                            <  2:57865/71‑1 (MQ=255)                                                                      gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGc                                                                                                                            >  1:57865/1‑71 (MQ=255)                                                                      gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      <  1:91245/93‑1 (MQ=255)                                                                      gcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      >  2:91245/1‑93 (MQ=255)                                                                            ggCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACg                                                  <  1:551420/139‑1 (MQ=255)                                                                             gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      >  1:693680/1‑86 (MQ=255)                                                                             gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                      <  2:693680/86‑1 (MQ=255)                                                                             gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCACAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                                <  2:953031/124‑1 (MQ=255)                                                                             gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCACAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                                >  1:953031/1‑124 (MQ=255)                                                                                    gCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTa                                          <  2:1215642/139‑1 (MQ=255)                                                                                    gCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTa                                          <  2:1209041/139‑1 (MQ=255)                                                                                          aTTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCt                                    >  2:145927/1‑139 (MQ=255)                                                                                                  aTGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGAt                                                                                   >  2:126866/1‑84 (MQ=255)                                                                                                  aTGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGAt                                                                                   <  1:126866/84‑1 (MQ=255)                                                                                                                   ccGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc           <  1:145927/139‑1 (MQ=255)                                                                                                                     ggCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCg                                                                              >  2:1062036/1‑70 (MQ=255)                                                                                                                     ggCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCg                                                                              <  1:1062036/70‑1 (MQ=255)                                                                                                                      gcgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAct        >  1:430827/1‑139 (MQ=255)                                                                                                                       cgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActc       >  1:333740/1‑139 (MQ=255)                                                                                                                         cAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCaa                  <  1:977616/126‑1 (MQ=255)                                                                                                                         cAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCaa                  >  2:977616/1‑126 (MQ=255)                                                                                                                            aaaaTTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACACTCAAACAGCGTAATGCCAAATCCTTCACTCTTaaa  <  2:1028361/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       GCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAA  >  NC_000913/1145884‑1146145

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