Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	972,866	Δ1 bp	coding (536/780 nt)	elyC ←	envelope biogenesis factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	972,863	0	G	.	100.0%	77.8 / NA	19	coding (539/780 nt)	elyC	envelope biogenesis factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (10/9);  total (10/9)

TCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTT  >  NC_000913/972819‑972905                                             |                                           tcatcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAgggcg                                        >  1:2167510/1‑47 (MQ=255)    tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCa                                    <  1:858/49‑1 (MQ=39)    tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCa                                    <  1:381735/49‑1 (MQ=39)    tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCa                                    <  1:330008/49‑1 (MQ=39)      gCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCAc                                   >  1:2805513/1‑48 (MQ=39)      gCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCAc                                   >  1:2697579/1‑48 (MQ=39)       cccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCa                                 <  1:828202/49‑1 (MQ=39)           cgGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCg                             <  1:137943/49‑1 (MQ=39)            ggCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGc                            >  1:1344040/1‑49 (MQ=39)                 aTGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGtt                       >  1:1884858/1‑49 (MQ=39)                    tgATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTc                     <  1:2615949/48‑1 (MQ=39)                      aTGTGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCAct                  >  1:1559683/1‑49 (MQ=38)                      aTGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCAct                  >  1:1559694/1‑49 (MQ=39)                      aTGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCAc                   <  1:1827642/48‑1 (MQ=39)                         cGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgct               <  1:1285505/49‑1 (MQ=39)                          ggAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgctg              >  1:763611/1‑49 (MQ=39)                           gAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgctg              >  1:1725103/1‑48 (MQ=39)                            aGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgctg              >  1:2306291/1‑47 (MQ=39)                               tCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGct         <  1:2694272/49‑1 (MQ=39)                                      cagGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGcttcttctt  >  1:2634642/1‑49 (MQ=39)                                             |                                           TCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTT  >  NC_000913/972819‑972905

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