Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,389,242	C→T	Q105* (CAA→TAA)	ycjG →	L‑Ala‑D/L‑Glu epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,389,242	0	C	T	100.0%	127.4 / NA	41	Q105* (CAA→TAA)	ycjG	L‑Ala‑D/L‑Glu epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

TCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCG  >  NC_000913/1389123‑1389364                                                                                                                        |                                                                                                                           tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                         >  2:842707/1‑139 (MQ=255) tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                         >  1:1477979/1‑139 (MQ=255) tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                         >  1:95924/1‑139 (MQ=255)    tGCCGCAACTAGAGAAAGGGCTTACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAAGGCGCTGTATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCgg                                                                                                      >  1:830452/1‑139 (MQ=255)           aCTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACg                                                                                               <  2:830452/139‑1 (MQ=255)            cTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTCGATTGTGCAATGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGc                                                                                              >  2:406355/1‑139 (MQ=255)                    aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                      <  1:1196788/139‑1 (MQ=255)                    aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                      <  2:1143333/139‑1 (MQ=255)                    aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                      <  2:1414122/139‑1 (MQ=255)                         tGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGt                                                                                 <  1:1483285/139‑1 (MQ=255)                            cacGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                              <  2:735103/139‑1 (MQ=255)                                     agTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcac                                                                      <  2:82377/138‑1 (MQ=255)                                     agTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcac                                                                      >  1:82377/1‑138 (MQ=255)                                       ttGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACg                                                                                               <  2:1409710/111‑1 (MQ=255)                                       ttGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACg                                                                                               >  1:1409710/1‑111 (MQ=255)                                               aaTTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTc                                                           <  2:1477979/139‑1 (MQ=255)                                               aaTTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTc                                                           <  2:95924/139‑1 (MQ=255)                                                    tcCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCgg                                                      <  1:842707/139‑1 (MQ=255)                                                    tcCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCgg                                                      <  1:17104/139‑1 (MQ=255)                                                                cagcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgaga                                                                                     <  2:1131277/96‑1 (MQ=255)                                                                cagcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCCTCCCgaga                                                                                     >  1:1131277/1‑96 (MQ=255)                                                                   caCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGAt                                       >  1:752793/1‑139 (MQ=255)                                                                   caCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGATATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGAt                                       >  1:237853/1‑139 (MQ=255)                                                                              cTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTa                                                    >  1:157766/1‑115 (MQ=255)                                                                              cTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTa                                                    <  2:157766/115‑1 (MQ=255)                                                                                ggATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCa                          >  2:698592/1‑139 (MQ=255)                                                                                 gATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                         <  1:385643/59‑1 (MQ=255)                                                                                 gATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                         >  2:385643/1‑59 (MQ=255)                                                                                      tgCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAaca                    >  1:532192/1‑139 (MQ=255)                                                                                          ttGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGc                                                 <  2:379625/106‑1 (MQ=255)                                                                                          ttGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGc                                                 >  1:379625/1‑106 (MQ=255)                                                                                           tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGc                           >  2:1268270/1‑127 (MQ=255)                                                                                           tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGc                           <  1:1268270/127‑1 (MQ=255)                                                                                             tgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTggc             >  1:204501/1‑139 (MQ=255)                                                                                                gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTggcagg          >  2:663776/1‑139 (MQ=255)                                                                                                       ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                              >  2:1344775/1‑64 (MQ=255)                                                                                                       ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                              <  1:1344775/64‑1 (MQ=255)                                                                                                       ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGgcgc   >  2:1489721/1‑139 (MQ=255)                                                                                                        cgcgcgAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGcgcg  <  1:1528266/139‑1 (MQ=255)                                                                                                          cgcgAAGGCAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTAc                                                   <  2:1525300/88‑1 (MQ=255)                                                                                                          cgcgAAGGCAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTAc                                                   >  1:1525300/1‑88 (MQ=255)                                                                                                                     cagtaATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                     <  2:711026/59‑1 (MQ=255)                                                                                                                        |                                                                                                                           TCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCG  >  NC_000913/1389123‑1389364

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