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%	188.3 / NA	59	*531* (TAG→TAA)	lsrK	autoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/28);  total (31/28)

GTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTG  >  NC_000913/1598496‑1598741                                                                                                                          |                                                                                                                             gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                             >  2:1162984/1‑139 (MQ=255)    aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCa                                                                                                          >  2:1437894/1‑139 (MQ=255)           gATTAATCTGCAAGCTGGCTATTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGTCGCTTTCCATAacgac                                                                                                   >  1:1177795/1‑139 (MQ=255)            aTGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgacg                                                                                                  >  2:1611986/1‑139 (MQ=255)                    gCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                          <  2:1177795/139‑1 (MQ=255)                         cTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATg                                                                                     <  1:1437894/139‑1 (MQ=255)                           ggCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAATCCAGgcg                                                                                                                 >  2:654664/1‑109 (MQ=255)                           ggCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 <  1:654664/109‑1 (MQ=255)                                      aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                        >  1:678689/1‑139 (MQ=255)                                        ccAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAATCAGCCCCAg                                                                      <  1:1162984/139‑1 (MQ=255)                                          aGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCt                                                                    <  1:1611986/139‑1 (MQ=255)                                            ctctAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 >  2:1472273/1‑92 (MQ=255)                                            ctctAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 <  1:1472273/92‑1 (MQ=255)                                                 aaaGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              <  1:1435148/106‑1 (MQ=255)                                                 aaaGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              >  2:1435148/1‑106 (MQ=255)                                                         gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                <  2:775924/80‑1 (MQ=255)                                                         gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                >  1:775924/1‑80 (MQ=255)                                                         gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              >  1:1186153/1‑98 (MQ=255)                                                         gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              <  2:1186153/98‑1 (MQ=255)                                                         gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActg                                                     >  1:427884/1‑139 (MQ=255)                                                               aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                  <  2:873892/120‑1 (MQ=255)                                                               aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                  >  1:873892/1‑120 (MQ=255)                                                                     aCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTAt                                         <  2:678689/139‑1 (MQ=255)                                                                       tataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                     >  2:1588006/1‑109 (MQ=255)                                                                       tataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                     <  1:1588006/109‑1 (MQ=255)                                                                           cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                   <  2:1437614/75‑1 (MQ=255)                                                                           cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                   >  1:1437614/1‑75 (MQ=255)                                                                           cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                   >  1:1282133/1‑139 (MQ=255)                                                                           cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                   >  1:1432692/1‑139 (MQ=255)                                                                              aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAAt                                <  2:1432692/139‑1 (MQ=255)                                                                                aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAACTTATAACCCAGgcgc                                                                                                                >  1:1065598/1‑57 (MQ=255)                                                                                aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                <  2:1065598/57‑1 (MQ=255)                                                                                 aaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                >  2:1549352/1‑56 (MQ=255)                                                                                 aaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                <  1:1549352/56‑1 (MQ=255)                                                                                      aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 >  1:471306/1‑50 (MQ=255)                                                                                      aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 <  2:471306/50‑1 (MQ=255)                                                                                      aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              <  1:1008763/69‑1 (MQ=255)                                                                                      aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                              >  2:1008763/1‑69 (MQ=255)                                                                                              ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                 <  1:1001230/42‑1 (MQ=255)                                                                                              ggTCATCATGATGTGGCTGTCAATGAAATAATAACCCAGgcg                                                                                                                 >  2:1001230/1‑42 (MQ=255)                                                                                                         tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                    >  1:1511233/1‑92 (MQ=255)                                                                                                         tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                    <  2:1511233/92‑1 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGAt                                                                 >  1:1228251/1‑78 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                    >  2:220092/1‑91 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                    <  1:220092/91‑1 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                   >  1:385564/1‑108 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                   <  2:385564/108‑1 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                 >  1:1588780/1‑110 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                 <  2:1588780/110‑1 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc          >  2:727181/1‑133 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc          <  1:727181/133‑1 (MQ=255)                                                                                                          gtgGCTGTCAATGAAATTATAACCCAGGCGCTTCCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGAt                                                                 <  2:1228251/78‑1 (MQ=255)                                                                                                            ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGTTAAAGTTCATGCTTTTCCGGGTCTg  >  2:1068842/1‑139 (MQ=255)                                                                                                             gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                   <  2:587869/41‑1 (MQ=255)                                                                                                             gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                   >  1:587869/1‑41 (MQ=255)                                                                                                             gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                 <  2:570913/123‑1 (MQ=255)                                                                                                             gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                 >  1:570913/1‑123 (MQ=255)                                                                                                                 tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAc                                                                               >  1:1112696/1‑57 (MQ=255)                                                                                                                 tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAc                                                                               <  2:1112696/57‑1 (MQ=255)                                                                                                                          |                                                                                                                             GTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTG  >  NC_000913/1598496‑1598741

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