Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,425,688	C→A	H57N (CAC→AAC)	priB →	primosomal replication protein N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,425,688	0	C	A	100.0%	61.6 / NA	19	H57N (CAC→AAC)	priB	primosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/11);  total (8/11)

CCGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAA‑CCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAA  >  NC_000913/4425444‑4425929                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ccGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCGCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTGGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTGCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGaaa                                                                                                                                                                                                                                              <  2:168614/251‑1 (MQ=255)              gatgatTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTAc                                                                                                                                                                                                                                <  1:195751/251‑1 (MQ=255)                   tttCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTcac                                                                                                                                                                                                                            >  2:109906/1‑250 (MQ=255)                                      gatgatGCTGAAGCTGGGGATGCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCg                                                                                                                                                                                                        <  2:215058/251‑1 (MQ=255)                                               gAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCg                                                                                                                                                                                                <  2:55359/250‑1 (MQ=255)                                                                 gagTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATg                                                                                                                                                                             <  2:102336/251‑1 (MQ=255)                                                                         tCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGTCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACCAG‑CCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACCAgg                                                                                                                                                                     >  1:31742/1‑251 (MQ=255)                                                                                   aaCCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGTGTTCATTTCATGCCACAAGGCAAAGAACgg                                                                                                                                                           >  1:17649/1‑251 (MQ=255)                                                                                       gTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAATCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTg                                                                                                                                                       <  2:64828/251‑1 (MQ=255)                                                                                                                 aGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                                              >  2:13884/1‑250 (MQ=255)                                                                                                                       ccccTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGattg                                                                                                                        <  1:53860/250‑1 (MQ=255)                                                                                                                             cGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGattgaattga                                                                                                                  <  2:198573/250‑1 (MQ=255)                                                                                                                                       gTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTgg                                                                                                        >  1:27215/1‑250 (MQ=255)                                                                                                                                            tCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACt                                                                                                   >  2:33231/1‑250 (MQ=255)                                                                                                                                            gcAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTAAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTa                                                                                                  <  1:75226/250‑1 (MQ=255)                                                                                                                                                                                                  gCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGCTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAg                                            <  2:27215/251‑1 (MQ=255)                                                                                                                                                                                                   ccGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAgg                                           >  2:227209/1‑251 (MQ=255)                                                                                                                                                                                                               ggCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAag                                   >  2:11722/1‑247 (MQ=255)                                                                                                                                                                                                                                            gTTAGCGGAAACGAAAA‑ACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCTAGAGATCGACGATAAAGATATCGCTACGCTGaaaa  <  2:251575/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CCGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAA‑CCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAA  >  NC_000913/4425444‑4425929

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