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

ATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCG  >  NC_000913/966453‑966719                                                                                                                                    |                                                                                                                                        aTTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGc                                                                                                                                  >  2:548111/1‑139 (MQ=255)  ttGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCg                                                                                                                                 <  2:907204/139‑1 (MQ=255)     aTACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCt                                                                                                                              >  2:642763/1‑139 (MQ=255)     aTACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCt                                                                                                                              >  1:224466/1‑139 (MQ=255)          ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCg                                                                                                                                 <  1:193176/131‑1 (MQ=255)          ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCg                                                                                                                                 >  2:193176/1‑131 (MQ=255)           cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTca                                                                                                                        >  2:197167/1‑139 (MQ=255)           cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTca                                                                                                                        >  2:1425205/1‑139 (MQ=255)            acaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTcaa                                                                                                                       <  1:642763/139‑1 (MQ=255)            acaTCAAGGTCAAACAACGGTCAATTTTTAAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTcaa                                                                                                                       >  1:191899/1‑139 (MQ=255)             caTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTcatc                                                                                                                      >  2:403178/1‑137 (MQ=255)              aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTcaaca                                                                                                                     <  2:1422616/139‑1 (MQ=255)                cAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACaaa                                                                                                                   >  1:680987/1‑139 (MQ=255)                        aacaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGtt                                                                                                           >  2:881435/1‑139 (MQ=255)                                    tttttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt                                                                                               <  1:548111/139‑1 (MQ=255)                                     ttttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATa                                                                                              <  1:197167/139‑1 (MQ=255)                                          aGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc                                                                                         <  1:1425205/139‑1 (MQ=255)                                             gcAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa                                                                                      <  2:1280430/139‑1 (MQ=255)                                              cAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAg                                                                                     <  2:1150943/139‑1 (MQ=255)                                              cAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAg                                                                                     >  2:457135/1‑139 (MQ=255)                                                                   cTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACAt                                                                <  1:881435/139‑1 (MQ=255)                                                                      aGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTca                                                                                                                        >  2:948541/1‑80 (MQ=255)                                                                      aGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTca                                                                                                                        <  1:948541/80‑1 (MQ=255)                                                                        ggAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccgc                                                           <  2:680987/139‑1 (MQ=255)                                                                            cAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGAc                                                       >  1:1230129/1‑139 (MQ=255)                                                                                  ttACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa                                                                                             <  1:711544/95‑1 (MQ=255)                                                                                  ttACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATaa                                                                                             >  2:711544/1‑95 (MQ=255)                                                                                         gTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCAt                                          >  1:481819/1‑139 (MQ=255)                                                                                         gTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCAt                                          >  2:796154/1‑139 (MQ=255)                                                                                          ttGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccg                                                            <  2:1224278/120‑1 (MQ=255)                                                                                          ttGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTccg                                                            >  1:1224278/1‑120 (MQ=255)                                                                                                   tcaGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                        <  1:517620/83‑1 (MQ=255)                                                                                                   tcaGTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                        >  2:517620/1‑83 (MQ=255)                                                                                                      gTGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCaaa                             <  2:481819/139‑1 (MQ=255)                                                                                                       tGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                            >  1:693606/1‑139 (MQ=255)                                                                                                       tGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                            >  1:1307398/1‑139 (MQ=255)                                                                                                       tGGTTTGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                            <  1:457135/139‑1 (MQ=255)                                                                                                          tttGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCa                                                                                                >  2:340074/1‑68 (MQ=255)                                                                                                          tttGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCa                                                                                                <  1:340074/68‑1 (MQ=255)                                                                                                          tttGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCgg                         >  1:149541/1‑139 (MQ=255)                                                                                                            tGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGt                                                                         <  1:1401986/89‑1 (MQ=255)                                                                                                            tGGCGTGCCAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGt                                                                         >  2:1401986/1‑89 (MQ=255)                                                                                                                   ccAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTgg                                                  >  1:171675/1‑105 (MQ=255)                                                                                                                   ccAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTgg                                                  <  2:171675/105‑1 (MQ=255)                                                                                                                   ccAGTGGATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCg                >  1:1202758/1‑139 (MQ=255)                                                                                                                        ggATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGT‑‑CAGACATTCCGCCGACtgtg                                                   <  2:680843/97‑1 (MQ=255)                                                                                                                        ggATAACGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGT‑‑CAGACATTCCGCCGACtgtg                                                   >  1:680843/1‑97 (MQ=255)                                                                                                                              cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                        >  1:141721/1‑56 (MQ=255)                                                                                                                              cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa                                                                                        <  2:141721/56‑1 (MQ=255)                                                                                                                              cGGGTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTa     >  1:359599/1‑139 (MQ=255)                                                                                                                                ggTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGc   >  2:633005/1‑139 (MQ=255)                                                                                                                                 gTAAAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCg  >  2:390236/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        ATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCG  >  NC_000913/966453‑966719

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