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%	193.6 / NA	54	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 (27/27);  total (27/27)

CAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCG  >  NC_000913/4609203‑4609454                                                                                                                                |                                                                                                                               cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                   <  1:200643/139‑6 (MQ=255)       cGGGCTTGCAACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGt                                                                                                              >  1:1541721/1‑139 (MQ=255)                   ctgctgGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCa                                                                                                  >  2:400934/1‑139 (MQ=255)                     gctggctgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cc                                                                                                                             <  2:581001/109‑1 (MQ=255)                     gctggctgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cc                                                                                                                             >  1:581001/1‑109 (MQ=255)                            gAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                   >  2:993657/1‑107 (MQ=255)                            gAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                   <  1:993657/112‑6 (MQ=255)                                  acacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaa                                                                                   <  1:1087118/139‑1 (MQ=255)                                   cacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCACAATTTTTCGTTTTcaac                                                                                  >  2:19988/1‑139 (MQ=255)                                    acCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaaca                                                                                 >  1:111058/1‑139 (MQ=255)                                                aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCa                                                                                                  >  2:918412/1‑110 (MQ=255)                                                aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCa                                                                                                  <  1:918412/110‑1 (MQ=255)                                                 aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    <  1:813458/139‑1 (MQ=255)                                                  ccAGTAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                                                                   >  1:1065893/1‑139 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    <  2:549686/128‑1 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    >  1:549686/1‑128 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                         >  2:1565963/1‑139 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                         >  1:616308/1‑139 (MQ=255)                                                             cgcgCCCACCTGGACCGTCTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                                                        >  2:1107022/1‑139 (MQ=255)                                                                gcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaaga                                                     <  1:1565963/139‑1 (MQ=255)                                                                   cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  1:286871/139‑1 (MQ=255)                                                                   cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  2:1065893/139‑1 (MQ=255)                                                                     ccTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCACAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaa                                                <  1:19988/139‑1 (MQ=255)                                                                         gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCCGTTGATGGGTAAAATAGc                                                                                                      >  1:506931/1‑82 (MQ=255)                                                                         gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCCGTTGATGGGTAAAATAGc                                                                                                      <  2:506931/82‑1 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        >  2:8677/1‑106 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        <  1:8677/106‑1 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCC‑CCCCCCCGTTGATGGGTAAAATAGcc                                                                                                     <  1:404417/77‑1 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCC‑CCCCCCCGTTGATGGGTAAAATAGcc                                                                                                     >  2:404417/1‑77 (MQ=255)                                                                                tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa                                                 >  1:1113907/1‑127 (MQ=255)                                                                                tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa                                                 <  2:1113907/127‑1 (MQ=255)                                                                                tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     <  1:1107022/139‑1 (MQ=255)                                                                                tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     <  2:1541721/139‑1 (MQ=255)                                                                                        gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                             <  2:616308/139‑1 (MQ=255)                                                                                              aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     <  1:1807136/125‑1 (MQ=255)                                                                                              aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     >  2:1807136/1‑125 (MQ=255)                                                                                                 ggAGCAGCTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                    <  2:111058/139‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGc                                                               >  1:186570/1‑89 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGc                                                               <  2:186570/89‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGAc                                      >  2:56995/1‑114 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGAc                                      <  1:56995/114‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                    >  1:1513060/1‑132 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                    <  2:1513060/132‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  1:1807379/1‑139 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  1:879398/1‑139 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  2:550329/1‑139 (MQ=255)                                                                                                          ccTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                          >  2:1722361/1‑92 (MQ=255)                                                                                                          ccTGTGTAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                          <  1:1722361/92‑1 (MQ=255)                                                                                                             gtgtAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGaa        >  2:1154710/1‑139 (MQ=255)                                                                                                               gtAAACACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAgc      <  1:550329/139‑1 (MQ=255)                                                                                                                  aaCACTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgcc   >  1:768367/1‑139 (MQ=255)                                                                                                                   acacTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCTCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgccg  <  1:1154710/139‑1 (MQ=255)                                                                                                                   acacTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgccg  <  2:1667643/139‑1 (MQ=255)                                                                                                                     acTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                         <  2:1579392/82‑1 (MQ=255)                                                                                                                     acTGATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                         >  1:1579392/1‑82 (MQ=255)                                                                                                                        gATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    >  1:273701/1‑68 (MQ=255)                                                                                                                        gATTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    <  2:273701/68‑1 (MQ=255)                                                                                                                         aTTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                              >  2:567389/1‑73 (MQ=255)                                                                                                                         aTTGCCTC‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                              <  1:567389/73‑1 (MQ=255)                                                                                                                                |                                                                                                                               CAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCG  >  NC_000913/4609203‑4609454

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