Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	531,988	A→G	S78S (AGT→AGC)	allS ←	DNA‑binding transcriptional activator AllS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	531,988	0	A	G	100.0%	88.2 / NA	26	S78S (AGT→AGC)	allS	DNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/13);  total (13/13)

GGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATACTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACTCCGGTATTCTCTTCCAGAAGTT  >  NC_000913/531863‑532109                                                                                                                              |                                                                                                                          gggTAACGCTCTTCAAGCCATGCCAGCAACTGGCCGATGGCCTGGGGGTTGTCGAGCAGGTTGTTGATGACATTATTCACCTGGCGTTCCACGCTGTCTTTCATCTGTTGTAGCTCGCTTGGCAGGCTTTcta                                                                                                                    >  2:510248/1‑131 (MQ=255)    tAACGCTCATTCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTccagcc                                                                                                                 >  2:417317/1‑133 (MQ=255)      aCGCTCATTCTGCCACGCCAGCAACTGGTCGAAGGCCTGGGTGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTccagccag                                                                                                               <  2:200536/133‑1 (MQ=255)            aTTCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGATGTAGAGCAGGCTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGc                                                                                                         >  2:297548/1‑133 (MQ=255)             ttCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGcc                                                                                                        <  1:267828/133‑1 (MQ=255)             ttCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGcc                                                                                                        <  2:365095/133‑1 (MQ=255)                      gCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTgg                                                                                               <  1:417317/133‑1 (MQ=255)                      gCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTgg                                                                                               <  1:236570/133‑1 (MQ=255)                                gggCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAg                                                                                     >  1:474320/1‑133 (MQ=255)                                 ggCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGt                                                                                    <  2:474320/133‑1 (MQ=255)                                      cGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATg                                                                               >  1:494895/1‑133 (MQ=255)                                                   gTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCt                                                                  >  1:459658/1‑133 (MQ=255)                                                           ggTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGt                                                          >  2:304874/1‑133 (MQ=255)                                                                         atatTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTct                                            <  2:494895/133‑2 (MQ=255)                                                                               aCCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCCCTGTCAAAGTCGCGATGCGAGGCGTACGGa                                      >  1:228449/1‑133 (MQ=255)                                                                                   ggCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGaaaaa                                  >  1:437893/1‑133 (MQ=255)                                                                                    gCGTTCCACGCCATCATTGACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAAc                                 >  1:33306/1‑133 (MQ=255)                                                                                        tCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAgcg                             <  1:304874/133‑1 (MQ=255)                                                                                         ccACGCCATCATTGACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAgcgc                            <  2:33306/133‑1 (MQ=255)                                                                                         ccACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAgcgc                            >  2:361949/1‑133 (MQ=255)                                                                                         ccACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGAAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAccgc                            >  2:291454/1‑133 (MQ=255)                                                                                           aCGCCATCATTCACCAGTTGCAGCTCACTTGGTATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTa                          <  2:228449/133‑1 (MQ=255)                                                                                            cGCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTAc                         <  2:25611/133‑1 (MQ=255)                                                                                             gCCATCATTCACCTGTTGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACt                        <  2:459658/133‑1 (MQ=255)                                                                                                   atTCACCTGTTGCAGCTCTCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACTCCGGTa                  <  1:198286/133‑1 (MQ=255)                                                                                                            ttGCAGCTCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACTCCGGTATTCTCTTcc         >  1:395734/1‑133 (MQ=255)                                                                                                                   tCGCTTGGCATGCTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACTCCGGTATTCTCTTCCAGAAGtt  >  1:93889/1‑133 (MQ=255)                                                                                                                              |                                                                                                                          GGGTAACGCTCATTCAGCCACGCCAGCAACTGGGCGACGGCCTGGGGGTTGTAGAGCAGGTTGTTGATGACAATATTCACCTGGCGTTCCACGCCATCATTCACCTGTTGCAGCTCGCTTGGCATACTTTCCAGCCAGCTCAGCCAGTCTCTGGCCTGGGAAAGTAGATGCTCGCCAGCCGCTGTCAACGTCACGCTGCGAGTCGTACGGAAAAACAGCGCTACTCCGGTATTCTCTTCCAGAAGTT  >  NC_000913/531863‑532109

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