Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	531,122	G→T	P35P (CCC→CCA)	selU ←	tRNA 2‑selenouridine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	531,122	0	G	T	100.0%	191.3 / NA	55	P35P (CCC→CCA)	selU	tRNA 2‑selenouridine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/31);  total (24/31)

TGCTGACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGC  >  NC_000913/530994‑531251                                                                                                                                 |                                                                                                                                  tgctgACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGcc                                                                                                                        <  1:858907/140‑1 (MQ=255) tgctgACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGcc                                                                                                                        <  2:981945/140‑1 (MQ=255)  gctgACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCg                                                                                                                       <  1:753020/140‑1 (MQ=255)   ctgACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGt                                                                                                                      >  2:1045356/1‑140 (MQ=255)   ctgACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGt                                                                                                                      >  1:568304/1‑140 (MQ=255)     gACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGc                                                                                                                    >  1:182528/1‑140 (MQ=255)      aCGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCt                                                                                                                   >  2:978746/1‑140 (MQ=255)                cccGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGAt                                                                                                         >  1:1078389/1‑140 (MQ=255)                 ccGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATa                                                                                                        <  2:884024/140‑1 (MQ=255)                 ccGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATa                                                                                                        <  1:929253/140‑1 (MQ=255)                         aGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGcgcg                                                                                                <  1:1045356/140‑1 (MQ=255)                         aGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGcgcg                                                                                                <  1:213954/140‑1 (MQ=255)                         aGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGcgcg                                                                                                <  2:1078389/140‑1 (MQ=255)                           tttATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGaa                                                                                              <  2:182528/140‑1 (MQ=255)                               tGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATc                                                                                          >  2:594319/1‑140 (MQ=255)                                gTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGcgcg                                                                                                >  2:919241/1‑133 (MQ=255)                                gTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGcgcg                                                                                                <  1:919241/133‑1 (MQ=255)                                 tCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACAt                                                                                           >  2:953287/1‑137 (MQ=255)                                 tCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACAt                                                                                           <  1:953287/137‑1 (MQ=255)                                        gccagcgcTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATggg                                                                                 <  1:9559/140‑1 (MQ=255)                                           agcgcTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGt                                                                              >  1:102539/1‑140 (MQ=255)                                             cgcTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTat                                                                            <  2:568304/140‑1 (MQ=255)                                             cgcTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTat                                                                            >  1:936404/1‑140 (MQ=255)                                                  cGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTatcagca                                                                       >  1:621642/1‑140 (MQ=255)                                                    tCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTatcagcaat                                                                     <  1:302454/140‑1 (MQ=255)                                                       gAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTatcagcaatcag                                                                  <  2:170262/140‑1 (MQ=255)                                                          ccTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATcagcag                                                               >  1:987848/1‑140 (MQ=255)                                                            ttGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAggg                                                             <  1:680183/140‑1 (MQ=255)                                                              gCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCa                                                           <  2:1021550/140‑1 (MQ=255)                                                                  tttATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGAt                                                       >  1:405835/1‑140 (MQ=255)                                                                    tataGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAg                                                     <  1:594319/140‑1 (MQ=255)                                                                        gCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCt                                                 <  2:987848/140‑1 (MQ=255)                                                                         cAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTg                                                >  2:347985/1‑140 (MQ=255)                                                                            gTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTTGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTc                                             >  2:862685/1‑140 (MQ=255)                                                                            gTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGtt                                              >  1:258939/1‑139 (MQ=255)                                                                            gTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGtt                                              <  2:258939/139‑1 (MQ=255)                                                                            gTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTc                                             >  2:588975/1‑140 (MQ=255)                                                                              gCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATATTGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCg                                           >  2:104324/1‑140 (MQ=255)                                                                              gCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCg                                           >  2:974399/1‑140 (MQ=255)                                                                                   cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCgtgtgt                                      <  1:978746/140‑1 (MQ=255)                                                                                   cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCgtgtgt                                      <  1:347985/140‑1 (MQ=255)                                                                                   cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCgtgtgt                                      <  1:564946/140‑1 (MQ=255)                                                                                   cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCgtgtgt                                      <  2:1069036/140‑1 (MQ=255)                                                                                   cggcggCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCgtgtgt                                      <  1:621255/140‑1 (MQ=255)                                                                                            gTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGc                             <  2:102539/140‑1 (MQ=255)                                                                                              tCATCGTTATTCATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCAt                           <  1:8095/140‑1 (MQ=255)                                                                                                         cATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATATTGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTg                <  1:104324/140‑1 (MQ=255)                                                                                                         cATTAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTg                <  2:41882/140‑1 (MQ=255)                                                                                                            tAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGc                             >  1:387701/1‑124 (MQ=255)                                                                                                            tAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGc                             <  2:387701/124‑1 (MQ=255)                                                                                                            tAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAAc             >  2:937476/1‑140 (MQ=255)                                                                                                            tAACGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAAc             >  1:453166/1‑140 (MQ=255)                                                                                                             aaCGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAAcc            <  2:936404/140‑1 (MQ=255)                                                                                                              aCGGCAGATTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCt           >  2:308748/1‑140 (MQ=255)                                                                                                                     aTTGATAGCGGCTGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCa    >  1:343462/1‑140 (MQ=255)                                                                                                                            gcggctGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGc  <  1:621874/135‑1 (MQ=255)                                                                                                                            gcggctGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGc  >  2:621874/1‑135 (MQ=255)                                                                                                                                 |                                                                                                                                  TGCTGACGAATTTCACCCGCCACCAGTTTATGTCCCAGCGCCAGCGCTGCGTCTGAGCCTTGCTGTTTATAGCAGGTGCCAACGGCGGCGCGTTCATCGTTATTCATTAACGGCAGATTGATAGCGGCGGGCATTGCGCCGTGCTCAAACTCGATAGGGGCGCGAACATCAATAATGGGCGTATCAGCAATCAGCAGGGCACGATAGTCCTGTTCCGTGTGTCTCTCTTGCATAGTTAAAAGTGAACCTCAAATCAGC  >  NC_000913/530994‑531251

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