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%	134.1 / NA	42	*172* (TAG→TAA)	luxS	S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/21);  total (21/21)

GTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCAC  >  NC_000913/2814083‑2814351                                                                                                                                        |                                                                                                                                      gTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAg                                                                                                                                    >  2:1273015/1‑139 (MQ=255)      gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGAtgtg                                                                                                                               >  2:776593/1‑139 (MQ=255)          gggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGt                                                                                                                           <  1:296590/139‑1 (MQ=255)                               ccGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                      <  2:995814/139‑1 (MQ=255)                                  ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAAc                                                                                                                  >  2:776887/1‑124 (MQ=255)                                  ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAAc                                                                                                                  <  1:776887/124‑1 (MQ=255)                                          cgcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc                                                                                           <  2:23609/139‑1 (MQ=255)                                               gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGt                                                                                      >  1:1056743/1‑139 (MQ=255)                                                   cggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAg                                                                                                                                    <  2:88357/89‑1 (MQ=255)                                                   cggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAg                                                                                                                                    >  1:88357/1‑89 (MQ=255)                                                         aGGCGGTGCGCACTAGGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                                                            >  2:703966/1‑139 (MQ=255)                                                                  gcACTAGGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                                   <  1:703966/139‑1 (MQ=255)                                                                   cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                                                  <  2:251732/139‑1 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTTTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt                                                    <  2:798444/138‑1 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTTTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt                                                    >  1:798444/1‑138 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                                                                 <  2:1067745/93‑1 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                                                                 >  1:1067745/1‑93 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt                                                    <  2:1138633/138‑1 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt                                                    >  1:1138633/1‑138 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  1:972129/1‑139 (MQ=255)                                                                                  aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  2:1078795/1‑139 (MQ=255)                                                                                  aGCCAGTACATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  2:848899/1‑139 (MQ=255)                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                      <  2:1027273/86‑1 (MQ=255)                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                                                      >  1:1027273/1‑86 (MQ=255)                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                                                  <  2:998901/122‑1 (MQ=255)                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                                                  >  1:998901/1‑122 (MQ=255)                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcg                                                 >  2:384743/1‑139 (MQ=255)                                                                                       gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACgcgc                                              >  2:197423/1‑139 (MQ=255)                                                                                        ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCa                                             >  2:1261599/1‑139 (MQ=255)                                                                                                aaCTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTg                                     <  1:186141/139‑1 (MQ=255)                                                                                                      cTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                                                            <  1:11020/94‑1 (MQ=255)                                                                                                      cTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                                                            >  2:11020/1‑94 (MQ=255)                                                                                                               attaattGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGa                      >  1:1274524/1‑139 (MQ=255)                                                                                                               attaattGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGa                      >  1:384305/1‑139 (MQ=255)                                                                                                                          aGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGt           <  2:1274524/139‑1 (MQ=255)                                                                                                                               gTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTg      <  1:1261599/139‑1 (MQ=255)                                                                                                                                tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCa                                                           >  2:342562/1‑85 (MQ=255)                                                                                                                                tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCa                                                           <  1:342562/85‑1 (MQ=255)                                                                                                                                  tACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCca   <  1:124783/139‑1 (MQ=255)                                                                                                                                  tACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCca   <  1:1078795/139‑1 (MQ=255)                                                                                                                                   aCTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcg                                   >  2:244302/1‑106 (MQ=255)                                                                                                                                   aCTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcg                                   <  1:244302/106‑1 (MQ=255)                                                                                                                                   aCTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCcac  <  2:972129/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      GTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCAC  >  NC_000913/2814083‑2814351

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