Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,256,613	(T)5→4	coding (723/1251 nt)	psuT ←	putative pseudouridine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,256,609	0	T	.	100.0%	136.1 / NA	33	coding (727/1251 nt)	psuT	putative pseudouridine transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (19/14);  total (19/14)

AAACGCCATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCTTTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAGAATAC  >  NC_000913/2256529‑2256683                                                                                 |                                                                           aaaCGCCATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCt                                                                            >  2:399412/1‑81 (MQ=255)   aCGCTATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttt                                                                          <  2:399376/81‑1 (MQ=255)     gCCATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttttg                                                                        <  2:399377/81‑2 (MQ=255)     gCCATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttttg                                                                        <  2:399378/81‑2 (MQ=255)      ccATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttttgg                                                                       >  1:399437/1‑79 (MQ=255)       cATTAACACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttttggc                                                                      >  1:399413/1‑78 (MQ=255)        aTTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCttttggcg                                                                     >  1:399457/1‑77 (MQ=255)          tACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTggcggc                                                                  >  1:399415/1‑81 (MQ=255)          tACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTggcggc                                                                  >  1:399423/1‑81 (MQ=255)           accaccGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTggcggc                                                                  <  1:399379/80‑1 (MQ=255)             caccGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGtt                                                               >  1:399399/1‑81 (MQ=255)              accGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGttt                                                              <  2:399380/81‑1 (MQ=255)               ccGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTc                                                             <  2:399381/81‑1 (MQ=255)                 gTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGc                                                           >  2:399430/1‑81 (MQ=255)                  tCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCt                                                          <  1:399382/81‑1 (MQ=255)                    gCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTg                                                         <  1:399383/80‑1 (MQ=255)                          ccAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGa                                                    >  1:399427/1‑79 (MQ=255)                          ccAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACa                                                  <  2:399384/81‑1 (MQ=255)                            aGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGa                                                    >  2:399447/1‑77 (MQ=255)                             gcggcgATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGAt                                               <  1:399385/81‑1 (MQ=255)                                gcgATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGAtttt                                            >  1:399394/1‑81 (MQ=255)                                   aTTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCaa                                         <  1:399386/81‑1 (MQ=255)                                        tAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAAtgt                                     >  1:399421/1‑80 (MQ=255)                                        tAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAAtgt                                     >  1:399420/1‑80 (MQ=255)                                        tAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAAtgt                                     <  2:399387/80‑1 (MQ=255)                                            ccGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCt                                <  1:399388/81‑1 (MQ=255)                                                  aTCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAg                          >  2:399424/1‑81 (MQ=255)                                                  aTCGCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAg                          >  2:399406/1‑81 (MQ=255)                                                     gCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCt                       >  1:399448/1‑81 (MQ=255)                                                     gCACCGCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCt                       >  1:399425/1‑81 (MQ=255)                                                          gCTCGCCGCCGCTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGtt                  <  2:399389/81‑1 (MQ=255)                                                             cgccgccgcTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGcc               >  1:399443/1‑81 (MQ=255)                                                              gccgccgcTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCg              >  1:399456/1‑81 (MQ=255)                                                               ccgccgcTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCgg             <  2:399390/81‑1 (MQ=255)                                                                    gcTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCt          >  1:399433/1‑79 (MQ=255)                                                                     cTTCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAg       <  2:399391/81‑1 (MQ=255)                                                                      ttCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAg       >  1:399441/1‑80 (MQ=255)                                                                       tCGATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAGaa     >  1:399431/1‑81 (MQ=255)                                                                         gATAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAGAAt    <  1:399392/80‑1 (MQ=255)                                                                          aTAAAGC‑TTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAGAATAc  >  1:399435/1‑81 (MQ=255)                                                                                 |                                                                           AAACGCCATTACCACCGTCGCCACACCAGCGGCGATTTTTAGCCCGGTCATCGCACCGCTCGCCGCCGCTTCGATAAAGCTTTTTGGCGGCGTTTCGCTGAACGACAGATTTTCAAATGTGACCTGCGAAGGCTCGGTTGCCGGGCTAAGAATAC  >  NC_000913/2256529‑2256683

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