Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,128,386	(C)5→6	coding (503/549 nt)	wcaF ←	putative acyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,128,381	1	.	C	100.0%	150.1 / NA	43	N170D (AAT→GAT)	wcaF	putative acyl transferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (26/17);  total (26/17)

CCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATT‑CCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAA  >  NC_000913/2128260‑2128515                                                                                                                           |                                                                                                                                       ccAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACAt                                                                                                                       >  2:107275/1‑140 (MQ=255)    ggTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTg                                                                                                                    >  2:217840/1‑140 (MQ=255)           aCCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAg                                                                                                             >  1:415882/1‑140 (MQ=255)               tCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGAt                                                                                                         >  1:7253/1‑140 (MQ=255)               tCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGAt                                                                                                         >  1:575020/1‑140 (MQ=255)                                 gTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGt                                                                                       <  2:675152/140‑1 (MQ=255)                                        gagCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCa                                                                                <  1:169215/140‑1 (MQ=255)                                         agCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCac                                                                               >  1:346593/1‑140 (MQ=255)                                           cGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCacga                                                                             >  2:250210/1‑140 (MQ=255)                                           cGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCacga                                                                             >  2:1981/1‑140 (MQ=255)                                            gACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCacgac                                                                            <  1:250210/140‑1 (MQ=255)                                              cTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACgg                                                                          >  1:331023/1‑140 (MQ=255)                                                       tGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGc                                                                 <  2:346593/140‑1 (MQ=255)                                                              ttCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtgt                                                          >  2:285316/1‑140 (MQ=255)                                                               tCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATtgtg                                                         <  2:7253/140‑1 (MQ=255)                                                                ccTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGa                                                        >  2:672419/1‑140 (MQ=255)                                                                  tctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACg                                                      >  1:328446/1‑140 (MQ=255)                                                                   ctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGc                                                     >  2:200718/1‑140 (MQ=255)                                                                   ctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGc                                                     >  2:212944/1‑140 (MQ=255)                                                                   ctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGc                                                     >  2:526406/1‑140 (MQ=255)                                                                   ctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGc                                                     >  1:533665/1‑140 (MQ=255)                                                                   ctGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGc                                                     >  2:418953/1‑140 (MQ=255)                                                                    tGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGcc                                                    <  2:533665/140‑1 (MQ=255)                                                                    tGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGcc                                                    <  2:575020/140‑1 (MQ=255)                                                                    tGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGcc                                                    >  2:382640/1‑140 (MQ=255)                                                                         ttttGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTgggg                                               <  1:627392/140‑1 (MQ=255)                                                                            tGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGcaa                                            <  1:670751/140‑1 (MQ=255)                                                                                           tCAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTg                             <  2:415882/140‑1 (MQ=255)                                                                                            cAACGCGTTCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTgc                            >  1:696455/1‑140 (MQ=255)                                                                                                   ttCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGccagccag                     >  1:681431/1‑140 (MQ=255)                                                                                                   ttCGCGTATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGccagccag                     >  2:22670/1‑140 (MQ=255)                                                                                                         tATCACCACTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTtc               <  1:645075/140‑1 (MQ=255)                                                                                                            caccacTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCg            <  2:737307/140‑1 (MQ=255)                                                                                                            caccacTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCg            <  1:526406/140‑1 (MQ=255)                                                                                                              ccacTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGcc          >  1:163534/1‑140 (MQ=255)                                                                                                                 cTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAAt       <  2:328446/140‑1 (MQ=255)                                                                                                                 cTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAAt       <  1:494601/140‑1 (MQ=255)                                                                                                                 cTGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAAt       <  1:369228/140‑1 (MQ=255)                                                                                                                  tGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATc      >  2:387702/1‑140 (MQ=255)                                                                                                                  tGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATc      >  2:544070/1‑140 (MQ=255)                                                                                                                  tGCGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATc      >  1:481389/1‑140 (MQ=255)                                                                                                                    cGGGATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATcac    >  1:211843/1‑140 (MQ=255)                                                                                                                      ggATTCCCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACaa  <  1:672419/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                       CCAGGTAAGAACCGTCTTGTCCGGTTACACCGGTGATGAGAGCGACTTTTGACATGTATTATTCCTCTGTATTTTTGAATTTATTCAGTTTCAACGCGTTCGCGTATCACCACTGCGGGATT‑CCCCCGGCAAACCACATTTGCCGGAAGCGATTTAAAAACACTGCTTCGTGCACCCACGACGGTGCCGTCGCCGATTGTGACGCCTGGGGCAACAAAGACATCGGTTGCCAGCCAGCATTTCTCGCCAATCACAA  >  NC_000913/2128260‑2128515

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