Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,473,615	C→T	W156* (TGG→TAG)	rpsG ←	30S ribosomal subunit protein S7

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,473,615	0	C	T	100.0%	45.9 / NA	21	W156* (TGG→TAG)	rpsG	30S ribosomal subunit protein S7
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/12);  total (9/12)

GTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCC  >  NC_000913/3473499‑3473734                                                                                                                     |                                                                                                                        tcgttgtGCGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGtgtg                                                                                                             <  2:352126/127‑1 (MQ=255)        aCGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAATCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGc                                                                                                      <  1:325367/129‑1 (MQ=255)         tgatccATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCTCAAAGCGGGCTGCTTATTGGAAGCGCCCGTCTGTTGATTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGcc                                                                                                     <  1:407497/125‑1 (MQ=255)         cGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAATCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGcc                                                                                                     <  1:319285/129‑1 (MQ=255)                  tGTTTCCTCGTTTATCTTTTAGGCATTCAATTTAAGTAGCCCAACGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGCCTAAACCTCCGACGGGATAACTACCGGTAGTGTGCGAACGCCTCGTTGGCt                                                                                            >  2:629010/1‑129 (MQ=255)                    tttCCTCGTTTATGCTTTAGGGCGTCAGTTTAAGTAGCCCAGAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAGCTCCGAAGGGATAACTAACGGTAGTGGGCGAACGCCTTGTTGGCTTc                                                                                          <  2:610212/129‑1 (MQ=255)                    cttCCTCGTTTATCTTTTCGGCTTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACTGTAGTGTGCGAACGCCTTGTTGGCTTc                                                                                          <  2:37000/128‑1 (MQ=255)                          cGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCATCGGCCaa                                                                                    >  2:493825/1‑128 (MQ=255)                                          tttCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTATCGTTATTTTTTGTTTTTCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCa                                                                    <  1:11607/128‑1 (MQ=255)                                                tttAAGTAGCCCAAAGCGGGCTGCTTACTGTACGCGCCCGCCTGGTGACTAAAACTCCGAAGGTATATCTATCTGTAGTGTTTTTACGTTTTGTTGGCTTAGGCCATACGGTGAACGTCTTCACGTTTc                                                              <  1:96426/129‑1 (MQ=255)                                                      tAGCCCAAAGCGGGCTGCTTACTGGAAGCCCCCGCCTGATAACTAAAACTCCGAAGGAATAACTAACAGTAGTGTGCGAACGCCTTGTTGACTTCGGCCATACGGTGAAAGTCTTACAGTTACTTAACt                                                        >  1:649015/1‑129 (MQ=255)                                                            aaaGCCGGCTGCCTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTTCTTTTGTTTTTTTTGCTTTGGCTATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTa                                                  <  1:20634/129‑1 (MQ=255)                                                              aGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTAcc                                                >  1:323910/1‑129 (MQ=255)                                                                        ttACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGACGGGATAACTAAAAAAAAAATAAAAACACAAAAATGACTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCt                                      >  1:589568/1‑129 (MQ=255)                                                                        ttACTGGAAGCGCCCGCCTGATGACTAAAACTCCGAAAGGATAACTAACGGTAGTGTGAGAACGCCTTATTGACTTCGAACATACGGTAAACGTATACACGATTCTTAACTGCAATAACATTGTTTTCt                                      >  2:609123/1‑129 (MQ=255)                                                                              gAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTTAACGTTTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTgcagca                                <  1:590439/129‑1 (MQ=255)                                                                                  cgcccgccTGGTGACTAAAACTCCGACGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAg                            <  1:180912/129‑1 (MQ=255)                                                                                              gACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACATTTGTTTTCTGCAGCATCACAAAGttcgttcg                >  2:463147/1‑129 (MQ=255)                                                                                                    aaCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGAATCAGAAAGTTCGTTCGCCAgac          >  2:74686/1‑127 (MQ=255)                                                                                                          gAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTAAACGAATTCGAGGATGTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCagag    >  2:360059/1‑129 (MQ=255)                                                                                                            aGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACTTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGcc  <  1:99249/129‑1 (MQ=255)                                                                                                                     |                                                                                                                        GTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCC  >  NC_000913/3473499‑3473734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑