Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,737,495	G→A	*114* (TAG→TAA)	raiA →	cold shock protein associated with 30S ribosomal subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,737,495	0	G	A	100.0%	97.5 / NA	32	*114* (TAG→TAA)	raiA	cold shock protein associated with 30S ribosomal subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/18);  total (14/18)

ATGTACACCGCAATTAACGAATTGATCAACAAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATT  >  NC_000913/2737361‑2737627                                                                                                                                       |                                                                                                                                     atGTACACCGCAATTAACGAATTGATCAACAAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCt                                                                                                                                  <  2:186132/139‑1 (MQ=255)  tGTACACCGCAATTAACGAATTGATCAACAAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCtt                                                                                                                                 <  2:1236654/139‑1 (MQ=255)   gTACACCGCAATTAACGAATTGATCAACAAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCttt                                                                                                                                <  2:1504107/139‑1 (MQ=255)                              caaGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCtt                                                                                                     >  2:965090/1‑139 (MQ=255)                                   tGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATGTTGAGTGTATCGCCAACGCGCCTTCGGgc                                                                                                >  2:1131252/1‑139 (MQ=255)                                   tGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGgc                                                                                                >  1:843391/1‑139 (MQ=255)                                      aaCGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCttt                                                                                                                                >  2:638078/1‑104 (MQ=255)                                      aaCGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCttt                                                                                                                                <  1:638078/104‑1 (MQ=255)                                                cAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTttgttg                                                                                   >  2:1250620/1‑139 (MQ=255)                                                          cagcaCAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGATCAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGt                                                                                            <  1:44443/120‑1 (MQ=255)                                                          cagcaCAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGt                                                                                            >  2:44443/1‑120 (MQ=255)                                                            gcaCAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGAc                                                                                 >  2:856691/1‑129 (MQ=255)                                                            gcaCAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGAc                                                                                 <  1:856691/129‑1 (MQ=255)                                                             cacaAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAAc                                                                      <  2:1622220/139‑1 (MQ=255)                                                                     cGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTCTTTTGTTGACAGCGTGAAAACAGTACggg                                                              <  1:965090/139‑1 (MQ=255)                                                                      gAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCG‑TTTTTGTTGACAGCGTGAAAACAGTACGGgta                                                            <  2:843391/139‑1 (MQ=255)                                                                       aaGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCtt                                                                                                                                 <  2:603697/70‑1 (MQ=255)                                                                       aaGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCtt                                                                                                                                 >  1:603697/1‑70 (MQ=255)                                                                       aaGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGaaaa                                                                       <  2:1207108/128‑1 (MQ=255)                                                                       aaGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGaaaa                                                                       >  1:1207108/1‑128 (MQ=255)                                                                              tcgtGTCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGAc                                                                                 >  1:1021147/1‑111 (MQ=255)                                                                              tcgtGTCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTATGTTGAc                                                                                 <  2:1021147/111‑1 (MQ=255)                                                                                  gCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAg                                                 <  1:1250620/139‑1 (MQ=255)                                                                                             ggTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACggg                                                              >  1:1353703/1‑115 (MQ=255)                                                                                             ggTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACggg                                                              <  2:1353703/115‑1 (MQ=255)                                                                                              gTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACggg                                                              >  2:1059/1‑114 (MQ=255)                                                                                              gTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACggg                                                              <  1:1059/114‑1 (MQ=255)                                                                                                                 tcgAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTttgtt                                                                                    >  2:1592509/1‑73 (MQ=255)                                                                                                                 tcgAAGAAGTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTttgtt                                                                                    <  1:1592509/73‑1 (MQ=255)                                                                                                                       aagTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAAc                                                                      <  2:1329539/81‑1 (MQ=255)                                                                                                                       aagTTGAAGAAGAGTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAAc                                                                      >  1:1329539/1‑81 (MQ=255)                                                                                                                                 agagTAATCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCAtt  <  1:1131252/139‑1 (MQ=255)                                                                                                                                   agtaatCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCa    >  2:1484270/1‑135 (MQ=255)                                                                                                                                   agtaatCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCa    <  1:1484270/135‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     ATGTACACCGCAATTAACGAATTGATCAACAAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATT  >  NC_000913/2737361‑2737627

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