Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,192	G→A	Q447Q (CAG→CAA)	ydgA →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,192	0	G	A	100.0%	88.0 / NA	28	Q447Q (CAG→CAA)	ydgA	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/15);  total (13/15)

ACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGG  >  NC_000913/1691070‑1691273                                                                                                                           |                                                                                  aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   >  1:198427/1‑139 (MQ=255) aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   >  2:356832/1‑139 (MQ=255) aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   <  1:827896/139‑1 (MQ=255) aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   <  1:606275/139‑1 (MQ=255)                  ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGAGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATAcc                                                  >  1:442334/1‑139 (MQ=255)                  ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGAGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATAcc                                                  <  2:442334/139‑1 (MQ=255)                     aaCTGAGTTTATGACTCAGGTAGCGAAGGTGGAAGGTTATCCGGAAGATCACGCGACAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCTGTCTGACCACCTTGCAGGACAATACCATc                                               >  2:556411/1‑139 (MQ=255)                          aGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                          <  2:198427/139‑1 (MQ=255)                            tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGTGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                <  2:515345/131‑1 (MQ=255)                            tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                >  1:515345/1‑131 (MQ=255)                             ttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAg                                       <  2:460284/139‑1 (MQ=255)                                    tCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                              >  2:712501/1‑109 (MQ=255)                                    tCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                              <  1:712501/109‑1 (MQ=255)                                          aGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCt                          <  1:575412/139‑1 (MQ=255)                                             gAAGCTGGAAGGTTATTAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                          <  1:506271/120‑1 (MQ=255)                                             gAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAGTGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                          >  2:506271/1‑120 (MQ=255)                                                gCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGt                    <  1:356832/139‑1 (MQ=255)                                                                gAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                <  2:406418/95‑1 (MQ=255)                                                                gAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                >  1:406418/1‑95 (MQ=255)                                                                  aGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATggg  >  2:28351/1‑139 (MQ=255)                                                                        aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATg                                                                                <  1:575238/55‑1 (MQ=255)                                                                        aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATg                                                                                >  2:575238/1‑55 (MQ=255)                                                                                                   tGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   >  2:32520/1‑41 (MQ=255)                                                                                                   tGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                   <  1:32520/41‑1 (MQ=255)                                                                                                      aGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTg                                   >  1:485186/1‑70 (MQ=255)                                                                                                      aGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTg                                   <  2:485186/70‑1 (MQ=255)                                                                                                            aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                          >  2:296940/1‑57 (MQ=255)                                                                                                            aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                          <  1:296940/57‑1 (MQ=255)                                                                                                                           |                                                                                  ACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGG  >  NC_000913/1691070‑1691273

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