Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,258,247	A→G	F41L (TTC→CTC)	psuG ←	pseudouridine 5'‑phosphate glycosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,258,247	0	A	G	100.0%	172.4 / NA	50	F41L (TTC→CTC)	psuG	pseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/31);  total (19/31)

ATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATAT  >  NC_000913/2258101‑2258390                                                                                                                                                   |                                                                                                                                                aTGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCgg                                                                                                                                              <  1:279721/150‑1 (MQ=255)          aCGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTga                                                                                                                                    <  2:191144/151‑1 (MQ=255)            gACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTTCGAATAGTTTTTTCAACTTCTAATGCGGTCTGGGCATTTTGTGGGAGCGGGATTCCg                                                                                                                                       >  1:604993/1‑146 (MQ=255)                cAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAaat                                                                                                                               >  1:215737/1‑150 (MQ=255)                  gTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAaataat                                                                                                                            <  2:33493/151‑1 (MQ=255)                    aaTTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATgg                                                                                                                          <  2:592969/151‑1 (MQ=255)                       tCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTc                                                                                                                        <  1:132424/150‑1 (MQ=255)                            ttcttcTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTc                                                                                                                   <  2:440545/150‑1 (MQ=255)                              cttcttTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGGCTGGGCATTTTGTGGGAGCGGCATCCCGCGAGAAATAATGGTCGATTcc                                                                                                                  >  1:259602/1‑149 (MQ=255)                                tcttTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCa                                                                                                                 <  1:211114/148‑1 (MQ=255)                                    tGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCca                                                                                                           <  1:129691/150‑1 (MQ=255)                                      cTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAACGCCCTGTTTACGAATAGGTTCTTCAACTTCAATTGCGGGCTGGGCATTTTGTGGGGGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCcaca                                                                                                         >  1:305757/1‑150 (MQ=255)                                       ttAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACaa                                                                                                        >  2:289846/1‑150 (MQ=255)                                        tAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAAc                                                                                                       <  1:337463/150‑1 (MQ=255)                                                 tttCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgtttttt                                                                                              >  2:452061/1‑150 (MQ=255)                                                  ttCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttg                                                                                             <  1:542727/150‑1 (MQ=255)                                                   tcatcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgt                                                                                            <  1:378433/150‑1 (MQ=255)                                                   tcatcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgt                                                                                            <  1:433696/150‑1 (MQ=255)                                                    catcaCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCGgttttttgttt                                                                                          >  2:135178/1‑151 (MQ=255)                                                         cgccgcCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTTTTCAACTTCAATTGCGGGCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCGgttttttgtttttt                                                                                       >  1:31746/1‑149 (MQ=255)                                                          gccgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTaaa                                                                                    <  2:31745/151‑1 (MQ=255)                                                           ccgccAATAATGGCGATCGTGGCAGGTACAGCGCCCCGTTTACGAATAGTTTTTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgtttttt                                                                                       >  1:632833/1‑147 (MQ=255)                                                            cgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTaaa                                                                                    <  1:432587/149‑1 (MQ=255)                                                            cgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGc                                                                                  >  2:503773/1‑151 (MQ=255)                                                                 ataatGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCt                                                                             <  2:613452/151‑1 (MQ=255)                                                                    atGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCa                                                                          <  1:539120/151‑1 (MQ=255)                                                                     tGGCGATCGTGGCAGGTACAGCGCCCTGTTTTCGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCGGGTTTTTGTTTTTTAAAGCGTCCTGCAc                                                                         >  1:388728/1‑151 (MQ=255)                                                                          aTCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCg                                                                    <  1:527138/151‑1 (MQ=255)                                                                            cGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgg                                                                   <  1:413091/150‑1 (MQ=255)                                                                                           gcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTaataat                                                   <  2:51078/151‑1 (MQ=255)                                                                                            cccTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAAtt                                                  <  1:327246/151‑1 (MQ=255)                                                                                                 tttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAggg                                             >  2:652836/1‑151 (MQ=255)                                                                                                  ttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCGCTTCCGGGGAAATTTGTAATAATTCAggg                                             <  1:666380/150‑1 (MQ=255)                                                                                                   tACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAgggg                                            >  1:460101/1‑150 (MQ=255)                                                                                                       aaTAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAtt                                       >  1:310904/1‑151 (MQ=255)                                                                                                        aTAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAttt                                      >  1:202828/1‑151 (MQ=255)                                                                                                           gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTa                                   <  2:202829/151‑1 (MQ=255)                                                                                                              tcttcAACTTCAATTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAAt                                 <  2:135323/150‑1 (MQ=255)                                                                                                                ttcaacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTc                               >  2:114021/1‑150 (MQ=255)                                                                                                                  caacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCTGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAg                             <  1:51626/150‑1 (MQ=255)                                                                                                                   aacttcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGAc                           >  2:29313/1‑151 (MQ=255)                                                                                                                       tcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACAtt                        <  1:289846/150‑1 (MQ=255)                                                                                                                       tcaaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACAt                         <  2:496534/149‑1 (MQ=255)                                                                                                                        caaTTGCGGTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCt                      <  2:49349/151‑1 (MQ=255)                                                                                                                                gTCTGGGCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTtctc              <  1:29314/151‑1 (MQ=255)                                                                                                                                    gggCATTTTGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCAc           <  2:477942/150‑1 (MQ=255)                                                                                                                                          tttGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGat    >  1:486921/1‑151 (MQ=255)                                                                                                                                           ttGTGGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGat    >  1:205471/1‑150 (MQ=255)                                                                                                                                            tgtgGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGata   <  1:212662/150‑1 (MQ=255)                                                                                                                                             gtgGGAGCGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCTCTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGatat  <  2:388732/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                ATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATAT  >  NC_000913/2258101‑2258390

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