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%	172.1 / NA	49	D135G (GAC→GGC)	wcaC	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/25);  total (24/25)

TGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGC  >  NC_000913/2131533‑2131799                                                                                                                                       |                                                                                                                                     tGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTg                                                                                                                                  >  2:1738422/1‑139 (MQ=255)           caGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTccag                                                                                                                        <  2:647480/139‑1 (MQ=255)                    tgCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTg                                                                                                               >  2:808314/1‑139 (MQ=255)                     gcgcGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTcca                                                                                                                         <  1:1550245/128‑1 (MQ=255)                     gcgcGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTcca                                                                                                                         >  2:1550245/1‑128 (MQ=255)                        cgATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                           <  1:120777/139‑1 (MQ=255)                           tCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCg                                                                                                        <  1:1633731/139‑1 (MQ=255)                                     ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGtt                                                                                              >  1:50847/1‑139 (MQ=255)                                         cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTc                                                                                          >  1:1121632/1‑139 (MQ=255)                                         cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTc                                                                                          >  2:1529355/1‑139 (MQ=255)                                                     tAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCaa                                                                              >  1:629435/1‑139 (MQ=255)                                                        ggTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAg                                                                           <  2:629435/139‑1 (MQ=255)                                                        ggTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAg                                                                           <  2:1121632/139‑1 (MQ=255)                                                          tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg                                                                         >  2:785120/1‑139 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCATCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTccagac                                                                                                                      <  1:157240/91‑1 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCa                                                                                                                                >  2:370677/1‑81 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCa                                                                                                                                <  1:370677/81‑1 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTccagac                                                                                                                      >  2:157240/1‑91 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  1:1105384/1‑139 (MQ=255)                                                             ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                      >  1:97220/1‑139 (MQ=255)                                                                 tttCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  2:97220/139‑1 (MQ=255)                                                                 tttCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  2:951900/139‑1 (MQ=255)                                                                 tttCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                  <  1:1325478/139‑1 (MQ=255)                                                                            cTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGAtt                                                       <  2:1105384/139‑1 (MQ=255)                                                                               ttttCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACAGTCCAGTGGCCGTGCAGCGTcca                                                                                                                         <  2:265760/70‑1 (MQ=255)                                                                               ttttCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTcca                                                                                                                         >  1:265760/1‑70 (MQ=255)                                                                                   ccAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAg                                                <  1:1188937/139‑1 (MQ=255)                                                                                      gCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGCGTGACATCCGGTTTATGGTTATTCACTTTTTCGCAACAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             <  1:1738422/139‑1 (MQ=255)                                                                                      gCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             >  2:468731/1‑139 (MQ=255)                                                                                       ccTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGc                                            >  2:1533149/1‑139 (MQ=255)                                                                                           cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                    >  2:36319/1‑111 (MQ=255)                                                                                           cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                    <  1:36319/111‑1 (MQ=255)                                                                                           cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCtgtg                                        <  1:808314/139‑1 (MQ=255)                                                                                           cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCtgtg                                        >  2:1683759/1‑139 (MQ=255)                                                                                                   tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcagca                                   >  2:166342/1‑136 (MQ=255)                                                                                                   tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcagca                                   <  1:166342/136‑1 (MQ=255)                                                                                                     ggTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                           >  1:366210/1‑62 (MQ=255)                                                                                                     ggTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                           <  2:366210/62‑1 (MQ=255)                                                                                                            gcgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCg                                                                                 >  1:181071/1‑81 (MQ=255)                                                                                                            gcgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCg                                                                                 <  2:181071/81‑1 (MQ=255)                                                                                                              gcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAgg                     <  2:335533/139‑1 (MQ=255)                                                                                                                              cTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCtg     >  2:145800/1‑139 (MQ=255)                                                                                                                                 cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             >  2:1677786/1‑96 (MQ=255)                                                                                                                                 cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAg                                             <  1:1677786/96‑1 (MQ=255)                                                                                                                                 cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGc  <  1:1683759/139‑1 (MQ=255)                                                                                                                                 cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGc  >  2:1280998/1‑139 (MQ=255)                                                                                                                                   gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac                                                                   <  1:1208411/72‑1 (MQ=255)                                                                                                                                   gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac                                                                   >  2:1208411/1‑72 (MQ=255)                                                                                                                                    tGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCt      <  2:1495988/132‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     TGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGC  >  NC_000913/2131533‑2131799

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