Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	532,918	A→G	E156E (GAA→GAG)	allA →	ureidoglycolate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	532,918	0	A	G	100.0%	126.1 / NA	36	E156E (GAA→GAG)	allA	ureidoglycolate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

CGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCGGCAGGCCAGGACAGGCGGAG  >  NC_000913/532783‑533052                                                                                                                                        |                                                                                                                                       cGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGctct                                                                                                                                    >  1:218940/1‑140 (MQ=255)        cAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGc                                                                                                                             <  2:340361/140‑1 (MQ=255)             aGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGAc                                                                                                                        >  2:478131/1‑140 (MQ=255)                     taccacGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCaa                                                                                                                <  2:454941/135‑1 (MQ=255)                     aCCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCaa                                                                                                                >  2:87925/1‑140 (MQ=255)                     aCCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCaa                                                                                                                <  1:294677/140‑1 (MQ=255)                          cGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                                           <  2:66561/140‑1 (MQ=255)                          cGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                                           <  1:136506/140‑1 (MQ=255)                          cGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                                           <  1:426738/140‑1 (MQ=255)                          cGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                                           <  2:290241/140‑1 (MQ=255)                          cGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                                           <  1:478131/140‑1 (MQ=255)                           gTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACt                                                                                                          >  2:130532/1‑140 (MQ=255)                               cGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCa                                                                                                      >  2:291245/1‑140 (MQ=255)                                  cTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATaa                                                                                                   <  1:461042/140‑1 (MQ=255)                                    ggCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGa                                                                                                 >  1:391632/1‑140 (MQ=255)                                     gCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGAt                                                                                                >  1:213885/1‑140 (MQ=255)                                         cACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAAc                                                                                            <  2:281263/140‑1 (MQ=255)                                          aCCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCa                                                                                                      <  1:134989/129‑1 (MQ=255)                                          aCCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCa                                                                                                      >  2:134989/1‑129 (MQ=255)                                           cccACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTa                                                                                          >  1:505523/1‑140 (MQ=255)                                                          gCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTa                                                                           >  1:475670/1‑140 (MQ=255)                                                            aGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATg                                                                         <  2:342177/140‑1 (MQ=255)                                                                         ttttCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGt                                                            <  2:475670/140‑1 (MQ=255)                                                                               gACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGTAGt                                                      <  2:218940/140‑1 (MQ=255)                                                                               gACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGt                                                      <  2:505523/140‑1 (MQ=255)                                                                               gACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCCGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGt                                                      >  1:43341/1‑140 (MQ=255)                                                                                aCCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGtt                                                     >  1:343279/1‑140 (MQ=255)                                                                                              cggcAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGAc                                       <  1:130532/140‑1 (MQ=255)                                                                                                                tGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCgg                     <  1:225643/140‑1 (MQ=255)                                                                                                                tGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCgg                     <  1:124908/140‑1 (MQ=255)                                                                                                                tGTTGAAAGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCgg                     <  2:89366/140‑1 (MQ=255)                                                                                                                     aaaGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACgg                          >  2:459181/1‑130 (MQ=255)                                                                                                                     aaaGTATTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACgg                          <  1:459181/130‑1 (MQ=255)                                                                                                                          aTTCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCGGCAGGCCAGGa           <  1:311453/140‑1 (MQ=255)                                                                                                                            tCCTGAACAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCGGCAGGCCAGGACa         >  1:256364/1‑140 (MQ=255)                                                                                                                                   cAGGAGCTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCGGCAGGCCAGGACAGGCGGAg  >  2:173253/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       CGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGACGGCGCGGCAGGCCAGGACAGGCGGAG  >  NC_000913/532783‑533052

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