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

TATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGT  >  NC_000913/1598486‑1598743                                                                                                                                    |                                                                                                                               tatCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAAcc                                                                                                                         <  2:505584/139‑1 (MQ=255)      tACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCAAGgc                                                                                                                    >  2:157249/1‑139 (MQ=255)       aCTGGTGATGGAACGATCAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGCAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                   <  2:622374/139‑1 (MQ=255)           gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                   >  1:205464/1‑135 (MQ=255)           gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                   <  2:205464/135‑1 (MQ=255)           gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                               >  2:37965/1‑139 (MQ=255)                ggAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                          >  1:293248/1‑139 (MQ=255)                     gATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                     <  1:157249/139‑1 (MQ=255)                          aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  2:943897/139‑1 (MQ=255)                          aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  1:484578/139‑1 (MQ=255)                          aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAACAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  2:293248/139‑1 (MQ=255)                                       cTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                          >  1:196517/1‑116 (MQ=255)                                       cTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                          <  2:196517/116‑1 (MQ=255)                                                aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                          >  2:449756/1‑139 (MQ=255)                                                                 gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAc                                                         >  1:123690/1‑139 (MQ=255)                                                                     aaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                               <  2:468165/81‑1 (MQ=255)                                                                     aaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                               >  1:468165/1‑81 (MQ=255)                                                                         aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCTTGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacga                                                                                                      <  2:136890/86‑1 (MQ=255)                                                                         aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacga                                                                                                      >  1:136890/1‑86 (MQ=255)                                                                                    aCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  1:312930/81‑1 (MQ=255)                                                                                    aCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                >  2:312930/1‑81 (MQ=255)                                                                                     cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                     >  1:523808/1‑139 (MQ=255)                                                                                         aaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATc                                 <  2:523808/139‑1 (MQ=255)                                                                                                          tcatcaTGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                 >  1:883292/1‑42 (MQ=255)                                                                                                          tcatcaTGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                 <  2:883292/42‑1 (MQ=255)                                                                                                          tcatcaTGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                 <  2:787050/42‑1 (MQ=255)                                                                                                          tcatcaTGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                 >  1:787050/1‑42 (MQ=255)                                                                                                              catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                >  1:202369/1‑55 (MQ=255)                                                                                                              catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  2:202369/55‑1 (MQ=255)                                                                                                                  atgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                            <  2:320660/71‑1 (MQ=255)                                                                                                                  atgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                            >  1:320660/1‑71 (MQ=255)                                                                                                                   tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                >  1:484752/1‑50 (MQ=255)                                                                                                                   tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                <  2:484752/50‑1 (MQ=255)                                                                                                                   tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGt       <  2:123690/139‑1 (MQ=255)                                                                                                                   tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGt       <  1:37965/139‑1 (MQ=255)                                                                                                                    gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                   <  2:8606/110‑1 (MQ=255)                                                                                                                    gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                   >  1:8606/1‑110 (MQ=255)                                                                                                                    gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc            <  2:78712/133‑1 (MQ=255)                                                                                                                    gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc            >  1:78712/1‑133 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGc                                                                             >  2:844432/1‑66 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGc                                                                             <  1:844432/66‑1 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                       <  1:997743/104‑1 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                       >  2:997743/1‑104 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTc                     >  2:785675/1‑122 (MQ=255)                                                                                                                      ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTc                     <  1:785675/122‑1 (MQ=255)                                                                                                                       gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTgg   >  1:826127/1‑139 (MQ=255)                                                                                                                        cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATcc                                                                 >  2:762809/1‑76 (MQ=255)                                                                                                                        cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATcc                                                                 <  1:762809/76‑1 (MQ=255)                                                                                                                        cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                          >  2:512066/1‑83 (MQ=255)                                                                                                                        cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                          <  1:512066/83‑1 (MQ=255)                                                                                                                        cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgt  >  2:239495/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                               TATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGT  >  NC_000913/1598486‑1598743

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