Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,603,562	T→C	F182L (TTT→CTT)	lsrC →	autoinducer 2 import system permease protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,603,562	0	T	C	100.0%	117.4 / NA	32	F182L (TTT→CTT)	lsrC	autoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/15);  total (17/15)

GTGGACTGGCGGCAAATGGATTGAAGGGTTACCCGCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATT  >  NC_000913/1603420‑1603705                                                                                                                                               |                                                                                                                                                gtgGACTGGCGGCAAATGGATTGAAGGGTTACCCGCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGAc                                                                                                                                           >  1:645325/1‑149 (MQ=255)          cggcAAATGGATTGAAGGGTTACCCGCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTAt                                                                                                                                >  1:593805/1‑151 (MQ=255)            gcAAATGGATTGAAGGGTTACCCGCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGc                                                                                                                              >  2:515063/1‑151 (MQ=255)                                  cgccgAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGgc                                                                                                        >  2:412141/1‑151 (MQ=255)                                   gccgAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGgcg                                                                                                       >  2:139850/1‑151 (MQ=255)                                     cgAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCt                                                                                                     <  2:77273/151‑1 (MQ=255)                                      gAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTc                                                                                                    <  2:553855/151‑1 (MQ=255)                                      gAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTc                                                                                                    >  2:446262/1‑151 (MQ=255)                                           gAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCa                                                                                                >  1:191443/1‑150 (MQ=255)                                            aaaCAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAAc                                                                                              >  2:193858/1‑151 (MQ=255)                                                 gCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGCCGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGCTAATTTACAGGGCGCTCGTCAACTGGGc                                                                                         >  2:668387/1‑151 (MQ=255)                                                      ccGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGtt                                                                                      <  2:589750/149‑1 (MQ=255)                                                       cGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGt                                                                                   >  2:489779/1‑151 (MQ=255)                                                        gCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTa                                                                                  <  1:560227/151‑1 (MQ=255)                                                             gctgctgctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGa                                                                              >  1:617402/1‑150 (MQ=255)                                                              ctgctgctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGa                                                                              >  2:657723/1‑149 (MQ=255)                                                               tgctgctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGc                                                                           >  2:384544/1‑151 (MQ=255)                                                                  tgctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCAt                                                                        >  1:341218/1‑151 (MQ=255)                                                                   gctTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCAtt                                                                       <  2:316603/151‑1 (MQ=255)                                                                      tGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTc                                                                      <  1:499914/149‑1 (MQ=255)                                                                       ggCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCa                                                                   <  1:5701/151‑1 (MQ=255)                                                                            ttCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATtgt                                                               <  2:521328/150‑1 (MQ=255)                                                                              cAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTgg                                                             <  1:405867/150‑1 (MQ=255)                                                                                   aTTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCAtttt                                                       <  2:315080/151‑1 (MQ=255)                                                                                             tgACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACg                                              <  1:207822/150‑1 (MQ=255)                                                                                               aCGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCt                                           <  1:498469/151‑1 (MQ=255)                                                                                                cGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTg                                          <  1:653845/151‑1 (MQ=255)                                                                                                       tCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATggc                                    <  1:640943/150‑1 (MQ=255)                                                                                                                           ggCTGCTGGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATtgtg                 <  2:120049/149‑1 (MQ=255)                                                                                                                                 tgGCAAAGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCtt          >  1:397084/1‑150 (MQ=255)                                                                                                                                       aGACGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTTCTTCGCCGAt   >  1:481384/1‑151 (MQ=255)                                                                                                                                          cGGCGCTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTTTGTTTTCTTTGCAGGtt  >  2:207950/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                GTGGACTGGCGGCAAATGGATTGAAGGGTTACCCGCCGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATT  >  NC_000913/1603420‑1603705

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