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%	106.3 / NA	34	*166* (TAG→TAA)	rzpQ	Rz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

TAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTGGT  >  NC_000913/1638899‑1639158                                                                                                                                    |                                                                                                                                 taataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACAttt                                                                                                                           >  2:1347710/1‑139 (MQ=255)      aaaCGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGc                                                                                                                      >  1:1275114/1‑139 (MQ=255)           tgtCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttttt                                                                                                                 <  1:1347710/139‑1 (MQ=255)                   tGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGAt                                                                                                                        >  2:286197/1‑124 (MQ=255)                   tGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGAt                                                                                                                        <  1:286197/124‑1 (MQ=255)                    gATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTaa                                                                                                        <  2:1275114/139‑1 (MQ=255)                      ttAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAg                                                                                                      <  2:83971/139‑1 (MQ=255)                          gACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGAt                                                                                                                        <  1:1128368/117‑1 (MQ=255)                          gACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGAt                                                                                                                        >  2:1128368/1‑117 (MQ=255)                               tttATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGa                                                                                             >  1:1172523/1‑139 (MQ=255)                                  aTCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCt                                                                                          <  2:1172523/139‑1 (MQ=255)                                       aCCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttttt                                                                                                                 <  2:1072563/111‑1 (MQ=255)                                       aCCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttttt                                                                                                                 >  1:1072563/1‑111 (MQ=255)                                            tAAGATTTCAGGAATTTTGGCTCATGGAAGGGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGt                                                                                <  1:919991/139‑1 (MQ=255)                                            tAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGt                                                                                <  1:761639/139‑1 (MQ=255)                                             aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTCTTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTc                                                                               >  2:1209907/1‑139 (MQ=255)                                             aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGa                                                                                                                         <  2:372260/97‑1 (MQ=255)                                             aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGa                                                                                                                         >  1:372260/1‑97 (MQ=255)                                             aaGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTc                                                                               >  2:1355561/1‑139 (MQ=255)                                                  ttCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGc                                                                                                                      <  1:882037/95‑1 (MQ=255)                                                  ttCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGc                                                                                                                      >  2:882037/1‑95 (MQ=255)                                                   tCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCt                                                                         >  2:260611/1‑139 (MQ=255)                                                    cAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTg                                                                        <  1:260611/139‑1 (MQ=255)                                                    cAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTg                                                                        <  1:1355561/139‑1 (MQ=255)                                                                               tttttATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACATTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAg                                             >  1:328562/1‑139 (MQ=255)                                                                                 tttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATa                                                                                               >  1:743530/1‑87 (MQ=255)                                                                                 tttatttaAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATa                                                                                               <  2:743530/87‑1 (MQ=255)                                                                                            tttACATTCCGCGATGTAAATGTTCAGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGc                                <  1:1209907/139‑1 (MQ=255)                                                                                                        gATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAat                    >  2:908699/1‑139 (MQ=255)                                                                                                          tGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAatat                  >  2:1473930/1‑139 (MQ=255)                                                                                                                    ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTctgc                          >  1:259246/1‑121 (MQ=255)                                                                                                                    ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTctgc                          <  2:259246/121‑1 (MQ=255)                                                                                                                      gATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGt      >  2:191381/1‑139 (MQ=255)                                                                                                                          tAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTGGt  <  1:908699/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 TAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGTTGGT  >  NC_000913/1638899‑1639158

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