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%	94.5 / NA	27	H57N (CAC→AAC)	priB	primosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/15);  total (12/15)

GTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAC  >  NC_000913/4425446‑4425904                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         gTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGaaaa                                                                                                                                                                                                                   >  2:42405/1‑250 (MQ=255)                atTTCGCAAACGAAACCGCTGATGATGCTGCAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTcac                                                                                                                                                                                                   >  1:184504/1‑251 (MQ=255)                     gCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCGCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATTCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGtat                                                                                                                                                                                              <  1:12062/251‑1 (MQ=255)                                    gatgatGCTTAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAATGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCGTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTc                                                                                                                                                                                <  2:168980/250‑1 (MQ=255)                                        atgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATa                                                                                                                                                                           <  2:4061/251‑1 (MQ=255)                                                gCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTc                                                                                                                                                                    >  1:53240/1‑250 (MQ=255)                                                                          gatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCa                                                                                                                                          <  1:104388/250‑1 (MQ=255)                                                                                cAACCGTCTGGTGTTGTCCGGCACCGTGTTCAGGTCTCCCCTTCGAAAGTTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTATCGGAAACGAAAAACAGGCCATTTCTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACg                                                                                                                                   <  2:218320/251‑1 (MQ=255)                                                                                        tGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGc                                                                                                                           >  1:21798/1‑251 (MQ=255)                                                                                               gtCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATgg                                                                                                                    >  2:95502/1‑251 (MQ=255)                                                                                               gtCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATgg                                                                                                                    >  1:215267/1‑251 (MQ=255)                                                                                                 ccGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCGGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCTTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtt                                                                                                                  <  1:130030/251‑1 (MQ=255)                                                                                                 ccGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGACCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCACGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTAATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtt                                                                                                                  >  2:193658/1‑251 (MQ=255)                                                                                                  cGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTTAGCAAAATGGttt                                                                                                                 >  2:169700/1‑251 (MQ=255)                                                                                                   ggCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtttt                                                                                                                <  1:94129/251‑1 (MQ=255)                                                                                                           gtgCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGcc                                                                                                         <  2:81586/250‑1 (MQ=255)                                                                                                             gCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAg                                                                                                      <  2:57087/251‑1 (MQ=255)                                                                                                                gggCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                     >  2:86845/1‑249 (MQ=255)                                                                                                                 ggCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGGTCATTTCATTCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                     >  2:4746/1‑248 (MQ=255)                                                                                                                         ttCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGattgaattg                                                                                          >  2:183985/1‑251 (MQ=255)                                                                                                                                   cAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTg                                                                                <  1:47996/251‑1 (MQ=255)                                                                                                                                   cAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTg                                                                                <  1:47999/251‑1 (MQ=255)                                                                                                                                       ccATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTGGCATAACCGTTCAGGGGTTGATTTCATGCCGCAAGGCAAAGAGCGGACTGATCAAAATGGTTTTGCATTCCGAGCAGATTGAATTGATAGATTCTGgaga                                                                            <  1:15401/251‑1 (MQ=255)                                                                                                                                            aGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGAGGGAATTGATAGATTCTGGAGACTAGc                                                                       <  2:215267/251‑1 (MQ=255)                                                                                                                                                            cAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTAtt                                                       >  2:100461/1‑251 (MQ=255)                                                                                                                                                                                              aaGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                     <  1:95502/251‑1 (MQ=255)                                                                                                                                                                                                                 ggCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAc  <  1:173202/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         GTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAC  >  NC_000913/4425446‑4425904

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