Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,609,336	(C)7→8	intergenic (+13/‑78)	yjjG → / → prfC	pyrimidine nucleotidase/peptide chain release factor RF3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,609,329	1	.	C	100.0%	111.8 / NA	32	intergenic (+6/‑85)	yjjG/prfC	pyrimidine nucleotidase/peptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/11);  total (21/11)

GTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTT  >  NC_000913/4609196‑4609461                                                                                                                                       |                                                                                                                                      gTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cccccc                                                                                                                                >  2:459333/1‑140 (MQ=255)      aTCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgtt                                                                                                                           >  2:779939/1‑136 (MQ=255)       tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                          >  1:792247/1‑135 (MQ=255)       tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                          <  1:779939/140‑6 (MQ=255)           cGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATgg                                                                                                                       >  1:781884/1‑140 (MQ=255)              cGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGCTGGGTa                                                                                                                    >  1:884085/1‑140 (MQ=255)                                 ctgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCg                                                                                                 <  1:459333/140‑1 (MQ=255)                                              aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                                                                    >  2:71206/1‑140 (MQ=255)                                               tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAGCACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                                                                   >  2:936543/1‑140 (MQ=255)                                               tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                                                                   >  1:461280/1‑140 (MQ=255)                                                cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcg                                                                                  >  2:369741/1‑140 (MQ=255)                                                      cAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                            <  2:461280/140‑1 (MQ=255)                                                       aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                           >  1:125910/1‑140 (MQ=255)                                                             aaGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                     <  1:369741/140‑1 (MQ=255)                                                                    cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                              >  1:383271/1‑140 (MQ=255)                                                                    cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                              >  2:900707/1‑140 (MQ=255)                                                                    cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                              >  2:270315/1‑140 (MQ=255)                                                                     gcgcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAGCACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaag                                                             <  1:936543/140‑1 (MQ=255)                                                                         ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                         <  2:781884/140‑1 (MQ=255)                                                                           aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaa                                                       >  1:889028/1‑140 (MQ=255)                                                                                 aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTa                                                 >  1:270973/1‑140 (MQ=255)                                                                                 aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTa                                                 >  2:276900/1‑140 (MQ=255)                                                                                    gTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                              >  1:101103/1‑140 (MQ=255)                                                                                     tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGAc                                             >  2:2154/1‑140 (MQ=255)                                                                                     tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGAc                                             >  1:513310/1‑140 (MQ=255)                                                                                      ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                            <  1:900707/140‑1 (MQ=255)                                                                                      ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                            <  1:71206/140‑1 (MQ=255)                                                                                        cttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGtt                                          <  1:2154/140‑1 (MQ=255)                                                                                                            cagcTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAgg                      <  2:125910/140‑1 (MQ=255)                                                                                                               cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTgg                   >  2:61046/1‑140 (MQ=255)                                                                                                               cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTgg                   >  1:564335/1‑140 (MQ=255)                                                                                                                         aaCACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgccg         >  1:476767/1‑140 (MQ=255)                                                                                                                            acTGATTGCCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCa       >  1:279286/1‑140 (MQ=255)                                                                                                                               gATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACttt   <  1:270315/140‑1 (MQ=255)                                                                                                                               gATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACttt   >  1:343401/1‑140 (MQ=255)                                                                                                                                aTTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtttt  <  2:884085/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      GTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTT  >  NC_000913/4609196‑4609461

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