Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	964,634	T→C	R105R (CGT→CGC)	ycaI →	ComEC family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	964,634	0	T	C	100.0%	150.7 / NA	40	R105R (CGT→CGC)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/22);  total (18/22)

GTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGTATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTAT  >  NC_000913/964491‑964775                                                                                                                                                |                                                                                                                                              gTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTcc                                                                                                                                        >  1:501155/1‑151 (MQ=255)   ttGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCt                                                                                                                                      >  2:125901/1‑151 (MQ=255)    tGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTg                                                                                                                                     >  1:51058/1‑151 (MQ=255)     ggggCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTg                                                                                                                                     <  2:160635/150‑1 (MQ=255)      gggCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCg                                                                                                                                   <  2:197380/151‑1 (MQ=255)                   gcaaagcaaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGa                                                                                                                      <  1:110751/151‑1 (MQ=255)                         caaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATg                                                                                                                >  1:19998/1‑151 (MQ=255)                          aaaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATgg                                                                                                               <  2:187625/151‑1 (MQ=255)                           aaTTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGc                                                                                                              >  2:191798/1‑151 (MQ=255)                              tttGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGa                                                                                                            >  1:632830/1‑150 (MQ=255)                                        ggAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGc                                                                                                 <  2:170383/151‑1 (MQ=255)                                         gagaAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCa                                                                                                <  1:281175/151‑1 (MQ=255)                                            aaaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGc                                                                                             <  2:341108/151‑1 (MQ=255)                                             aaCCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTATAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCg                                                                                            >  1:277216/1‑151 (MQ=255)                                                   aaCTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTtgt                                                                                      >  2:537029/1‑151 (MQ=255)                                                           cgACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGaca                                                                              <  2:400308/151‑1 (MQ=255)                                                             aCGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGacaac                                                                            >  1:279064/1‑151 (MQ=255)                                                               gCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGacaaca                                                                           >  2:86845/1‑150 (MQ=255)                                                                cAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAAt                                                                         >  2:418741/1‑151 (MQ=255)                                                                  ggATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAg                                                                                              <  1:614593/128‑1 (MQ=255)                                                                  ggATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAg                                                                                              >  2:614576/1‑128 (MQ=255)                                                                      gCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCa                                                                    <  1:22744/150‑1 (MQ=255)                                                                                    cACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGtt                                                     <  1:412420/151‑1 (MQ=255)                                                                                     aCTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTc                                                    <  1:305903/151‑1 (MQ=255)                                                                                        gCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCg                                                   <  1:499674/149‑1 (MQ=255)                                                                                          aaCTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGc                                                >  2:67726/1‑150 (MQ=255)                                                                                                cGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCAt                                         >  1:576671/1‑151 (MQ=255)                                                                                                cGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCAt                                         >  1:322853/1‑151 (MQ=255)                                                                                                      gACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCa                                    <  1:67726/150‑1 (MQ=255)                                                                                                             cATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAatg                            <  2:197416/151‑1 (MQ=255)                                                                                                             cATTACGGTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAat                             <  2:54728/150‑1 (MQ=255)                                                                                                                   ggTCAAATTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATggcg                      <  1:191803/151‑1 (MQ=255)                                                                                                                     tCAAATTACGCATCTACAAGGTCGACGCATATTCCGTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATggcg                      <  2:95888/149‑1 (MQ=255)                                                                                                                       aaaTTACTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCtt                  <  2:109413/151‑1 (MQ=255)                                                                                                                            aCTCATCTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATa             <  1:215855/151‑1 (MQ=255)                                                                                                                                  cTACAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGc       >  1:535468/1‑151 (MQ=255)                                                                                                                                    aCAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCg      <  1:175747/150‑1 (MQ=255)                                                                                                                                    aCAAGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGt     >  1:266692/1‑151 (MQ=255)                                                                                                                                       aGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTAt  >  2:359664/1‑151 (MQ=255)                                                                                                                                       aGGTCGACGCATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTAGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTAt  <  1:480805/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                              GTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAACCTTAACTGGCGCGACGCAGGATGCAATTGTTGAGATCACTGCAACTGACGGCATGACCACTCATTACGGTCAAATTACTCATCTACAAGGTCGACGTATATTCCCTGCGTCAGGTCTCGTGATGTATGGCGAATATCTTCCGCAAGCGGTTTGTGCTGGACAACAATGGTCAATGAAACTCAAAGTTCGTGCAGTTCATGGCCAACTTAATGATGGCGGCTTTGATAGCCAGCGTTAT  >  NC_000913/964491‑964775

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