Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,533,223	A→G	intergenic (‑40/+16)	yjhP ← / ← yjhQ	putative methyltransferase/GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,533,223	0	A	G	100.0%	192.1 / NA	50	intergenic (‑40/+16)	yjhP/yjhQ	putative methyltransferase/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/21);  total (29/21)

GCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGCTCCTCCGGGATGGGATA  >  NC_000913/4533094‑4533358                                                                                                                                  |                                                                                                                                        gcgACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGAt                                                                                                                                <  2:196415/139‑1 (MQ=255)   gACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTg                                                                                                                              <  1:262061/139‑1 (MQ=255)               cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGcaacaa                                                                                                                                   >  1:294650/1‑122 (MQ=255)               cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGcaacaa                                                                                                                                   <  2:294650/122‑1 (MQ=255)               cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCa                                                                                                                  >  1:752011/1‑139 (MQ=255)               cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATGGGGTTAAGCCTAAAAAAGAACGCAACAAGATTGTTATTCACGTCa                                                                                                                  >  2:274630/1‑139 (MQ=255)                        ctTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCgg                                                                                                         >  1:482576/1‑139 (MQ=255)                           cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGttt                                                                                                      >  2:21058/1‑139 (MQ=255)                              tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGa                                                                                                                                 >  1:696641/1‑109 (MQ=255)                              tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGa                                                                                                                                 <  2:696641/109‑1 (MQ=255)                              tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcat                                                                                                   >  2:203726/1‑139 (MQ=255)                                aaCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACgcaac                                                                                                                                     <  2:179899/103‑1 (MQ=255)                                aaCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACgcaac                                                                                                                                     >  1:179899/1‑103 (MQ=255)                                              gCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCg                                                                                   >  1:234412/1‑139 (MQ=255)                                               cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCt                                                                                                             >  1:704708/1‑112 (MQ=255)                                               cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCt                                                                                                             <  2:704708/112‑1 (MQ=255)                                               cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCa                                                                                      >  1:315264/1‑135 (MQ=255)                                               cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCa                                                                                      <  2:315264/135‑1 (MQ=255)                                               cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCgg                                                                                  >  1:557853/1‑139 (MQ=255)                                                    ttCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgtg                                                                             >  2:120509/1‑139 (MQ=255)                                                     tCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                                                            >  2:695851/1‑139 (MQ=255)                                                                 aTGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTc                                                                >  2:529234/1‑139 (MQ=255)                                                                     tAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATgg                                                            >  1:398413/1‑139 (MQ=255)                                                                       aaaaTACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGc                                                          <  1:203726/139‑1 (MQ=255)                                                                           tACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                                                      <  2:482576/139‑1 (MQ=255)                                                                             cGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCg                                                    >  1:161419/1‑139 (MQ=255)                                                                             cGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCg                                                    >  1:396907/1‑139 (MQ=255)                                                                              gTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCgg                                                   >  2:520675/1‑139 (MQ=255)                                                                               tGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                                                  >  1:382943/1‑139 (MQ=255)                                                                                   atatCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagt                                              <  2:557853/139‑1 (MQ=255)                                                                                      tCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAg                                                                                                                 <  2:434358/69‑1 (MQ=255)                                                                                      tCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAg                                                                                                                 >  1:434358/1‑69 (MQ=255)                                                                                         attattTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc                                        >  1:551318/1‑139 (MQ=255)                                                                                                aaTTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCa                                                                                      >  2:598494/1‑86 (MQ=255)                                                                                                aaTTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCa                                                                                      <  1:598494/86‑1 (MQ=255)                                                                                                aaTTCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATcc                                 >  2:666407/1‑139 (MQ=255)                                                                                                  ttCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcag                               >  1:484454/1‑139 (MQ=255)                                                                                                   tCCTGTTATTGGGTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcagc                              <  2:752011/139‑1 (MQ=255)                                                                                                    tcctttttaaggTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                                                      >  1:378377/11‑114 (MQ=255)                                                                                                    tcctttttaaggTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                                                      <  2:378377/104‑1 (MQ=255)                                                                                                             gggTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                                                 <  2:94371/110‑1 (MQ=255)                                                                                                             gggTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                                                 >  1:94371/1‑110 (MQ=255)                                                                                                             gggTTAAGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTtgtg                    <  2:396907/139‑1 (MQ=255)                                                                                                                  aaGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                                                            >  2:321309/1‑78 (MQ=255)                                                                                                                  aaGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                                                            <  1:321309/78‑1 (MQ=255)                                                                                                                  aaGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                                              >  1:506514/1‑92 (MQ=255)                                                                                                                  aaGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTATGTGctcct               >  1:501032/1‑139 (MQ=255)                                                                                                                  aaGCCTTAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAAGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                                              <  2:506514/92‑1 (MQ=255)                                                                                                                        tAAAAAGGACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGCTCCTCCggga         <  1:695851/139‑1 (MQ=255)                                                                                                                        tAAAAAGGACGCAACAAGATTGTTATTCAAGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATAGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGCTCCTCCggga         <  1:529234/139‑1 (MQ=255)                                                                                                                               gACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGc                   <  2:376089/122‑1 (MQ=255)                                                                                                                               gACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGc                   >  1:376089/1‑122 (MQ=255)                                                                                                                               gACGCAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGCTCCTCCGGGATGGGATa  <  2:234412/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        GCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAGGCTTTGTGCTCCTCCGGGATGGGATA  >  NC_000913/4533094‑4533358

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