Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,814,218	C→T	*172* (TAG→TAA)	luxS ←	S‑ribosylhomocysteine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,814,218	0	C	T	100.0%	125.2 / NA	40	*172* (TAG→TAA)	luxS	S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

CTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTG  >  NC_000913/2814089‑2814347                                                                                                                                  |                                                                                                                                  ctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGc                                                                                                                          >  1:40104/1‑139 (MQ=255)                    cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGc                                                                                                                          <  1:909406/120‑1 (MQ=255)                    cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGc                                                                                                                          >  2:909406/1‑120 (MQ=255)                    cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCtt                                                                                                       <  2:1228644/139‑1 (MQ=255)                         ccGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                  <  2:524162/139‑1 (MQ=255)                            ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                                                      >  2:424637/1‑132 (MQ=255)                            ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                                                      <  1:424637/132‑1 (MQ=255)                            ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                                                                                               >  2:507731/1‑139 (MQ=255)                                 tGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAg                                                                                                                                <  1:880418/101‑1 (MQ=255)                                 tGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAg                                                                                                                                >  2:880418/1‑101 (MQ=255)                                           ggcggcGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                <  2:40104/139‑1 (MQ=255)                                                         tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                  <  2:222773/107‑1 (MQ=255)                                                         tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                  >  1:222773/1‑107 (MQ=255)                                                          gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                  <  2:681696/106‑1 (MQ=255)                                                          gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                  >  1:681696/1‑106 (MQ=255)                                                           cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGa                                                                         >  2:682146/1‑130 (MQ=255)                                                           cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGa                                                                         <  1:682146/130‑1 (MQ=255)                                                               cTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACg                                                            <  1:507731/139‑1 (MQ=255)                                                                   gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTcc                                                        >  1:335847/1‑139 (MQ=255)                                                                          tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATg                                                  <  1:770278/138‑1 (MQ=255)                                                                          tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATg                                                  >  2:770278/1‑138 (MQ=255)                                                                            aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGc                                                                               >  1:25474/1‑107 (MQ=255)                                                                            aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGc                                                                               <  2:25474/107‑1 (MQ=255)                                                                            aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                               >  1:1271601/1‑139 (MQ=255)                                                                                 gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACgcgc                                          >  1:379982/1‑139 (MQ=255)                                                                                           aCTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgc                                <  2:379982/139‑1 (MQ=255)                                                                                           aCTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTCGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgc                                >  1:344999/1‑139 (MQ=255)                                                                                            cTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcg                               >  1:1439907/1‑139 (MQ=255)                                                                                               gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgc                                >  2:378589/1‑135 (MQ=255)                                                                                               gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgc                                <  1:378589/135‑1 (MQ=255)                                                                                                          ttaattGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                 <  2:725714/139‑1 (MQ=255)                                                                                                                    aGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                <  2:198977/66‑1 (MQ=255)                                                                                                                    aGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                >  1:198977/1‑66 (MQ=255)                                                                                                                        aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcgc                              >  2:1122723/1‑112 (MQ=255)                                                                                                                        aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcgc                              <  1:1122723/112‑1 (MQ=255)                                                                                                                         gTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTg  >  1:167407/1‑139 (MQ=255)                                                                                                                          tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACg                                                                 <  2:805655/75‑1 (MQ=255)                                                                                                                          tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACg                                                                 >  1:805655/1‑75 (MQ=255)                                                                                                                           ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                               <  2:1127094/76‑1 (MQ=255)                                                                                                                           ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                               >  1:1127094/1‑76 (MQ=255)                                                                                                                                  |                                                                                                                                  CTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTG  >  NC_000913/2814089‑2814347

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