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

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%	98.9 / NA	29	intergenic (+6/‑85)	yjjG/prfC	dUMP phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/15);  total (14/15)

TCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609193‑4609459                                                                                                                                          |                                                                                                                                   tcGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑cc                                                                                                                                  >  1:236249/1‑139 (MQ=255)   ggTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTAGAGCAGCTCCTGTGTAAACACTGATTGCCT‑cccc                                                                                                                                >  1:305794/1‑139 (MQ=255)          tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgtt                                                                                                                         >  2:938212/1‑135 (MQ=255)           cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  2:298573/139‑6 (MQ=255)           cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTAGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  2:305794/139‑6 (MQ=255)                               gctggctgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTACCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAt                                                                                                     >  1:649062/1‑139 (MQ=255)                                       aaTGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                                                             <  1:938212/139‑1 (MQ=255)                                        aTGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGttt                                                                                            <  2:649062/139‑1 (MQ=255)                                                 aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 >  1:45719/1‑125 (MQ=255)                                                 aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 <  2:45719/125‑1 (MQ=255)                                                       aGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                             <  1:191868/139‑1 (MQ=255)                                                            ccAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTaaaa                                                                                                               >  1:138526/1‑100 (MQ=255)                                                            ccAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTaaaa                                                                                                               <  2:138526/100‑1 (MQ=255)                                                                aaGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGc                                                                    <  2:236249/139‑1 (MQ=255)                                                                 aGGCATCGCGCCCACCTGGACTGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   <  1:631131/139‑1 (MQ=255)                                                                    cATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 <  1:203900/106‑1 (MQ=255)                                                                    cATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 >  2:203900/1‑106 (MQ=255)                                                                       cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  2:134454/79‑6 (MQ=255)                                                                       cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        >  1:134454/1‑74 (MQ=255)                                                                        gcgcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                            <  1:795643/139‑1 (MQ=255)                                                                                   gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCg                                                                     <  1:338005/119‑1 (MQ=255)                                                                                   gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCg                                                                     >  2:338005/1‑119 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGa                                                                                                                        <  1:482495/63‑1 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGa                                                                                                                        >  2:482495/1‑63 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   <  1:550645/116‑1 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   <  2:309124/116‑1 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   >  1:309124/1‑116 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   >  2:550645/1‑116 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            >  2:125499/1‑139 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  1:218441/62‑6 (MQ=255)                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        >  2:218441/1‑57 (MQ=255)                                                                                          tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            >  1:50723/1‑137 (MQ=255)                                                                                          tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            <  2:50723/137‑1 (MQ=255)                                                                                                  gCACGAACTGGAGCAGCTCCTGTTTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                  <  1:623695/139‑1 (MQ=255)                                                                                                  gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                  >  1:266187/1‑139 (MQ=255)                                                                                                                  cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg                  >  1:605581/1‑139 (MQ=255)                                                                                                                                  gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtt  >  2:30504/1‑139 (MQ=255)                                                                                                                                  gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtt  >  2:1004187/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                   TCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609193‑4609459

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: ‑