Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	543,407	T→C	N214D (AAC→GAC)	allE ←	S‑ureidoglycine aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	543,407	0	T	C	100.0%	165.1 / NA	46	N214D (AAC→GAC)	allE	S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/23);  total (23/23)

GTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAA  >  NC_000913/543273‑543539                                                                                                                                       |                                                                                                                                     gTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATaa                                                                                                                                  >  1:981956/1‑139 (MQ=255)     cgGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAAccc                                                                                                                              >  2:369476/1‑139 (MQ=255)               cTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGa                                                                                                                    <  1:121345/139‑1 (MQ=255)               cTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGa                                                                                                                    <  2:96278/139‑1 (MQ=255)                tttCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGGTGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGaa                                                                                                                   <  2:675932/139‑1 (MQ=255)                tttCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGaa                                                                                                                   <  2:879771/139‑1 (MQ=255)                                   aaCGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGtt                                                                                                >  1:320067/1‑139 (MQ=255)                                     cGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGaa                                                                                                                   <  1:250722/118‑1 (MQ=255)                                     cGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGaa                                                                                                                   >  2:250722/1‑118 (MQ=255)                                       cTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTg                                                                                            <  1:886141/139‑1 (MQ=255)                                          caccacGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAAc                                                                                         <  2:981956/139‑1 (MQ=255)                                                gCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACgtgtgt                                                                                   <  2:90334/139‑1 (MQ=255)                                                 cccTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGtt                                                                                  >  2:640928/1‑139 (MQ=255)                                                  ccTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGttt                                                                                 <  1:345862/139‑1 (MQ=255)                                                        ccATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGatat                                                                           >  2:682774/1‑139 (MQ=255)                                                           tAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGTGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAAc                                                                        >  1:109627/1‑139 (MQ=255)                                                            aaCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGttt                                                                                 <  1:262844/129‑1 (MQ=255)                                                            aaCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGttt                                                                                 >  2:262844/1‑129 (MQ=255)                                                            aaCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAAcc                                                                       <  2:751147/139‑1 (MQ=255)                                                            aaCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAAcc                                                                       >  2:606645/1‑139 (MQ=255)                                                              ccAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGt                                                                     >  2:315351/1‑139 (MQ=255)                                                                  ccTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGtggct                                                                 >  1:221322/1‑139 (MQ=255)                                                                   cTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGtggctg                                                                >  2:171633/1‑139 (MQ=255)                                                                   cTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGtggctg                                                                <  1:862439/139‑1 (MQ=255)                                                                      tAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCa                                                             >  2:814843/1‑139 (MQ=255)                                                                      tAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCa                                                             >  1:397291/1‑139 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCAcc                                                                                                     >  1:967670/1‑89 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCAcc                                                                                                     <  2:967670/89‑1 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCaa                                                   <  2:76124/139‑1 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCaa                                                   <  2:545229/139‑1 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCaa                                                   <  1:96240/139‑1 (MQ=255)                                                                                gcgcCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCaa                                                   >  1:261910/1‑139 (MQ=255)                                                                                    ccATAAAGATGTAATCGCCTTTTTTCACCGTGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAaga                                               <  2:109627/139‑1 (MQ=255)                                                                                      aTAAAGATGTAATCTCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGTCCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGAGGCTGGCACCTGGTGCCAAagaga                                             <  1:646576/139‑1 (MQ=255)                                                                                         aaGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGa                                          <  1:682774/139‑1 (MQ=255)                                                                                              gTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGc                                     >  1:849204/1‑139 (MQ=255)                                                                                              gTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTCGTGCAAAAGAGAGGATATGc                                     >  2:874640/1‑139 (MQ=255)                                                                                                       ttttttCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGAcc                                                                                                                       >  1:935847/1‑48 (MQ=255)                                                                                                       ttttttCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGAcc                                                                                                                       <  2:935847/48‑1 (MQ=255)                                                                                                          tttCACCGGGATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCATCTGGTGCAAAAGAGAGGATATGCATGTTCATATCa                         <  1:846294/139‑1 (MQ=255)                                                                                                                  ggATCCAGTTATTGTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTaa                 >  1:57505/1‑139 (MQ=255)                                                                                                                               gTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAaga                                               <  2:835682/96‑1 (MQ=255)                                                                                                                               gTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAaga                                               >  1:835682/1‑96 (MQ=255)                                                                                                                               gTCGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGa    >  2:610793/1‑139 (MQ=255)                                                                                                                                 cGAGGTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGaaa  >  1:364395/1‑139 (MQ=255)                                                                                                                                    ggTCATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAgg                                           <  2:469451/95‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     GTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAA  >  NC_000913/543273‑543539

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