Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	576,410	G→A	pseudogene (416/1040 nt)	nmpC ←	putative outer membrane porin NmpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	576,410	0	G	A	100.0%	25.7 / NA	10	pseudogene (416/1040 nt)	nmpC	putative outer membrane porin NmpC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/6);  total (4/6)

CATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGG  >  NC_000913/576222‑576626                                                                                                                                                                                             |                                                                                                                                                                                                                         cATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGTATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAg                                                                                                                                                            <  1:344956/251‑1 (MQ=255)                                                            cTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCTAATTCTGTCAGTACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCg                                                                                                                                                   <  1:293298/201‑1 (MQ=255)                                                                           ccTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGCACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTgg                                                                                   >  1:324477/1‑250 (MQ=255)                                                                                gTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGGAAc                                                                             >  1:25099/1‑251 (MQ=255)                                                                                                  cAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACCCCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCAc                                                                                                                             >  1:64270/1‑185 (MQ=255)                                                                                                  cAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCgg                                                            <  1:28332/250‑1 (MQ=255)                                                                                                                             gttgttACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGTAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGAc                                <  1:235908/251‑1 (MQ=255)                                                                                                                                                          cACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTgg   <  1:166849/251‑1 (MQ=255)                                                                                                                                                           ccGTTGAGTCATGAACACGGCGGTTTGAGTCCAAATGTCACCACCGAATTCGGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGTAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGTAGGACCCGTGAGAGTCAGCGCGGTTGCGTGTGAATTCATATTCCCACTGACCGAATCCAGTCAGTTGATCGTTGATTTggg  <  1:290350/250‑1 (MQ=255)                                                                                                                                                                                       gTCCAAATGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGcc                                                                                                                >  1:266406/1‑113 (MQ=255)                                                                                                                                                                                             |                                                                                                                                                                                                                         CATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGTCAGTCCACGCACCGATGTCGTATGCTACACCGTAGTTACGGCCGTAATCGATGGAGCCGTAGTCACCGAATTTCAGGCCAGCGAAGGCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGG  >  NC_000913/576222‑576626

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