Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	994,959	A→G	L344L (TTG→CTG)	ssuD ←	FMNH2‑dependent alkanesulfonate monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	994,959	0	A	G	100.0%	100.2 / NA	29	L344L (TTG→CTG)	ssuD	FMNH2‑dependent alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/14);  total (15/14)

CCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGGC  >  NC_000913/994881‑995031                                                                               |                                                                         ccACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCt                                                                        <  2:185240/81‑1 (MQ=255)  cACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTc                                                                       <  1:185241/81‑1 (MQ=255)      gCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGcc                                                                    >  2:185277/1‑80 (MQ=255)      gCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGcc                                                                    >  2:185273/1‑80 (MQ=255)       cTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCaa                                                                  >  2:185278/1‑81 (MQ=255)          cGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAAccc                                                               >  2:185276/1‑81 (MQ=255)          cGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAAccc                                                               >  1:185285/1‑81 (MQ=255)           gCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAAccc                                                               >  2:185275/1‑80 (MQ=255)           gCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAAccc                                                               <  2:185242/80‑1 (MQ=255)             cTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGa                                                             >  1:185281/1‑80 (MQ=255)               tGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATAc                                                          >  1:185268/1‑81 (MQ=255)                  ggATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGct                                                       >  2:185283/1‑81 (MQ=255)                    aTTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGcttc                                                     <  2:185243/81‑1 (MQ=255)                       cAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGcttctt                                                   <  2:185244/80‑1 (MQ=255)                                cggctgGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGc                                           <  1:185245/79‑1 (MQ=255)                                     ggggAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGc                                    >  1:185262/1‑81 (MQ=255)                                      gggAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAg                                     <  2:185246/79‑1 (MQ=255)                                         aaTTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGa                                <  2:185247/81‑1 (MQ=255)                                         aaTTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGa                                <  2:185248/81‑1 (MQ=255)                                              ccGGGATGGCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGca                           >  1:185288/1‑81 (MQ=255)                                               cGGGATGGCGACATCCAGAAGCGGGAACAGCGGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGcac                          <  2:185249/81‑1 (MQ=255)                                                      gCGACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACt                    <  1:185250/80‑1 (MQ=255)                                                        gACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTc                 <  1:185251/81‑1 (MQ=255)                                                        gACATCCAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTc                 >  2:185299/1‑81 (MQ=255)                                                             ccAGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGc            >  2:185282/1‑81 (MQ=255)                                                               aGAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGcc           >  1:185293/1‑80 (MQ=255)                                                                gAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGcc           >  1:185287/1‑79 (MQ=255)                                                                gAAGCGGGAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCaa         <  1:185252/81‑1 (MQ=255)                                                                       gAACAGCAGCTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGGc  <  1:185253/81‑1 (MQ=255)                                                                               |                                                                         CCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGGC  >  NC_000913/994881‑995031

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