Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	538,138	A→G	N169S (AAC→AGC)	ybbW →	putative allantoin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	538,138	0	A	G	100.0%	157.1 / NA	46	N169S (AAC→AGC)	ybbW	putative allantoin transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/25);  total (21/25)

CGCGGGGTCACTGGCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAACGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTTTGAC  >  NC_000913/538007‑538271                                                                                                                                    |                                                                                                                                      cgcgGGGTCACTGGCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGta                                                                                                                               >  1:64741/1‑140 (MQ=255)      ggTCACTGGCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGt                                                                                                                          >  1:187168/1‑140 (MQ=255)       ggtcgctgCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGtt                                                                                                                         >  2:190622/8‑140 (MQ=255)            tGGCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGc                                                                                                                    <  2:513523/140‑1 (MQ=255)                  gCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTGCTGGCTGATCAGCGGTGGTATGGGTTTTGGCGGTGGc                                                                                                              >  1:241380/1‑140 (MQ=255)                   cTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCa                                                                                                             <  1:190622/140‑1 (MQ=255)                   cTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTCTCTTAATCCTCTGGCTTGTCAGCGTTGGTAGAGGTTTTGGCGGTGGCa                                                                                                             >  2:262572/1‑140 (MQ=255)                         tCTGATTGGCAAAATCTGGCCGGGATTTTTAACGCACGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGttt                                                                                                       <  2:353707/140‑1 (MQ=255)                                  cAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTAGGTATAGCTTTTGGCGGTGGCTAAGTTTTaaataaat                                                                                              >  1:283493/1‑140 (MQ=255)                                         tGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGc                                                                                       <  2:40483/140‑1 (MQ=255)                                         tGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGc                                                                                       <  2:6971/140‑1 (MQ=255)                                         tGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCACTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGc                                                                                       <  1:414881/140‑1 (MQ=255)                                                 aTTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTa                                                                               <  1:633771/140‑1 (MQ=255)                                                   ttttAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAAc                                                                             <  1:170001/140‑1 (MQ=255)                                                      tAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCg                                                                          <  1:429313/140‑1 (MQ=255)                                                           ctcGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCAt                                                                     <  2:720429/140‑1 (MQ=255)                                                              ggtggtGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCta                                                                  >  1:76374/1‑140 (MQ=255)                                                              ggtggtGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCta                                                                  >  1:113979/1‑140 (MQ=255)                                                                          aCTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTcgg                                                      >  2:532396/1‑140 (MQ=255)                                                                          aCTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTAGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTcgg                                                      <  2:758467/140‑1 (MQ=255)                                                                                ttAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTAt                                                >  2:710115/1‑140 (MQ=255)                                                                                        ttCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAAACCGTGCAGCTATATGGTTTTCCGCGGTATGGCGAttt                                        >  1:231862/1‑140 (MQ=255)                                                                                          ttctACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCTGTGCATCAATATCGATTTCGGTGGTATGGCGATTTgg                                      >  1:227955/2‑140 (MQ=255)                                                                                              aCCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGa                                  <  1:710115/140‑1 (MQ=255)                                                                                               ccGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGAt                                 >  2:55584/1‑140 (MQ=255)                                                                                                  ggCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTc                              >  1:554466/1‑140 (MQ=255)                                                                                                   gCTTAATTACTTTCTTAATCTTCTGGCTGTTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  1:280390/140‑1 (MQ=255)                                                                                                   gCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  2:90587/140‑1 (MQ=255)                                                                                                   gCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  1:532396/140‑1 (MQ=255)                                                                                                   gCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  1:592252/140‑1 (MQ=255)                                                                                                   gCTTAATTACTTTCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  2:187168/140‑1 (MQ=255)                                                                                                   gCTTAATTACTTTCTCAATCTTCTGCCTGTTCAGCGTAGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTTACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCa                             <  1:262572/140‑1 (MQ=255)                                                                                                             tttCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGa                   >  2:285809/1‑140 (MQ=255)                                                                                                             tttCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGa                   <  1:55584/140‑1 (MQ=255)                                                                                                             tttCTTAATCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGa                   <  2:64741/140‑1 (MQ=255)                                                                                                                   aaTCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTc             <  2:554466/140‑1 (MQ=255)                                                                                                                   aaTCTTCTGGCTGGTCAGCGTTGGTATAGGTTTTGCCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTc             <  2:283493/140‑1 (MQ=255)                                                                                                                         ctggctggTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCt       >  1:598028/1‑140 (MQ=255)                                                                                                                         ctggctggTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCt       >  2:731611/1‑140 (MQ=255)                                                                                                                         ctggctggTCAGCGTTGGTATAGGGTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGGCCAATCt       >  2:569416/1‑140 (MQ=255)                                                                                                                          tggctggTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCtt      >  2:414887/1‑140 (MQ=255)                                                                                                                              tggTCAGCGTTGGTATATGTTTCGGCGGTGGAAAATTTTCAAATAAATCCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTTTGAc  <  2:231862/140‑1 (MQ=255)                                                                                                                               ggTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTTTGa   <  2:204521/138‑1 (MQ=255)                                                                                                                               ggTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTTTGa   >  1:204521/1‑138 (MQ=255)                                                                                                                                gTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTcc            <  2:509805/128‑1 (MQ=255)                                                                                                                                gTCAGCGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTcc            >  1:509805/1‑128 (MQ=255)                                                                                                                                    |                                                                                                                                      CGCGGGGTCACTGGCCTGCTTGATTCTGATTGGCAAAATCTGGCCGGGATTTTTAACTCTCGGTGGTGATTTCACTCTGTTAGGCCTTTCTCTACCGGGCTTAATTACTTTCTTAATCTTCTGGCTGGTCAACGTTGGTATAGGTTTTGGCGGTGGCAAAGTTTTAAATAAATTCACTGCCATTCTTAACCCGTGCATCTATATCGTTTTCGGCGGTATGGCGATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTTTGAC  >  NC_000913/538007‑538271

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