Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,351,623	G→A	Y46Y (TAC→TAT)	ycjD ←	DUF559 domain‑containing protein YcjD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,351,623	0	G	A	100.0%	75.6 / NA	26	Y46Y (TAC→TAT)	ycjD	DUF559 domain‑containing protein YcjD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

GGGAGAGGGTGACCGGCGGTTCAGTTCCTGCAGAATATTCTCCAGCACCGTCTCCTCATTACAATCAATCTCGTTATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGC  >  NC_000913/1351410‑1351868                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gggAGAGGGTGACCGGCGGTTCAGTTCCTGCAGAATATTCTCCAGCACCGTCTCCTCATTACAATCAATCTCGTTATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGa                                                                                                                                                                                                                                  <  1:239160/235‑1 (MQ=255)                   gTTCAGTTCCTGCAGAATATTCTCCAGCGCCGTCTCCTCATTACAATCAATCTCGTGATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCGAGGGCGCCTTTAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTg                                                                                                                                                                                                                                   <  1:358670/216‑1 (MQ=255)                                                        tcATTACAATCAATCTCGTTATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCTAGCAGCAAGCAAAATCGAGAATATAGCTccc                                                                                                                                                                                                                                               <  1:166308/167‑1 (MQ=255)                                                                atcaatcTCGTTATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCTGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTc                                                                                                                                                   <  1:244051/251‑1 (MQ=255)                                                                        cGTTAGTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAg                                                                                                                                                                                            <  1:467592/202‑1 (MQ=255)                                                                               ccAGANNCGNNNCANNGTCCAGCCCTGCGACTCAAGCCANCTAGTGCGCCTGGAATCATAGGCAACTGCNANANCATGNNGCCNACCATCCAGCNCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCNCNCTGGATGTTGACgg                                                                                                                                                                                                                               >  2:377140/1‑160 (MQ=255)                                                                               ccAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTTCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACgg                                                                                                                                                                                                                               <  1:377140/160‑1 (MQ=255)                                                                                          gcaCGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTcccc                                                                                                                                                                                                                                              >  1:188881/1‑134 (MQ=255)                                                                                                  cAGCCCTGNGNNNCNNGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCNGCCCACCATCCAGCTNAANGACTNCACNCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGa                                                                                                                                                                                                                                  >  2:359113/1‑138 (MQ=255)                                                                                                  cAGCCCTGCTACTCAAGCCAGCTAGTGCGCCTGGAATCATAGTCAACTGCTAAATCATGCTGCCCACCATCTAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGa                                                                                                                                                                                                                                  <  1:359113/138‑1 (MQ=255)                                                                                                      ncTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCAcc                                                                                                             >  1:133898/2‑251 (MQ=255)                                                                                                      ccTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGc                                                                                                                                                                                     >  1:415080/1‑179 (MQ=255)                                                                                                                                   gAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGCCTACACGCGCCGAGCAGCAAGCAAAATCGAGCATATAGCTCCCCCCTGGATGTTGACGGCGAAATTTGAAATCCCTAAAACGTCGGCTGCGAAg                                                                                                                                                                                            >  1:474981/1‑143 (MQ=255)                                                                                                                                                  gCTAAATCATGNTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCNNANNNGNGAATATAg                                                                                                                                                                                                                                                    <  2:358824/72‑1 (MQ=255)                                                                                                                                                  gCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGCGAATATAg                                                                                                                                                                                                                                                    >  1:358824/1‑72 (MQ=255)                                                                                                                                                                            cTCAACGACTACACGCGACGAGCAGCAAGCAAAATCGAGAATATAGCTCCAAACTGGATGTTGACGGc                                                                                                                                                                                                                              >  1:450667/1‑68 (MQ=255)                                                                                                                                                                              cAACGACTACACGCGCCGCGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGCTATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACACGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                     >  1:168369/1‑251 (MQ=255)                                                                                                                                                                              cAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGt                                                                                                                                                                         <  1:231055/119‑1 (MQ=255)                                                                                                                                                                                       cacGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAACTTTGACATCACTATAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAg                                                                                                                                                               >  1:423362/1‑120 (MQ=255)                                                                                                                                                                                       cacGCGCCGAGCAGCAAGCAAAAGCGAGAATATAGCTCCACACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGATGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCGTAAAGCACGAGCATTTGATGTAATTTTATCGATCACCTCATTCTGACAAGATTTAATCTGTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATGCAGAAAAGTAGCGATCTTTACg                            <  1:183121/251‑1 (MQ=255)                                                                                                                                                                                        acGCGCCGAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGc                           >  1:109609/1‑251 (MQ=255)                                                                                                                                                                                            gcCGAGCAGCAAGCAAA‑TCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTACTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                      >  1:324610/1‑251 (MQ=255)                                                                                                                                                                                               gAGCAGCAAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCAt                    >  1:65790/1‑251 (MQ=255)                                                                                                                                                                                                cgcagcaAGCAAAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATc                   <  1:137840/250‑1 (MQ=255)                                                                                                                                                                                                       agcaaAATCGAGAATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATc                                                                                                                                            <  1:78594/123‑1 (MQ=255)                                                                                                                                                                                                                 agaATATAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCTCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCTTCCGCCGTCAATGCCTCGc  <  1:385627/251‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GGGAGAGGGTGACCGGCGGTTCAGTTCCTGCAGAATATTCTCCAGCACCGTCTCCTCATTACAATCAATCTCGTTATTCCAGAAACGCAGCACGGTCCAGCCCTGCGACTCAAGCCAGCTAGTGCGCCTGGAATCATAGGCAACTGCTAAATCATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGC  >  NC_000913/1351410‑1351868

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