Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	966,368	A→G	A683A (GCA→GCG)	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,368	0	A	G	100.0%	206.3 / NA	59	A683A (GCA→GCG)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/28);  total (31/28)

GGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACT  >  NC_000913/966271‑966461                                                                                                  |                                                                                              ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCt                                                                                             <  2:1246718/100‑1 (MQ=255)   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                           <  2:1158143/100‑1 (MQ=255)   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                           <  1:127636/100‑1 (MQ=255)   cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                           >  2:458604/1‑100 (MQ=255)        aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGACGCTCgcatcg                                                                                      >  2:1236817/1‑100 (MQ=255)         gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                     >  1:917924/1‑100 (MQ=255)          tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                    >  2:1427512/1‑100 (MQ=255)           gCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggca                                                                                   >  2:595766/1‑100 (MQ=255)            cAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggca                                                                                   <  2:1457847/99‑1 (MQ=255)              ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatcg                                                                                >  1:1369769/1‑100 (MQ=255)                  aCATTGCTTGAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGGTCGCATCGGCATCGCGCt                                                                            <  2:13134/100‑1 (MQ=255)                      tGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATaa                                                                        >  1:1243450/1‑100 (MQ=255)                       gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAc                                                                       >  1:1272749/1‑100 (MQ=255)                        cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                      >  1:421832/1‑100 (MQ=255)                            aGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATg                                                                  >  1:1473514/1‑100 (MQ=255)                             ggTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATgg                                                                 >  1:1345594/1‑100 (MQ=255)                               tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAACGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCg                                                               >  1:545171/1‑100 (MQ=255)                                 ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                             >  1:185164/1‑100 (MQ=255)                                  cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                            >  1:1274947/1‑100 (MQ=255)                                  cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                            >  1:708776/1‑100 (MQ=255)                                      ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGGGCTCACATCGGCATCGCGCTATAACGCATGGCGACTGc                                                          >  2:1200050/1‑98 (MQ=255)                                      ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                        >  2:1466095/1‑100 (MQ=255)                                         tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                     >  1:161290/1‑100 (MQ=255)                                         tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                     <  2:1524747/100‑1 (MQ=255)                                         tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                     <  2:167777/100‑1 (MQ=255)                                          ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                      <  2:545170/98‑1 (MQ=255)                                          ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                    <  2:478905/100‑1 (MQ=255)                                          ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                    <  1:506805/100‑1 (MQ=255)                                            cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                  >  2:188696/1‑100 (MQ=255)                                                 aTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGt                                             <  2:1455753/100‑1 (MQ=255)                                                   aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                           <  2:1121295/100‑1 (MQ=255)                                                     cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                         <  1:96347/100‑1 (MQ=255)                                                      tcgtcagcGTTGCCATGAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCAGCGCGCTATAACGCATGGCGACTGCCCTCTAACaaa                                               <  2:690222/90‑1 (MQ=255)                                                       catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                       >  2:260074/1‑100 (MQ=255)                                                        atcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCAt                                      >  2:631161/1‑100 (MQ=255)                                                         tcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATc                                     <  2:1368165/100‑1 (MQ=255)                                                          catcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGTCATCGCGCTATAACGCATGGCGTCTGCCCTCTAACAAAGTTAAGCATCg                                    <  2:753431/100‑1 (MQ=255)                                                           atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                   <  1:1410265/100‑1 (MQ=255)                                                             cGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTa                                 >  1:1145433/1‑100 (MQ=255)                                                               ttGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATc                               >  1:375224/1‑100 (MQ=255)                                                                  ccATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAAc                            <  2:1243445/100‑1 (MQ=255)                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACCAAGTTAAGCATCGCTATCAACt                           <  2:917919/100‑1 (MQ=255)                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                           >  2:686062/1‑100 (MQ=255)                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                           >  2:1203711/1‑100 (MQ=255)                                                                   cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                           >  1:890609/1‑100 (MQ=255)                                                                      tAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCa                        <  1:458605/100‑1 (MQ=255)                                                                       aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                       >  2:81238/1‑100 (MQ=255)                                                                       aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                       >  2:1469598/1‑100 (MQ=255)                                                                        aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                      <  2:375224/100‑1 (MQ=255)                                                                        aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                      >  1:125503/1‑100 (MQ=255)                                                                        aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                      <  1:559902/100‑1 (MQ=255)                                                                        aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                      <  1:722807/100‑1 (MQ=255)                                                                         ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                     <  1:387014/100‑1 (MQ=255)                                                                              cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                >  2:590988/1‑100 (MQ=255)                                                                               gAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa               <  2:185164/100‑1 (MQ=255)                                                                                 gTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAAt             <  2:7237/100‑1 (MQ=255)                                                                                          aaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATa    <  2:1145428/100‑1 (MQ=255)                                                                                           aaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATAc   <  1:1571928/100‑1 (MQ=255)                                                                                            aGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACt  >  1:317574/1‑100 (MQ=255)                                                                                                  |                                                                                              GGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACT  >  NC_000913/966271‑966461

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