Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,609,352	A→G	intergenic (+29/‑62)	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,352	0	A	G	100.0%	154.1 / NA	43	intergenic (+29/‑62)	yjjG/prfC	pyrimidine nucleotidase/peptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/17);  total (26/17)

GCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCATTATTTCTCACCCGGA  >  NC_000913/4609216‑4609481                                                                                                                                         |                                                                                                                                  gCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGcc                                                                                                                                >  2:155856/1‑140 (MQ=255)             ctgAATGCACACCATCGCGAGC‑ACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGtt                                                                                                                   <  1:304991/140‑1 (MQ=255)              tgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGtt                                                                                                                   <  1:187510/140‑1 (MQ=255)                      cacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTcaacaa                                                                                                           >  1:857012/1‑140 (MQ=255)                          aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAgcgc                                                                                                       >  1:870352/1‑140 (MQ=255)                           tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAgcgcg                                                                                                      >  2:96696/1‑140 (MQ=255)                            cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAgcgcgg                                                                                                     >  2:852233/1‑140 (MQ=255)                            cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAgcgcgg                                                                                                     >  1:382046/1‑140 (MQ=255)                                 gCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                                                >  1:23624/1‑140 (MQ=255)                                   aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                                                                                              >  2:676903/1‑140 (MQ=255)                                       aGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGc                                                                                          <  1:130760/140‑1 (MQ=255)                                        gAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                                         <  2:870352/140‑1 (MQ=255)                                           ggCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTAc                                                                                      >  1:888806/1‑140 (MQ=255)                                             cATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                                                    <  2:480182/140‑1 (MQ=255)                                               tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                                                  >  1:980749/1‑140 (MQ=255)                                                   gcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaag                                                                              >  1:23016/1‑140 (MQ=255)                                                             aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt                                                                    >  1:887939/1‑140 (MQ=255)                                                                 tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                                                >  1:113771/1‑140 (MQ=255)                                                                 tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                                                >  1:800149/1‑140 (MQ=255)                                                                 tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                                                <  2:23016/140‑1 (MQ=255)                                                                   tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGtt                                                              >  2:811532/1‑140 (MQ=255)                                                                         ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                                        <  2:951464/140‑1 (MQ=255)                                                                         ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                                        <  2:782070/140‑1 (MQ=255)                                                                           gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                                        >  1:655551/1‑138 (MQ=255)                                                                           gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                                        <  2:655551/138‑1 (MQ=255)                                                                           gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCtt                                                      >  2:407991/1‑140 (MQ=255)                                                                                  cTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                                               <  2:272934/140‑1 (MQ=255)                                                                                  cTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                                               <  1:407991/140‑1 (MQ=255)                                                                                       gcagcTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAgg                                          <  2:10215/140‑1 (MQ=255)                                                                                           cTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGc                                      >  1:434734/1‑140 (MQ=255)                                                                                           cTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGc                                      >  1:728202/1‑140 (MQ=255)                                                                                           cTCCTGTGTAAACACTGATTGCCTCCCCACCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGc                                      >  1:442979/1‑140 (MQ=255)                                                                                                        acTGATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACt                         >  2:468712/1‑140 (MQ=255)                                                                                                           gATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtttt                      >  1:982264/1‑140 (MQ=255)                                                                                                           gATTGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGAAGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtttt                      <  2:888806/140‑1 (MQ=255)                                                                                                             ttGCCTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGc                    >  2:602101/1‑140 (MQ=255)                                                                                                                ccTCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTACTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCat                 <  1:155856/140‑1 (MQ=255)                                                                                                                  tCCCCCCCGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCatta               >  1:94416/1‑140 (MQ=255)                                                                                                                    ccccccGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCat                 <  2:31378/136‑1 (MQ=255)                                                                                                                    ccccccGTTGATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCat                 >  1:31378/1‑136 (MQ=255)                                                                                                                             gATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCATTATTTCTCACCCg    >  2:544028/1‑140 (MQ=255)                                                                                                                             gATGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCATTATTTCTCACCCg    <  1:493913/140‑1 (MQ=255)                                                                                                                               tGGGTAAAATGGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCATTATTTCTCACCCGGa  <  1:544028/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTTTTGCCATTATTTCTCACCCGGA  >  NC_000913/4609216‑4609481

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