Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,599,375	G→A	Q279* (CAA→TAA)	lsrK ←	autoinducer‑2 kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,599,375	0	G	A	100.0%	113.2 / NA	33	Q279* (CAA→TAA)	lsrK	autoinducer‑2 kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/14);  total (19/14)

GCATCGCGGAACCAGCGCATGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTGCTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGCCTTCAGACCGCAGAGTTCCGCCGCTT  >  NC_000913/1599239‑1599509                                                                                                                                          |                                                                                                                                       gCATCGCGGAACCAGCGCATGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTtactg                                                                                                                                     <  1:613583/140‑1 (MQ=255)      gcgGAACCAGCGCATGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGa                                                                                                                                >  1:437847/1‑140 (MQ=255)                cgcATGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTgccgcc                                                                                                                      >  1:431462/1‑140 (MQ=255)                    tGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCTAAAATTTACAACTTACTGCCAGAATGTGCCGCCAAg                                                                                                                   >  2:359880/1‑140 (MQ=255)                        gAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACcg                                                                                                              <  2:437847/140‑1 (MQ=255)                             cGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCTAAAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGt                                                                                                          <  1:359880/140‑1 (MQ=255)                               gTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTtg                                                                                                       >  1:621528/1‑140 (MQ=255)                                                      ttGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAAcc                                                                                >  2:817727/1‑140 (MQ=255)                                                         tACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCaa                                                                             >  1:412029/1‑140 (MQ=255)                                                           ccATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATATCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAgg                                                                           <  1:741280/140‑1 (MQ=255)                                                              tGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCaa                                                                        >  2:803125/1‑140 (MQ=255)                                                              tGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCaa                                                                        <  2:621528/140‑1 (MQ=255)                                                              tGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCaa                                                                        >  2:562781/1‑140 (MQ=255)                                                              tGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCaa                                                                        >  2:186086/1‑140 (MQ=255)                                                               gCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGc                                                                                      <  2:840453/125‑1 (MQ=255)                                                               gCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGc                                                                                      >  1:840453/1‑125 (MQ=255)                                                                ccAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAAcc                                                                      <  1:858963/140‑1 (MQ=255)                                                                ccAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAAcc                                                                      <  1:186086/140‑1 (MQ=255)                                                                ccAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAAcc                                                                      <  1:803125/140‑1 (MQ=255)                                                                          cATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACg                                                            >  2:552775/1‑140 (MQ=255)                                                                             gAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCg                                                         <  1:817727/140‑1 (MQ=255)                                                                               ggATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGcc                                                       >  2:254693/1‑140 (MQ=255)                                                                               ggATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGcc                                                       >  2:257794/1‑140 (MQ=255)                                                                               ggATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGcc                                                       <  1:48211/140‑1 (MQ=255)                                                                                gATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGcct                                                      >  1:618696/1‑140 (MQ=255)                                                                                     aCGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCtcctcc                                                 >  2:579702/1‑140 (MQ=255)                                                                                      cgcgCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCa                                                <  1:552775/140‑1 (MQ=255)                                                                                                       gggTCTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGcc                               <  2:618696/140‑1 (MQ=255)                                                                                                           cTGTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGcc                           >  1:444491/1‑140 (MQ=255)                                                                                                             gTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGcc                               >  1:61479/1‑134 (MQ=255)                                                                                                             gTCACCGGCGCGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGcc                               <  2:61479/134‑1 (MQ=255)                                                                                                                    gcgcGGCT‑AAATTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGCCTTCAGACCg                  >  1:659883/1‑140 (MQ=255)                                                                                                                             aaaTTTACAACTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGCCTTCAGACCGCAGAGTTc          <  2:444491/140‑1 (MQ=255)                                                                                                                                     aaCTTACTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGCCTTCAGACCGCAGAGTTCCGCCGCtt  >  2:270556/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       GCATCGCGGAACCAGCGCATGGTGAGTCCGGTAAAAAAGCTTATAGATTCAGCTTGTACCATGCCAGGAATAACATGAGGATTAACGCGCACGTTCATTTCTGGGTCTGTCACCGGCGCGGCT‑AAATTTACAACTTGCTGCCAGAATGTGCCGCCAAGAACCGCGGTTTGTGCCGGACGCACAACGCCTAACCCAAGGCAACCAAGCTGCACGTCGCCTCCTCCAACGACCACCGGAGTGCCCGCCTTCAGACCGCAGAGTTCCGCCGCTT  >  NC_000913/1599239‑1599509

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