Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,131,667	T→C	D135G (GAC→GGC)	wcaC ←	putative glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,131,667	0	T	C	100.0%	156.0 / NA	43	D135G (GAC→GGC)	wcaC	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/21);  total (22/21)

CCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGC  >  NC_000913/2131535‑2131799                                                                                                                                     |                                                                                                                                     cccGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCa                                                                                                                                <  1:621675/139‑1 (MQ=255)                  tgCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTg                                                                                                               <  1:737506/139‑1 (MQ=255)                              cTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTa                                                                                                   <  1:313634/139‑1 (MQ=255)                                   ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGtt                                                                                              >  1:108822/1‑139 (MQ=255)                                   ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGtt                                                                                              >  1:743942/1‑139 (MQ=255)                                       cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTc                                                                                          >  1:886013/1‑139 (MQ=255)                                                      ggTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAg                                                                           >  1:768260/1‑139 (MQ=255)                                                        tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg                                                                         >  1:108827/1‑139 (MQ=255)                                                        tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg                                                                         >  2:395428/1‑139 (MQ=255)                                                          gggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGca                                                                       >  2:137635/1‑139 (MQ=255)                                                           ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  2:168943/1‑139 (MQ=255)                                                           ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  1:101697/1‑139 (MQ=255)                                                           ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  1:781292/1‑139 (MQ=255)                                                           ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  1:543340/1‑139 (MQ=255)                                                             cATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                    <  2:543340/139‑1 (MQ=255)                                                               tttCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  2:1002333/139‑1 (MQ=255)                                                                   tGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGctct                                                              <  2:101697/139‑1 (MQ=255)                                                                    ggCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCtt                                                             <  1:331341/139‑1 (MQ=255)                                                                         ccTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCa                                                                                                                                >  1:51253/1‑67 (MQ=255)                                                                         ccTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCa                                                                                                                                <  2:51253/67‑1 (MQ=255)                                                                         ccTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGAt                                                        >  2:351878/1‑139 (MQ=255)                                                                          cTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGAtt                                                       <  2:18119/139‑1 (MQ=255)                                                                                  cAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGt                                               <  2:598415/139‑1 (MQ=255)                                                                                    gCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             >  2:757077/1‑139 (MQ=255)                                                                                    gCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             <  2:108827/139‑1 (MQ=255)                                                                                         cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCtgtg                                        <  1:168943/139‑1 (MQ=255)                                                                                         cGCAACCGTCGGTAAAGGCGCAGCGCCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGt                                                                                                                <  2:171323/67‑1 (MQ=255)                                                                                         cGCAACCGTCGGTAAAGGCGCAGCGCCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGt                                                                                                                >  1:171323/1‑67 (MQ=255)                                                                                            aaCCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcag                                     <  2:886013/139‑1 (MQ=255)                                                                                                 tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             >  2:148394/1‑126 (MQ=255)                                                                                                 tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             <  1:148394/126‑1 (MQ=255)                                                                                                 tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACAt                                >  1:886124/1‑139 (MQ=255)                                                                                                           cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCg                                                                                 >  1:615369/1‑80 (MQ=255)                                                                                                           cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCg                                                                                 <  2:615369/80‑1 (MQ=255)                                                                                                                         acacTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg        <  2:314871/139‑1 (MQ=255)                                                                                                                             tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  2:330048/77‑1 (MQ=255)                                                                                                                             tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  >  1:330048/1‑77 (MQ=255)                                                                                                                             tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGaa                             <  2:139319/114‑1 (MQ=255)                                                                                                                             tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGaa                             >  1:139319/1‑114 (MQ=255)                                                                                                                             tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCtgt    <  2:108822/139‑1 (MQ=255)                                                                                                                               cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  >  1:881456/1‑75 (MQ=255)                                                                                                                               cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  2:881456/75‑1 (MQ=255)                                                                                                                               cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGc  >  2:758467/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                     CCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGC  >  NC_000913/2131535‑2131799

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