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

TACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTG  >  NC_000913/1598491‑1598748                                                                                                                               |                                                                                                                                    tACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgc                                                                                                                         >  1:289459/1‑139 (MQ=255)  aCTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                        <  1:104482/139‑1 (MQ=255)   cTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcgc                                                                                                                       <  2:767901/139‑1 (MQ=255)         aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                      >  1:1061927/1‑134 (MQ=255)         aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                      <  2:1061927/134‑1 (MQ=255)                     aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                     <  2:878116/139‑1 (MQ=255)                     aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                     <  1:1133427/139‑1 (MQ=255)                     aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                     <  2:289459/139‑1 (MQ=255)                         gCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                                 <  1:102165/139‑1 (MQ=255)                                           aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                               >  2:641846/1‑139 (MQ=255)                                              cAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                            <  1:641846/139‑1 (MQ=255)                                                           agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTAGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGCTCCTGATaaa                                                               >  1:298104/1‑139 (MQ=255)                                                            gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAc                                                              >  2:544903/1‑139 (MQ=255)                                                            gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAc                                                              >  1:654295/1‑139 (MQ=255)                                                               aaaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                           <  1:1219617/139‑1 (MQ=255)                                                                    aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                     <  2:678754/92‑1 (MQ=255)                                                                    aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                     >  1:678754/1‑92 (MQ=255)                                                                        tAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAcc                                                                                     <  1:508226/104‑1 (MQ=255)                                                                        tAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAcc                                                                                     >  2:508226/1‑104 (MQ=255)                                                                             ataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcg                                                                                                      <  2:777493/82‑1 (MQ=255)                                                                             ataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcg                                                                                                      >  1:777493/1‑82 (MQ=255)                                                                                cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                       >  1:13176/1‑110 (MQ=255)                                                                                cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                       <  2:13176/110‑1 (MQ=255)                                                                                cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                          >  1:461307/1‑139 (MQ=255)                                                                                   aaaaaCTAACAGCCACGGTCATCATGATGTTGCTGTCAATGAAATtat                                                                                                                                  >  1:89439/1‑48 (MQ=255)                                                                                   aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATtat                                                                                                                                  <  2:89439/48‑1 (MQ=255)                                                                                   aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                        >  1:586965/1‑58 (MQ=255)                                                                                   aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                        <  2:586965/58‑1 (MQ=255)                                                                                     aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATcc                                     >  1:893850/1‑139 (MQ=255)                                                                                           aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                        <  1:65843/50‑1 (MQ=255)                                                                                           aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                        >  2:65843/1‑50 (MQ=255)                                                                                                  cGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCAt                        >  1:608183/1‑139 (MQ=255)                                                                                                       atcatGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAct                                                             >  2:139724/1‑97 (MQ=255)                                                                                                       atcatGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAct                                                             <  1:139724/97‑1 (MQ=255)                                                                                                       atcatGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCttt                   >  2:1126313/1‑139 (MQ=255)                                                                                                         catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 <  2:893850/139‑1 (MQ=255)                                                                                                         catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 <  2:608183/139‑1 (MQ=255)                                                                                                             atgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 >  1:673295/1‑135 (MQ=255)                                                                                                             atgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 <  2:673295/135‑1 (MQ=255)                                                                                                               gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                       <  1:295861/63‑1 (MQ=255)                                                                                                               gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                       >  2:295861/1‑63 (MQ=255)                                                                                                                 ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg         >  2:300445/1‑139 (MQ=255)                                                                                                                   cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgt       >  1:1179571/1‑139 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACgtcgtc                                                                                                    >  1:1083966/1‑43 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACgtcgtc                                                                                                    <  2:1083966/43‑1 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTcc                                                                                               <  1:500967/48‑1 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTcc                                                                                               >  2:500967/1‑48 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGGTCAACCAGCCCCAGCTGATCCTGATAAActgc                                                           <  1:524496/84‑1 (MQ=255)                                                                                                                      tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                           >  2:524496/1‑84 (MQ=255)                                                                                                                        aaTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGtgtgtg  <  2:654295/139‑1 (MQ=255)                                                                                                                           gAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 <  1:339252/121‑1 (MQ=255)                                                                                                                           gAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                 >  2:339252/1‑121 (MQ=255)                                                                                                                               |                                                                                                                                    TACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGTG  >  NC_000913/1598491‑1598748

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