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%	90.7 / NA	22	coding (539/780 nt)	elyC	envelope biogenesis factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (14/8);  total (14/8)

AAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGG  >  NC_000913/972814‑972908                                                  |                                              aaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAg                                               <  1:719126/50‑1 (MQ=255)  aaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAgg                                              <  1:1078857/50‑1 (MQ=255)  aaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAgg                                              <  1:1269980/50‑1 (MQ=255)  aaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAgg                                              <  1:1375204/50‑1 (MQ=255)         tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCAc                                      <  1:1079560/50‑1 (MQ=39)           gCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCAcc                                     >  1:1921070/1‑49 (MQ=39)             ccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAAt                                  >  1:1692549/1‑50 (MQ=39)             ccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAAt                                  >  1:1362173/1‑50 (MQ=39)             ccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAAt                                  >  1:1492089/1‑50 (MQ=39)                 ggCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGcc                              >  1:848335/1‑50 (MQ=39)                 ggCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGcc                              >  1:1405808/1‑50 (MQ=39)                 ggCAGATGTGATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGcc                              <  1:1344215/50‑1 (MQ=39)                       tgtgATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGttt                         >  1:788938/1‑49 (MQ=39)                        gtgATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCa                       >  1:771415/1‑50 (MQ=39)                         tgATGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCAc                      <  1:2247621/50‑1 (MQ=39)                            tGCGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgc                   <  1:604686/50‑1 (MQ=39)                              cGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgctg                 >  1:1334012/1‑50 (MQ=39)                              cGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgctg                 >  1:2001712/1‑50 (MQ=39)                              cGGAGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCActgct                  <  1:1076330/49‑1 (MQ=39)                                 aGGTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAg              <  1:1584329/50‑1 (MQ=39)                                   gTCACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGct            <  1:144984/50‑1 (MQ=39)                                     cACCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGcttc          >  1:1574331/1‑50 (MQ=39)                                      aCCAGCAGGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGcttct         >  1:1717592/1‑50 (MQ=39)                                         agcagGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGcttcttct      <  1:401452/50‑1 (MQ=39)                                           cagGAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCttcttcttc    >  1:1628703/1‑50 (MQ=39)                                             ggAAA‑GGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCgg  >  1:1037491/1‑50 (MQ=39)                                                  |                                              AAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGG  >  NC_000913/972814‑972908

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