Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,085	2 bp→TA	pseudogene (51‑52/213 nt)	renD →	protein RenD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,085	0	C	T	100.0%	34.3 / NA	11	pseudogene (51/213 nt)	renD	protein RenD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/7);  total (4/7)
*	NC_000913	568,086	0	G	A	100.0%	34.8 / NA	11	pseudogene (52/213 nt)	renD	protein RenD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/7);  total (4/7)

TCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCT  >  NC_000913/567862‑568165                                                                                                                                                                                                                                ||                                                                                ggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGgtg                                                       <  1:251491‑M2/77‑1 (MQ=255)  gcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagcgaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgacggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGTTCTGGCGAACGgtgt                                                      <  2:638803‑M2/78‑1 (MQ=255)               aaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACcggtttgc                                         >  2:227532‑M2/161‑248 (MQ=255)                aatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCt                                        >  1:729090‑M2/160‑251 (MQ=255)                      gcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCag                                   <  2:608027‑M2/97‑1 (MQ=255)                                                 atcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg       <  1:283532‑M2/125‑1 (MQ=255)                                                  tcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg       >  1:473267‑M2/126‑250 (MQ=255)                                                   caaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaa     <  1:419521‑M2/127‑1 (MQ=255)                                                   caaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGa      >  1:416149‑M2/125‑250 (MQ=255)                                                    aaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACgaac    <  2:199011‑M2/128‑2 (MQ=255)                                                      agctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCacgaacct  <  1:227532‑M2/130‑4 (MQ=255)                                                                                                                                                                                                                                ||                                                                                TCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCT  >  NC_000913/567862‑568165

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