Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,821,424	G→C	G196G (GGC→GGG)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,821,424	0	G	C	100.0%	116.7 / NA	32	G196G (GGC→GGG)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (19/13);  total (19/13)

ACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAC  >  NC_000913/2821292‑2821554                                                                                                                                     |                                                                                                                                   aCCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCcccccc                                                                                                                              >  1:370097/1‑139 (MQ=255)     aTACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCaaa                                                                                                                          >  2:999430/1‑139 (MQ=255)      tACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAt                                                                                                                         >  2:75932/1‑139 (MQ=255)      tACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAt                                                                                                                         >  1:261572/1‑139 (MQ=255)              tCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGc                                                                                                                 <  2:261572/139‑1 (MQ=255)                 aCAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGt                                                                                                              <  2:684013/139‑1 (MQ=255)                     aCGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATc                                                                                                          >  1:678938/1‑139 (MQ=255)                        gCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTa                                                                                                       >  2:541983/1‑139 (MQ=255)                             ggCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTagaaga                                                                                                  <  2:855052/139‑1 (MQ=255)                                   ggTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCgg                                                                                            >  2:478467/1‑139 (MQ=255)                                       ccATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCa                                                                                        <  1:206061/139‑1 (MQ=255)                                                cATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                               <  2:1051293/139‑1 (MQ=255)                                                     gCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGAc                                                                          <  2:1238506/139‑1 (MQ=255)                                                                         gAAGACGATGTTCCAGATCTCAATTTAGCGGTCGCCGGCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAg                                                      <  1:545629/139‑1 (MQ=255)                                                                                           ctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTAt                                    >  1:368946/1‑139 (MQ=255)                                                                                              aaTGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgc                                 >  1:1157646/1‑139 (MQ=255)                                                                                              aaTGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgc                                 >  2:679007/1‑139 (MQ=255)                                                                                                tGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgcgc                               >  1:1056286/1‑139 (MQ=255)                                                                                                 gTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCa                              <  1:1095379/139‑1 (MQ=255)                                                                                                     cGGTCGCCGTCTTCTTCCGGGCTTCCCGGGGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGc                                                                                >  1:99462/1‑85 (MQ=255)                                                                                                     cGGTCGCCGTCTTCTTCCGGGCTTCCCGGGGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGc                                                                                <  2:99462/85‑1 (MQ=255)                                                                                                     cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCtt                          >  2:418081/1‑139 (MQ=255)                                                                                                     cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCtt                          >  2:423742/1‑139 (MQ=255)                                                                                                     cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCtt                          >  1:1149043/1‑139 (MQ=255)                                                                                                                  cttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGAt             <  2:1056286/139‑1 (MQ=255)                                                                                                                  cttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGAt             <  1:478467/139‑1 (MQ=255)                                                                                                                     ccGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCg          <  1:423742/139‑1 (MQ=255)                                                                                                                          cTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaataa     >  2:1141566/1‑139 (MQ=255)                                                                                                                             cccGGAGGCCCCCCCCAAATGTGGTTGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAc  >  2:342390/1‑139 (MQ=255)                                                                                                                             cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                               <  2:509778/62‑1 (MQ=255)                                                                                                                             cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                               >  1:509778/1‑62 (MQ=255)                                                                                                                             cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAc  >  2:277347/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   ACCCATACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATAC  >  NC_000913/2821292‑2821554

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