Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,598,617	C→T	*531* (TAG→TAA)	lsrK ←	autoinducer‑2 (AI‑2) kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,598,617	0	C	T	100.0%	116.9 / NA	37	*531* (TAG→TAA)	lsrK	autoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/20);  total (17/20)

CGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGT  >  NC_000913/1598489‑1598751                                                                                                                                 |                                                                                                                                        cGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAg                                                                                                                               >  2:244228/1‑139 (MQ=255)    aCTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                            <  1:758850/139‑1 (MQ=255)     cTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                           <  1:58601/139‑1 (MQ=255)      tGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATATTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                          <  2:413788/139‑1 (MQ=255)        gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                                        >  2:121316/1‑139 (MQ=255)                aCGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTTTCAATGAAATTATAACCCAGGCGCTTTCCATAacg                                                                                                                <  2:462884/139‑1 (MQ=255)                       aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                         <  2:774250/139‑1 (MQ=255)                           gCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                                     <  2:281369/139‑1 (MQ=255)                            cAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                                    <  1:527791/139‑1 (MQ=255)                                                  gCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTg                                                                             >  2:310634/1‑139 (MQ=255)                                                   ctctAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGa                                                                            <  1:121316/139‑1 (MQ=255)                                                   ctctAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGa                                                                            >  2:328884/1‑139 (MQ=255)                                                   ctctAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGa                                                                            <  1:328884/139‑1 (MQ=255)                                                              gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAAc                                                                 >  2:464139/1‑139 (MQ=255)                                                                           aaCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTa                                                    >  2:345841/1‑139 (MQ=255)                                                                           aaCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTa                                                    >  1:281634/1‑139 (MQ=255)                                                                                 aCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGt                                              <  1:244228/139‑1 (MQ=255)                                                                                     aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAActgcctgc                                                          <  1:523313/123‑1 (MQ=255)                                                                                     aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAActgcctgc                                                          >  2:523313/1‑123 (MQ=255)                                                                                     aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAAt                                          >  2:742096/1‑139 (MQ=255)                                                                                     aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAAt                                          <  2:768804/139‑1 (MQ=255)                                                                                           tAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGAt                                    <  2:781367/139‑1 (MQ=255)                                                                                             aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATTCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATa                                   >  1:492160/1‑139 (MQ=255)                                                                                                               atgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCggg                <  2:696172/139‑1 (MQ=255)                                                                                                                   ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGa                                                                                               <  1:603152/57‑1 (MQ=255)                                                                                                                   ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGa                                                                                               >  2:603152/1‑57 (MQ=255)                                                                                                                   ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg            >  2:325480/1‑139 (MQ=255)                                                                                                                   ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg            <  1:152250/139‑1 (MQ=255)                                                                                                                    gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                         >  1:7685/1‑46 (MQ=255)                                                                                                                    gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                         <  2:7685/46‑1 (MQ=255)                                                                                                                     cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGTCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTGCCGGGTCTGgt          <  1:464139/139‑1 (MQ=255)                                                                                                                     cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgt          <  1:409277/139‑1 (MQ=255)                                                                                                                     cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgt          >  2:589761/1‑139 (MQ=255)                                                                                                                        tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg       >  1:406097/1‑139 (MQ=255)                                                                                                                        tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg       >  1:244351/1‑139 (MQ=255)                                                                                                                          aaTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGtgtgtg     >  2:702973/1‑139 (MQ=255)                                                                                                                             gAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGt  <  1:325480/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                        CGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGA‑TCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTGCGT  >  NC_000913/1598489‑1598751

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