Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,678,711	(A)5→4	coding (285/945 nt)	ydgH →	DUF1471 domain‑containing protein YdgH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,678,707	0	A	.	100.0%	74.6 / NA	18	coding (281/945 nt)	ydgH	DUF1471 domain‑containing protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/10);  total (8/10)

TTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGA  >  NC_000913/1678629‑1678778                                                                               |                                                                        tttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaa                                                                        <  2:303926/80‑1 (MQ=255)      gttgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaag                                                                     >  1:303971/1‑77 (MQ=255)           cGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCagaaga                                                             >  2:303967/1‑80 (MQ=255)           cGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCagaaga                                                             >  1:303961/1‑80 (MQ=255)            gACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAgaagaa                                                            <  2:303928/80‑1 (MQ=255)                   cTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTaa                                                    >  1:303963/1‑81 (MQ=255)                      aTTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAAt                                                   <  1:303929/79‑1 (MQ=255)                       ttttGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATcgc                                                >  1:303968/1‑81 (MQ=255)                           ggTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTa                                             <  1:303930/80‑1 (MQ=255)                            gTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTaa                                            <  1:303931/80‑1 (MQ=255)                                cAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAAc                                       >  1:303980/1‑81 (MQ=255)                                     gTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtgt                                  <  1:303932/81‑1 (MQ=255)                                      tAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGtgtt                                 >  1:303959/1‑81 (MQ=255)                                             cGTGTGGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAAc                          <  1:303933/81‑1 (MQ=255)                                                 tgGTCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGcc                      <  1:303934/81‑1 (MQ=255)                                                    tCGCTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGaa                   <  1:303935/81‑1 (MQ=255)                                                       cTGACCTCTATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGa                    >  1:303966/1‑77 (MQ=255)                                                            ctctATAAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAgg           >  1:303970/1‑81 (MQ=255)                                                                 taAAGCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCt         <  2:303936/78‑1 (MQ=255)                                                                     gCCGATGCTG‑AAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTgt  <  2:303937/81‑2 (MQ=255)                                                                               |                                                                        TTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGA  >  NC_000913/1678629‑1678778

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