Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	966,584	G→A	*755* (TAG→TAA)	ycaI →	ComEC family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	966,584	0	G	A	100.0%	54.3 / NA	19	*755* (TAG→TAA)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/9);  total (10/9)

TGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTA  >  NC_000913/966455‑966716                                                                                                                                  |                                                                                                                                     tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCgg                                                                                                                             <  2:776663/139‑1 (MQ=255)   aTACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCt                                                                                                                           >  1:746006/1‑139 (MQ=255)      cTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTAtt                                                                                                                        <  1:796215/139‑1 (MQ=255)                            ggTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                                                                                                  >  1:363547/1‑139 (MQ=255)                             gTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                                                                                                 <  2:475577/139‑1 (MQ=255)                                      tCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGt                                                                                                         >  1:332126/1‑122 (MQ=255)                                      tCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGt                                                                                                         <  2:332126/122‑1 (MQ=255)                                         gcgcAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                     <  2:827651/139‑1 (MQ=255)                                                    gCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATctctct                                                                          <  2:633262/139‑1 (MQ=255)                                                                  tCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACAtt                                                            >  2:227810/1‑139 (MQ=255)                                                                                                     tGGTTTGGCGTGCCAGTGGATAACGGGTAAGATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGa                                                                                 >  2:275237/1‑83 (MQ=255)                                                                                                     tGGTTTGGCGTGCCAGTGGATAACGGGTAAGATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGa                                                                                 <  1:275237/83‑1 (MQ=255)                                                                                                               tGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa                                                                                   >  1:79535/1‑71 (MQ=255)                                                                                                               tGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa                                                                                   <  2:79535/71‑1 (MQ=255)                                                                                                               tGCCAGTGGATAACGAGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAAGGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGAt               >  2:560587/1‑139 (MQ=255)                                                                                                                          aaCGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg    >  1:798880/1‑139 (MQ=255)                                                                                                                            cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                     <  2:302580/56‑1 (MQ=255)                                                                                                                            cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                     >  1:302580/1‑56 (MQ=255)                                                                                                                            cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTa  >  2:142412/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                     TGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTA  >  NC_000913/966455‑966716

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