Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,086,426	Δ1 bp	coding (225/252 nt)	yqgD ←	DUF2684 domain‑containing protein YqgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,086,423	0	C	.	100.0%	102.8 / NA	25	coding (228/252 nt)	yqgD	DUF2684 domain‑containing protein YqgD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/13);  total (12/13)

CAACCCGGCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATACCCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCAGCCACG  >  NC_000913/3086350‑3086496                                                                          |                                                                          cAACCCGGCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAAt                                                                   >  2:509393/1‑81 (MQ=255)  aaCCCGGCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAAt                                                                   <  2:509352/80‑1 (MQ=255)       ggCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAg                                                              >  1:509390/1‑80 (MQ=255)       ggCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAg                                                              >  1:509391/1‑80 (MQ=255)        gCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGa                                                             >  2:509384/1‑80 (MQ=255)         cTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGaaa                                                           <  2:509353/81‑1 (MQ=255)         cTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGaa                                                            >  2:509382/1‑80 (MQ=255)              aGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGc                                                      <  2:509354/81‑1 (MQ=255)              aGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGc                                                      <  2:509355/81‑1 (MQ=255)                   cAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTGTAATA‑CCCACAATGGCAGAAAATTGCAAAAg                                                 <  2:509356/81‑1 (MQ=255)                       acGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATa                                              >  1:509397/1‑80 (MQ=255)                        cGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATaaa                                            <  1:509358/81‑1 (MQ=255)                          aaagaaaCGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATa                                          >  2:509401/1‑81 (MQ=255)                             gaaaCGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGc                                       <  2:509358/81‑1 (MQ=255)                               aaCGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAg                                     <  2:509359/81‑1 (MQ=255)                               aaCGTCGTGCGCTTTTTATTTAAGCCGCGCGCCGCGTTTTACA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACg                                        >  1:509401/1‑78 (MQ=255)                                aCGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGa                                    <  2:509360/81‑1 (MQ=255)                                     gtgtGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGcc                               <  2:509361/81‑1 (MQ=255)                                            ttttATTTAAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTg                        >  2:509392/1‑81 (MQ=255)                                                atttaAGCCGCGCGCCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCa                     <  2:509362/80‑1 (MQ=255)                                                         gcgcgcCGCGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAAttt           >  1:509392/1‑81 (MQ=255)                                                            cgccgcGTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCa        <  2:509363/81‑1 (MQ=255)                                                                 cgTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCAGCCAc   <  2:509364/81‑1 (MQ=255)                                                                  gTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCAGCCACg  >  1:509395/1‑81 (MQ=255)                                                                  gTTTTATA‑CCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCAGCCACg  >  2:509408/1‑81 (MQ=255)                                                                          |                                                                          CAACCCGGCTGGAAGTGGCAACACGAAAGAAACGTCGTGTGCTTTTTATTTAAGCCGCGCGCCGCGTTTTATACCCCACAATGGCAGAAAATTGCAAAAGATAAATACGCAGAATGCCGGCATTGTCAGGAAAAATTTCCAGCCACG  >  NC_000913/3086350‑3086496

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