Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,948	G→C	R25R (CGG→CGC)	yicG →	PF03458 family inner membrane protein YicG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,818,948	0	G	C	100.0%	91.6 / NA	30	R25R (CGG→CGC)	yicG	PF03458 family inner membrane protein YicG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (18/12);  total (18/12)

TTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTTT  >  NC_000913/3818704‑3819173                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  ttttCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATg                                                                                                                                                                                                                             >  2:232032/1‑251 (MQ=255)    tCTTATGCATTTTCTCAGAAAGGAAGACCGAAAAATGCATCTACT‑CTACCGAACATAGCGATTTCCCCTAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACACATCGCTACCCTGCCAGACAGATTTTTAGGTAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGTCGGCCGGTCGTCGCCGCATGGAt                                                                                                                                                                                                                          >  2:424412/1‑251 (MQ=255)      ttATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATAc                                                                                                                                                                                                                        >  2:219269/1‑251 (MQ=255)          gCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACAttt                                                                                                                                                                                                                    >  1:45289/1‑251 (MQ=255)                            agTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCg                                                                                                                                                                                                  >  2:126295/1‑251 (MQ=255)                                agaagaATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCg                                                                                                                                                                                              >  1:88481/1‑251 (MQ=255)                                    gaaTGCATCT‑ACGCAACAGAAACTAGCGATTTCACATAACTACAAGATATCAATTTCCTCTCACCTTAAAAAATCTCAATCGAGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCAGGGCGACGCCGCATGGATACCTTTGGAGTGATTATTATTGCGACAGCCAc                                                                                                                                                                                          >  2:383755/1‑251 (MQ=255)                                                    aCAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg                                                                                                                                                                            <  2:45289/250‑1 (MQ=255)                                                        aaaaTAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCa                                                                                                                                                                       <  1:379568/251‑1 (MQ=255)                                                                    tCACATATCTACAAGTTATCAATTTCCTCTCACCTGAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTGGGCATTACTGCCGAAGCCAGGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                                                                                                                                                     <  2:140608/225‑1 (MQ=255)                                                                    tCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                                                                                                                                                     >  1:140608/1‑225 (MQ=255)                                                                          aaCTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTctgct                                                                                                                                                      <  2:120856/250‑1 (MQ=255)                                                                                     aTCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATcc                                                                                                                                          >  1:243762/1‑251 (MQ=255)                                                                                         aTTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTc                                                                                                                                      >  1:436836/1‑251 (MQ=255)                                                                                           ttCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAATTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGATGTTACACATTTTGTATATGGTTGGCATTACTGCCGAAGCCATGACGGGAGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGACACCGCAATTGTCGGCGGGTCGGTGCTGGTTATTCTTCTGGGCACCTATACGCTCgg                                                                                                                                    >  2:141833/1‑251 (MQ=255)                                                                                               cctcacctTAAAAAATCGCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTTCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCTCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGAGATTCTGCTGGGCCACTATCCGCTCGGCTgg                                                                                                                                <  2:415977/251‑1 (MQ=255)                                                                                                 cctcctTAAAAAATCTCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTCCTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGGCGTCGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGt                                                                                                                              <  2:126550/248‑1 (MQ=255)                                                                                                  cacctTAAAAAATCTCAAGCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTgg                                                                                                                                <  1:219269/246‑1 (MQ=255)                                                                                                   accttccaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCCCCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTc                                                                                                                             >  2:182965/8‑250 (MQ=255)                                                                                                        aaaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAacac                                                                                                                       <  1:232032/251‑1 (MQ=255)                                                                                                                     cGTGACAATACGCACAAATCGCTACCCTGCCAGACCGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATAtgtgc                                                                                                          >  1:102221/1‑250 (MQ=255)                                                                                                                        gACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGAtt                                                                                                        <  1:18705/250‑1 (MQ=255)                                                                                                                                gcACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCt                                                                                               >  2:324878/1‑251 (MQ=255)                                                                                                                                                   cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGCATATGTGATTATCGTCGCTACCGCCGCAGTACTcaccc                                                                            >  2:88530/1‑250 (MQ=255)                                                                                                                                                                         cacagggtgatACACATTTGGTATTGGGTTGGCATTCCGGCCGAAGCCATGGGGGGGGCTCTGGCGGCCGGGCGACGCCGCAGGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATAGGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGcc                                                      <  2:102221/241‑1 (MQ=255)                                                                                                                                                                                  ttACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTgcg                                             >  2:111470/1‑251 (MQ=255)                                                                                                                                                                                   tacacaaattGTATTTGAATGGCATTACTGCCGATGCCATGACGCGGGCTATGGAGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAACGCAAATGGCGGAGGGTCAGTGCGCGATATTCTGCTGAGCAACTATCCGCTCGGCTGGGTCCACCACCCGGCATATGTGCTTATCGTCGCTACCGCCGCAGTACTCACCCCTATCGTTGATCCTGTAATGCCTTACCTgcgc                                            >  2:196929/9‑251 (MQ=255)                                                                                                                                                                                                                   gAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCAc                                                                                                                                               <  2:137523/120‑1 (MQ=255)                                                                                                                                                                                                                   gAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCAc                                                                                                                                               >  1:137523/1‑120 (MQ=255)                                                                                                                                                                                                                             ggggggCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTAGCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTTGGACTGGTGGTCttt  <  1:182965/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  TTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTTT  >  NC_000913/3818704‑3819173

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