Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	965,917	A→G	D533G (GAT→GGT)	ycaI →	ComEC family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	965,917	0	A	G	100.0%	178.2 / NA	50	D533G (GAT→GGT)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

GCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCATATATGGCCGATGTTA  >  NC_000913/965783‑966047                                                                                                                                       |                                                                                                                                   gCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGt                                                                                                                                >  2:1143831/1‑139 (MQ=255)   cGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTgg                                                                                                                              >  1:253301/1‑139 (MQ=255)    gCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGACCTGAAGGGGGTAGTggg                                                                                                                             <  1:1401087/139‑1 (MQ=255)        tgGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGcaac                                                                                                                         <  2:172499/139‑1 (MQ=255)                aCCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGAtt                                                                                                                 <  1:1143831/139‑1 (MQ=255)                         tGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTgg                                                                                                        <  1:194788/139‑1 (MQ=255)                                                ttGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCAGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAAcc                                                                                 >  2:208948/1‑139 (MQ=255)                                                   aTGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTgg                                                                                                        <  2:958744/113‑1 (MQ=255)                                                   aTGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTgg                                                                                                        >  1:958744/1‑113 (MQ=255)                                                   aTGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGa                                                                              >  2:1531275/1‑139 (MQ=255)                                                                   ggCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGc                                                              >  2:967006/1‑139 (MQ=255)                                                                    gCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCa                                                                                             <  1:1329027/107‑1 (MQ=255)                                                                    gCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCa                                                                                             >  2:1329027/1‑107 (MQ=255)                                                                        tGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCAGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGa                                                         <  1:208948/139‑1 (MQ=255)                                                                         ggTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGaa                                                        >  1:1346626/1‑139 (MQ=255)                                                                          gTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCg                                                                         >  2:397996/1‑121 (MQ=255)                                                                          gTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCg                                                                         <  1:397996/121‑1 (MQ=255)                                                                          gTGATAACCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAAc                                                       >  1:988060/1‑139 (MQ=255)                                                                              tAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCt                                                   >  2:170482/1‑139 (MQ=255)                                                                                  cAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATg                                                          <  2:228223/128‑1 (MQ=255)                                                                                  cAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATg                                                          >  1:228223/1‑128 (MQ=255)                                                                                  cAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGAt                                               <  2:1346626/139‑1 (MQ=255)                                                                                          gCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCgggg                                       >  1:172481/1‑139 (MQ=255)                                                                                           cAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACAt                                                     >  1:912974/1‑124 (MQ=255)                                                                                           cAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACAt                                                     <  2:912974/124‑1 (MQ=255)                                                                                              aGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAggg                                   <  1:846611/139‑1 (MQ=255)                                                                                              aGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAggg                                   <  1:393370/139‑1 (MQ=255)                                                                                              aGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAggg                                   <  1:967006/139‑1 (MQ=255)                                                                                                cGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCt                                 <  1:1531275/139‑1 (MQ=255)                                                                                                        tATGACACAGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCaa                         >  1:1311924/1‑139 (MQ=255)                                                                                                                aGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGa                                                                              >  1:1351674/1‑78 (MQ=255)                                                                                                                aGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGa                                                                              <  2:1351674/78‑1 (MQ=255)                                                                                                                aGGACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCat                 >  1:819760/1‑139 (MQ=255)                                                                                                                  gACTGGCCTGGCCTGAAGGGGGTAGTGGGCAACAACTGATTATCCTCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACACCTGGATCACCGGGGAGGGCTGGATTAAATATTGCgtat               <  2:988060/139‑1 (MQ=255)                                                                                                                    ctggcctggcCTGAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCatatat             >  2:416662/1‑139 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGa                                                                              >  2:1062234/1‑71 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGa                                                                              <  1:1062234/71‑1 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAggg                                   <  1:376584/114‑1 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAggg                                   >  2:376584/1‑114 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCATATATGGc          >  2:453102/1‑139 (MQ=255)                                                                                                                       gcctggcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCATATATGGc          <  2:819760/139‑1 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                                   <  1:243224/45‑1 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCa                                                                                                   >  2:243224/1‑45 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAg                                                                            >  1:969743/1‑68 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAg                                                                            <  2:969743/68‑1 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGAtt                            <  2:50907/116‑1 (MQ=255)                                                                                                                            gcctgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGAtt                            >  1:50907/1‑116 (MQ=255)                                                                                                                               tgAAGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCATATATGGCCGATGTTa  <  2:172481/139‑1 (MQ=255)                                                                                                                                 aaGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAg                                                               <  2:1496205/76‑1 (MQ=255)                                                                                                                                 aaGGGGGTAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAg                                                               >  1:1496205/1‑76 (MQ=255)                                                                                                                                       |                                                                                                                                   GCCGCTGTGGCAAAAACCTCGACCTGACGAGTGGCAGCTGTACATGCTTGATGTCGGGCAAGGGCTGGCAATGGTGATAGCCAGAAACGGCAAAGCGATTCTCTATGACACAGGACTGGCCTGGCCTGAAGGGGATAGTGGGCAACAACTGATTATCCCCTGGCTCCACTGGCATAATCTTGAACCGGAAGGCGTTATTCTGAGCCATGAACATCTGGATCACCGGGGAGGGCTGGATTCAATATTGCATATATGGCCGATGTTA  >  NC_000913/965783‑966047

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: ‑