Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,639,030	C→T	*166* (TAG→TAA)	rzpQ ←	Rz‑like protein, Qin prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,639,030	0	C	T	100.0%	87.5 / NA	28	*166* (TAG→TAA)	rzpQ	Rz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/11);  total (17/11)

AACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTGG  >  NC_000913/1638905‑1639157                                                                                                                              |                                                                                                                                aaCGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCt                                                                                                                    >  2:403554/1‑139 (MQ=255)            cTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCtt                                                                                                         >  2:959631/1‑139 (MQ=255)                       cTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATa                                                                                              <  1:571841/139‑1 (MQ=255)                       cTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATa                                                                                              <  1:403554/139‑1 (MQ=255)                         tttATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGa                                                                                            >  2:1118951/1‑139 (MQ=255)                          ttATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGAt                                                                                           <  1:1118951/139‑1 (MQ=255)                                     ttAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAGGATTCATAGATCTGTTGACAAg                                                                                >  1:913270/1‑139 (MQ=255)                                       aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTc                                                                              >  1:323057/1‑139 (MQ=255)                                       aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCCTAAAGATTCATAGATCTGTTGACAAGTc                                                                              >  2:239430/1‑139 (MQ=255)                                            ttCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTcc                                                                         >  1:29154/1‑139 (MQ=255)                                                             aTGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTa                                                                                                        >  2:920057/1‑91 (MQ=255)                                                             aTGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTa                                                                                                        <  1:920057/91‑1 (MQ=255)                                                                   gagTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTc                                                                          <  1:705063/115‑1 (MQ=255)                                                                   gagTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTc                                                                          >  2:705063/1‑115 (MQ=255)                                                                      tCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAAttt                                               >  1:356120/1‑139 (MQ=255)                                                                      tCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAAttt                                               >  1:1176/1‑139 (MQ=255)                                                                        cTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCa                                             >  1:465469/1‑139 (MQ=255)                                                                           tttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCa                                                                                                <  1:1118128/85‑1 (MQ=255)                                                                           tttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCa                                                                                                >  2:1118128/1‑85 (MQ=255)                                                                           tttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACAATTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCa                                          >  2:870359/1‑139 (MQ=255)                                                                                   aaTTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTg                                  <  1:959631/139‑1 (MQ=255)                                                                                                 cgATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCaa                    >  1:401009/1‑139 (MQ=255)                                                                                                          tGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTTATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGa           <  1:870359/139‑1 (MQ=255)                                                                                                              ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTa                                                                                                        >  2:378708/1‑42 (MQ=255)                                                                                                              ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTa                                                                                                        <  1:378708/42‑1 (MQ=255)                                                                                                              ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCt                                                                        >  2:1007550/1‑74 (MQ=255)                                                                                                              ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCt                                                                        <  1:1007550/74‑1 (MQ=255)                                                                                                                   ttAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTgg  <  2:29154/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                AACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTGG  >  NC_000913/1638905‑1639157

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