Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,212,703	C→T	*135* (TAG→TAA)	ycgX ←	DUF1398 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,212,703	0	C	T	100.0%	80.2 / NA	26	*135* (TAG→TAA)	ycgX	DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

TGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGC  >  NC_000913/1212574‑1212822                                                                                                                                  |                                                                                                                        tGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAAt                                                                                                                >  1:233150/1‑139 (MQ=255)  ggTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATg                                                                                                               <  2:921080/139‑1 (MQ=255)   gTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATgg                                                                                                              >  2:225718/1‑139 (MQ=255)                          aaaaGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCaa                                                                                       <  1:1033588/139‑1 (MQ=255)                                       aaTCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGt                                                                          <  1:516933/139‑1 (MQ=255)                                             ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                    >  2:973434/1‑139 (MQ=255)                                             ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                    >  1:589315/1‑139 (MQ=255)                                                   agaTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaa                                                                                <  1:395596/121‑1 (MQ=255)                                                   agaTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaa                                                                                >  2:395596/1‑121 (MQ=255)                                                       aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCaa                                                                                       >  1:492211/1‑110 (MQ=255)                                                       aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCaa                                                                                       <  2:492211/110‑1 (MQ=255)                                                       aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAAtatata                                                                                 >  2:345850/1‑116 (MQ=255)                                                       aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAAtatata                                                                                 <  1:345850/116‑1 (MQ=255)                                                       aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtag                                                          >  1:156228/1‑139 (MQ=255)                                                           tctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAAc                                                      <  1:225718/139‑1 (MQ=255)                                                              gCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTAAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGc                                                   >  1:169650/1‑139 (MQ=255)                                                                    ttATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGt                                             <  2:1248949/139‑1 (MQ=255)                                                                    ttATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGt                                             >  1:867041/1‑139 (MQ=255)                                                                           cATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGAtt                                      <  2:589315/139‑1 (MQ=255)                                                                                     aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGt                                                 <  2:1055170/118‑1 (MQ=255)                                                                                     aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGt                                                 >  1:1055170/1‑118 (MQ=255)                                                                                                     tGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGc            >  2:327062/1‑139 (MQ=255)                                                                                                       aaGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGccc          >  1:1030255/1‑139 (MQ=255)                                                                                                               cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc  <  2:233150/139‑1 (MQ=255)                                                                                                               cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc  <  2:1030255/139‑1 (MQ=255)                                                                                                               cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc  <  1:973434/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                        TGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGC  >  NC_000913/1212574‑1212822

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