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%	124.1 / NA	36	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 (18/18);  total (18/18)

CAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609203‑4609459                                                                                                                                |                                                                                                                                   cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  2:426152/139‑6 (MQ=255)  aaCGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑cccccccccgttg                                                                                                                       <  2:1097573/139‑7 (MQ=255)   aCGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGAt                                                                                                                       <  1:451640/139‑1 (MQ=255)       cGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGt                                                                                                                   >  2:334735/1‑139 (MQ=255)                                      cATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgtt                                                                                                                         >  1:685352/1‑97 (MQ=255)                                      cATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgtt                                                                                                                         <  2:685352/101‑5 (MQ=255)                                                aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                          >  1:1422183/1‑139 (MQ=255)                                                 aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑c                                                                                                                                   <  1:1108005/80‑1 (MQ=255)                                                 aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑c                                                                                                                                   >  2:1108005/1‑80 (MQ=255)                                                         gCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTa                                                                 <  1:381525/139‑1 (MQ=255)                                                           aTCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                               <  1:14312/139‑1 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                          >  2:589492/1‑127 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                          <  1:589492/127‑1 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   >  2:993717/1‑134 (MQ=255)                                                            tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   <  1:993717/134‑1 (MQ=255)                                                                gcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaaga                                                          <  1:334735/139‑1 (MQ=255)                                                                 cccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaa                                                         >  2:398182/1‑139 (MQ=255)                                                                  ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaag                                                        <  1:398182/139‑1 (MQ=255)                                                                  ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        >  2:1250017/1‑69 (MQ=255)                                                                  ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  1:1250017/74‑6 (MQ=255)                                                                   cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAACTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                       <  1:164643/139‑1 (MQ=255)                                                                     ccTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaa                                                     <  2:1422183/139‑1 (MQ=255)                                                                          aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                                                                       <  1:1279211/100‑1 (MQ=255)                                                                          aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                                                                       >  2:1279211/1‑100 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            >  1:745218/1‑139 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        >  2:240564/1‑57 (MQ=255)                                                                              tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  1:240564/62‑6 (MQ=255)                                                                                tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                          <  1:656154/139‑1 (MQ=255)                                                                                        gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcgcg                                                                           >  2:1304432/1‑98 (MQ=255)                                                                                        gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcgcg                                                                           <  1:1304432/98‑1 (MQ=255)                                                                                        gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        <  2:761441/52‑6 (MQ=255)                                                                                        gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑ccccccccgttg                                                                                                                        >  1:761441/1‑47 (MQ=255)                                                                                        gCACGAACTGGAGCAGCGCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTAATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                  >  1:922462/1‑139 (MQ=255)                                                                                              aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                          <  1:275454/125‑1 (MQ=255)                                                                                              aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                          >  2:275454/1‑125 (MQ=255)                                                                                                 ggAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                         <  2:745218/139‑1 (MQ=255)                                                                                                 ggAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                         >  1:231660/1‑139 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 <  1:402177/60‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 >  2:402177/1‑60 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTa                                                                 >  2:300570/1‑92 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTa                                                                 <  1:300570/92‑1 (MQ=255)                                                                                                        cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg                  >  2:112337/1‑139 (MQ=255)                                                                                                          ccTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                             >  1:794825/1‑78 (MQ=255)                                                                                                          ccTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                             <  2:794825/78‑1 (MQ=255)                                                                                                                       tGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   <  2:1505471/75‑1 (MQ=255)                                                                                                                       tGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                   >  1:1505471/1‑75 (MQ=255)                                                                                                                        gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACTAGCGCGGCGCGATGCCGCTTACTCAAGCAGAACGAATTATGACGTTGTCTCCTTATTTGCAAGAAGTGGCGAAGCGCCGCACtt  >  2:602191/1‑139 (MQ=255)                                                                                                                        gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtt  >  1:511954/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                   CAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609203‑4609459

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