Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,117,002	T→C	T42T (ACA→ACG)	wcaL ←	putative glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,117,002	0	T	C	100.0%	101.5 / NA	28	T42T (ACA→ACG)	wcaL	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/14);  total (14/14)

GCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTC  >  NC_000913/2116860‑2117147                                                                                                                                               |                                                                                                                                                  gCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCg                                                                                                                                             <  2:596979/149‑1 (MQ=255)             tACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAgcgc                                                                                                                               >  1:288185/1‑151 (MQ=255)             tACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAgcg                                                                                                                                >  1:175144/1‑150 (MQ=255)               cGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGcgcg                                                                                                                              >  1:460930/1‑150 (MQ=255)                gATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCgag                                                                                                                            <  1:88185/151‑1 (MQ=255)                           cgcAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCt                                                                                                                 <  2:446571/151‑1 (MQ=255)                                     ggctggctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATccc                                                                                                         >  2:595150/1‑149 (MQ=255)                                               gTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAAc                                                                                             <  1:532004/151‑1 (MQ=255)                                                   cgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCgg                                                                                         >  1:578007/1‑151 (MQ=255)                                                      aGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTaa                                                                                      <  1:150860/151‑1 (MQ=255)                                                       gTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAATGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAAt                                                                                     <  1:520119/151‑1 (MQ=255)                                                                     cccGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGa                                                                       <  1:241377/151‑1 (MQ=255)                                                                       cGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAg                                                                     <  2:121155/151‑1 (MQ=255)                                                                          aGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGtt                                                                  >  1:338627/1‑151 (MQ=255)                                                                          aGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGt                                                                   <  2:201082/150‑1 (MQ=255)                                                                               cGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTga                                                             <  2:8952/151‑1 (MQ=255)                                                                                                       cagcaAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTaaa                                       >  1:540204/1‑149 (MQ=255)                                                                                                         gcaAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTaaaa                                      >  2:329945/1‑148 (MQ=255)                                                                                                                tGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGtt                                                                  <  1:599544/113‑1 (MQ=255)                                                                                                                tGTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGtt                                                                  >  2:599550/1‑113 (MQ=255)                                                                                                                 gTATTTCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGAcc                           >  2:249728/1‑151 (MQ=255)                                                                                                                     ttCGTCCATGCCGCGTGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTc                        <  1:654021/150‑1 (MQ=255)                                                                                                                          ccATGCCGCGCGGGTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTc                   >  2:629663/1‑150 (MQ=255)                                                                                                                                  cgTGGGAGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCg          >  1:135396/1‑151 (MQ=255)                                                                                                                                     gggTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGGCCTTCATTTCGCCGCCCCGATg       <  1:129642/151‑1 (MQ=255)                                                                                                                                     gggTGTTTTGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCCCCTCCCCGAt        <  1:10014/150‑1 (MQ=255)                                                                                                                                          ttttGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTc  >  1:669447/1‑151 (MQ=255)                                                                                                                                           tttGCGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTc  >  2:290584/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                  GCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGCCTTTTTGCAGCGCGAGAATCTCTACCTCAAATCCCATATCAATAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTC  >  NC_000913/2116860‑2117147

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