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

TCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCAC  >  NC_000913/2814084‑2814351                                                                                                                                       |                                                                                                                                      tCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGa                                                                                                                                   >  1:138004/1‑139 (MQ=255)                tgtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCa                                                                                                                    >  1:639855/1‑139 (MQ=255)                         cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCtt                                                                                                           >  2:1184497/1‑139 (MQ=255)                             aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCgg                                                                                                       <  1:1676079/139‑1 (MQ=255)                             aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCgg                                                                                                       <  2:639855/139‑1 (MQ=255)                               cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCa                                                                                                     >  2:583379/1‑139 (MQ=255)                                          gcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttct                                                                                          <  2:1782512/139‑1 (MQ=255)                                              gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGt                                                                                      >  2:805311/1‑139 (MQ=255)                                                          gCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGc                                                                          >  1:1436463/1‑139 (MQ=255)                                                                cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACg                                                                     >  1:605612/1‑138 (MQ=255)                                                                cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACg                                                                     <  2:605612/138‑1 (MQ=255)                                                                  cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                                                  <  1:604554/139‑1 (MQ=255)                                                                      aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGtt                                                              <  1:583379/139‑1 (MQ=255)                                                                              ctggcttaGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTa                                                                       >  2:509709/8‑122 (MQ=255)                                                                              ctggcttaGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTa                                                                       <  1:509709/115‑1 (MQ=255)                                                                               tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                    <  1:1755925/108‑1 (MQ=255)                                                                               tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                    >  2:1755925/1‑108 (MQ=255)                                                                                 aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATg                                                                           >  1:1854999/1‑115 (MQ=255)                                                                                 aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATg                                                                           <  2:1854999/115‑1 (MQ=255)                                                                                 aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  2:863942/1‑139 (MQ=255)                                                                                 aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  2:1583171/1‑139 (MQ=255)                                                                                 aGCCAGATCATTTGAACTGGCTTTTTTCAATTATTTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                                                   >  1:1740684/1‑139 (MQ=255)                                                                                  gCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                  <  2:1436463/139‑1 (MQ=255)                                                                                  gCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                  <  1:1583171/139‑1 (MQ=255)                                                                                   ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcg                                                 >  2:1683337/1‑139 (MQ=255)                                                                                      gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACgcgc                                              >  1:72875/1‑139 (MQ=255)                                                                                       ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                                                                                               <  1:157732/89‑1 (MQ=255)                                                                                       ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                                                                                               >  2:157732/1‑89 (MQ=255)                                                                                       ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCa                                             <  1:805311/139‑1 (MQ=255)                                                                                                 cTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcg                                   >  2:1765790/1‑139 (MQ=255)                                                                                                          ttCAATTAATTGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCa                          <  1:1859209/139‑1 (MQ=255)                                                                                                               ttaattGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     <  2:72875/139‑1 (MQ=255)                                                                                                                  attGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGa                                                         <  2:579553/100‑1 (MQ=255)                                                                                                                  attGTGAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGa                                                         >  1:579553/1‑100 (MQ=255)                                                                                                                    tgtgAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     >  1:1637359/1‑134 (MQ=255)                                                                                                                    tgtgAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     <  2:1637359/134‑1 (MQ=255)                                                                                                                     gtgAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATc               >  2:193217/1‑139 (MQ=255)                                                                                                                       gAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     <  1:110786/131‑1 (MQ=255)                                                                                                                       gAAGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     >  2:110786/1‑131 (MQ=255)                                                                                                                        aaGATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTgg            >  1:531110/1‑139 (MQ=255)                                                                                                                          gATAGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTa          <  2:138004/139‑1 (MQ=255)                                                                                                                             aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                                                            >  2:1194880/1‑70 (MQ=255)                                                                                                                             aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                                                            <  1:1194880/70‑1 (MQ=255)                                                                                                                             aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                                   <  2:1848243/79‑1 (MQ=255)                                                                                                                             aGTTTACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                                   >  1:1848243/1‑79 (MQ=255)                                                                                                                               tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                                                                 >  2:1013302/1‑47 (MQ=255)                                                                                                                               tttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                                                                 <  1:1013302/47‑1 (MQ=255)                                                                                                                                ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                  >  2:452732/1‑93 (MQ=255)                                                                                                                                ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                  <  1:452732/93‑1 (MQ=255)                                                                                                                                ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                     <  1:1143629/122‑1 (MQ=255)                                                                                                                                ttACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGcc    >  1:799083/1‑139 (MQ=255)                                                                                                                                ttACTGATTAGATGAGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTAAGTCACGTTCCAGAATGCTACGCGCAACATCCTGCGCTTCCTGCAACGAg                     >  2:1143629/1‑122 (MQ=255)                                                                                                                                 tACTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCca   <  1:863942/139‑1 (MQ=255)                                                                                                                                  aCTGATTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTGGCTGTTGATGCGTACGTCAGGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCcac  <  2:1740684/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCAC  >  NC_000913/2814084‑2814351

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