Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	566,728	G→T	pseudogene (45/93 nt)	renD →	protein RenD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,728	0	G	T	100.0%	58.9 / NA	18	pseudogene (45/93 nt)	renD	protein RenD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (8/10);  total (8/10)

AATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCT  >  NC_000913/566484‑566892                                                                                                                                                                                                                                                      |                                                                                                                                                                     aaTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTc                                                                                                                                                                >  1:29763/1‑251 (MQ=255)                                aaggccgcagatgagcttgtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTACGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                                                                                                                  >  1:70041/24‑251 (MQ=255)                                aaggccgcagatgagcttgtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                                                                                                                  <  1:107533/228‑1 (MQ=255)                                 aggccgcagatgagcttgtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGGGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgcg                                                                                                                                 <  2:70041/229‑1 (MQ=255)                                        agatgagcttgtccatATGACTGCGAGAATTACCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTACCAGGCGCAAcagt                                                                                                                          >  2:140430/16‑248 (MQ=255)                                              gcttgtccatATGACTGCGAGAATTAACCGTGGCGAGGCGATCCCTGAACCCGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACCGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGCCTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGCCGCATAATAGCTTCt                                                                                                                                                               >  1:277025/10‑208 (MQ=255)                                                              gcgagaattAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccg                                                                                                      <  1:340592‑M1/245‑19 (MQ=255)                                                              gcgagaattAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACGTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcg                                                                                                    <  2:128643‑M1/247‑21 (MQ=255)                                                                  gAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGAGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACTTTGCCGCGCTAACAGGCGCCACAGTAAccagcataaatcaggccgcggct                                                                                                 >  2:222820‑M1/1‑227 (MQ=255)                                                                                         ccTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctg                                                                         <  2:107889‑M1/251‑48 (MQ=255)                                                                                            gAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggtt                                                                      >  1:315790‑M1/1‑201 (MQ=255)                                                                                                               cTGTCATGGGCGGTGGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggc                                                   >  1:236851‑M1/1‑182 (MQ=255)                                                                                                                      gggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacgg                                             <  2:236851‑M1/250‑76 (MQ=255)                                                                                                                           gTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctgggtatcgaaggtaaggtctggcgaacggtgtatt                                       <  2:256688‑M1/251‑82 (MQ=255)                                                                                                                              gACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattacc                                    <  2:29763‑M1/251‑85 (MQ=255)                                                                                                                                          gTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctacca                        <  2:95780‑M1/251‑97 (MQ=255)                                                                                                                                                aGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccaggga                    >  2:11150‑M1/1‑149 (MQ=255)                                                                                                                                                                cGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagat  <  1:235474‑M1/251‑119 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                     AATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCT  >  NC_000913/566484‑566892

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