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%	135.9 / NA	42	*755* (TAG→TAA)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/20);  total (22/20)

AATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGC  >  NC_000913/966448‑966718                                                                                                                                         |                                                                                                                                       aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTaaaa                                                                                                                                      >  2:151496/1‑139 (MQ=255) aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTaaaa                                                                                                                                      >  2:1125635/1‑139 (MQ=255)        tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCgg                                                                                                                               <  1:885975/139‑1 (MQ=255)        tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCgg                                                                                                                               <  2:58295/139‑1 (MQ=255)         gATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGc                                                                                                                              >  1:65897/1‑139 (MQ=255)               ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTc                                                                                                                        >  1:909691/1‑139 (MQ=255)               ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTc                                                                                                                        >  1:568220/1‑139 (MQ=255)                   aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTcaaca                                                                                                                    <  2:909691/139‑1 (MQ=255)                        ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGc                                                                                                               >  2:157365/1‑139 (MQ=255)                                 acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGtttttt                                                                                                      <  1:157365/139‑1 (MQ=255)                                   ggTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                                                                                                    >  1:962186/1‑139 (MQ=255)                                                     aGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGAt                                                                                  <  2:962186/139‑1 (MQ=255)                                                           gCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATctctct                                                                            >  2:512608/1‑139 (MQ=255)                                                               aTTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGt                                                                        <  2:65897/139‑1 (MQ=255)                                                                ttAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                              >  2:91401/1‑116 (MQ=255)                                                                ttAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                              <  1:91401/116‑1 (MQ=255)                                                                        cTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACAt                                                               >  1:839982/1‑139 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                                                                                                   >  1:625734/1‑94 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                                                                                                   <  2:625734/94‑1 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                       >  2:1079254/1‑106 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                       <  1:1079254/106‑1 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGAc                                                      >  2:704909/1‑139 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCAAGTGGATAACGGGTAAAATATGCGGc                                                                                                                              >  2:640956/1‑67 (MQ=255)                                                                                 cAAATTTTACCTCGTTGGTATCATCAGTGGTCTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGc                                                                                                                              <  1:640956/67‑1 (MQ=255)                                                                                  aaaTTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGCTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGAc                                                      >  1:694882/1‑138 (MQ=255)                                                                                  aaaTTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGCTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGAc                                                      <  2:694882/138‑1 (MQ=255)                                                                                        tACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGcc                                               <  1:1125635/139‑1 (MQ=255)                                                                                              gTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccgcc                                                         >  2:647436/1‑123 (MQ=255)                                                                                              gTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccgcc                                                         <  1:647436/123‑1 (MQ=255)                                                                                                       atcaGTGGTTTGGCGTGTCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa                                                                                            <  2:732032/79‑1 (MQ=255)                                                                                                       atcaGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa                                                                                            >  1:732032/1‑79 (MQ=255)                                                                                                       atcaGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCttt                                >  2:5359/1‑139 (MQ=255)                                                                                                         caGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTTTGGCCAACCATTGCGCCTTTCa                              >  2:100449/1‑139 (MQ=255)                                                                                                         caGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCa                              <  1:151496/139‑1 (MQ=255)                                                                                                         caGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCa                              >  1:623526/1‑139 (MQ=255)                                                                                                          aGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACtgtg                                                  >  2:400801/1‑118 (MQ=255)                                                                                                          aGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACtgtg                                                  <  1:400801/118‑1 (MQ=255)                                                                                                            tGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                           <  2:568220/139‑1 (MQ=255)                                                                                                             ggTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGc                          >  1:793384/1‑139 (MQ=255)                                                                                                                       gCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATc                <  1:704909/139‑1 (MQ=255)                                                                                                                        ccAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCg               >  1:290847/1‑139 (MQ=255)                                                                                                                                   cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTa    >  1:831851/1‑139 (MQ=255)                                                                                                                                     ggTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGc  <  1:512608/139‑1 (MQ=255)                                                                                                                                      gTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTgcg                                     <  2:99108/103‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       AATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGC  >  NC_000913/966448‑966718

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