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

AATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCG  >  NC_000913/966448‑966714                                                                                                                                          |                                                                                                                                   aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTaaaa                                                                                                                                  >  1:710376/1‑139 (MQ=255)      aTTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGc                                                                                                                             >  2:170636/1‑139 (MQ=255)          aTACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCt                                                                                                                         <  2:397778/139‑1 (MQ=255)                cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAata                                                                                                                                >  2:1466500/1‑126 (MQ=255)                cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAata                                                                                                                                <  1:1466500/126‑1 (MQ=255)                cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTca                                                                                                                   >  1:1167005/1‑139 (MQ=255)                 acaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTcaa                                                                                                                  >  1:210115/1‑139 (MQ=255)                   aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTcaaca                                                                                                                <  2:761871/139‑1 (MQ=255)                    tCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCg                                                                                                                            <  2:581509/126‑1 (MQ=255)                    tCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCg                                                                                                                            >  1:581509/1‑126 (MQ=255)                     cAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACaaa                                                                                                              >  2:428906/1‑139 (MQ=255)                     cAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACaaa                                                                                                              >  2:161689/1‑139 (MQ=255)                     cAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACaaa                                                                                                              >  1:446003/1‑139 (MQ=255)                       aGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATg                                                                                                            <  2:1108178/139‑1 (MQ=255)                          tCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTg                                                                                                         >  1:951383/1‑139 (MQ=255)                             aacaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGtt                                                                                                      >  1:411633/1‑139 (MQ=255)                              acaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGttt                                                                                                     <  2:477424/139‑1 (MQ=255)                                   ggTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                                                                                                >  2:1329777/1‑139 (MQ=255)                                    gTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                                                                                               >  2:343938/1‑139 (MQ=255)                                    gTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGaa                                                                                               <  2:605411/139‑1 (MQ=255)                                        aTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCa                                                                                           <  2:463449/139‑1 (MQ=255)                                         tttttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                          <  1:1312848/139‑1 (MQ=255)                                               aGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc                                                                                    <  1:428906/139‑1 (MQ=255)                                               aGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc                                                                                    >  2:540252/1‑139 (MQ=255)                                                gcgcAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                   <  2:951383/139‑1 (MQ=255)                                                             ggATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTAc                                                                      <  1:540252/139‑1 (MQ=255)                                                                      gcCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGAc                                                             <  1:518052/139‑1 (MQ=255)                                                                         tCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACAtt                                                          >  2:612351/1‑139 (MQ=255)                                                                         tCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACAtt                                                          >  2:26458/1‑139 (MQ=255)                                                                              gAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccgc                                                      <  1:1216138/139‑1 (MQ=255)                                                                                  aaaTTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACt                                                 <  2:1167005/139‑1 (MQ=255)                                                                                        tACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAATTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGcc                                           <  2:210115/139‑1 (MQ=255)                                                                                        tACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGcc                                           <  2:474725/139‑1 (MQ=255)                                                                                               ttGGTATCATCAGTGGTTTGCCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTctgc                                                      >  2:902597/1‑118 (MQ=255)                                                                                               ttGATATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTctgc                                                      <  1:902597/121‑4 (MQ=255)                                                                                                        tcaGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACg                                                                                      >  1:380859/1‑80 (MQ=255)                                                                                                        tcaGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACg                                                                                      <  2:380859/80‑1 (MQ=255)                                                                                                          aGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                          >  1:371610/1‑74 (MQ=255)                                                                                                          aGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                          <  2:371610/74‑1 (MQ=255)                                                                                                          aGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa                                                                                 <  1:892132/83‑1 (MQ=255)                                                                                                          aGTGATTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa                                                                                 >  2:892132/1‑83 (MQ=255)                                                                                                           gTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCaaa                        <  1:26458/139‑1 (MQ=255)                                                                                                            tGGTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATACGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                       >  1:113333/1‑139 (MQ=255)                                                                                                              gTTTGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCg                     <  1:612351/139‑1 (MQ=255)                                                                                                                 tGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGt                  >  2:445233/1‑139 (MQ=255)                                                                                                                 tGGCGTGCCAGTGGATAAC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGt                  >  2:576757/1‑139 (MQ=255)                                                                                                                                 aaC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg  <  1:161689/139‑1 (MQ=255)                                                                                                                                 aaC‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg  >  1:904049/1‑139 (MQ=255)                                                                                                                                   c‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGAtct                                                                            >  1:1328373/1‑63 (MQ=255)                                                                                                                                   c‑GGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGAtct                                                                            <  2:1328373/63‑1 (MQ=255)                                                                                                                                       gTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCa                          >  1:305265/1‑110 (MQ=255)                                                                                                                                       gTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCa                          <  2:305265/110‑1 (MQ=255)                                                                                                                                       gTAAAATATGCGCCTATTTCAACAAATGCTGGTTTTTTGAATGCATa                                                                                         <  1:1261131/47‑1 (MQ=39)                                                                                                                                       gTAAAATATGCGCCTATTTCAACAAATGCTGGTTTTTTGAATGCATa                                                                                         >  2:1261131/1‑47 (MQ=38)                                                                                                                                          |                                                                                                                                   AATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAC‑GGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCG  >  NC_000913/966448‑966714

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