Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,923,308	T→C	G159G (GGA→GGG)	truC ←	tRNA(Ile1,Asp) pseudouridine(65) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,308	0	T	C	100.0%	183.3 / NA	54	G159G (GGA→GGG)	truC	tRNA(Ile1,Asp) pseudouridine(65) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/23);  total (31/23)

TCAATGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGCGAT  >  NC_000913/2923160‑2923454                                                                                                                                                     |                                                                                                                                                   tCAATGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGccc                                                                                                                                                  >  1:145655/1‑151 (MQ=255)   aaTGACAACTGACTGGTATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGTGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCg                                                                                                                                                 <  1:547664/150‑1 (MQ=255)   aaTGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCgg                                                                                                                                                <  2:902422/151‑1 (MQ=255)    aTGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGt                                                                                                                                               <  1:618184/151‑1 (MQ=255)        cAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTc                                                                                                                                           >  1:885863/1‑151 (MQ=255)         aaCTGACTGGCATGTAGCATTAACCGCTGGAGGGCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGAAGCAAATGGTCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTc                                                                                                                                           >  2:515082/1‑150 (MQ=255)             gactgGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGc                                                                                                                                        <  1:111719/149‑1 (MQ=255)                    aTGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTcc                                                                                                                               >  2:825774/1‑151 (MQ=255)                         gCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAg                                                                                                                          <  1:176210/151‑1 (MQ=255)                         gCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATATTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTccacc                                                                                                                            >  2:683538/1‑149 (MQ=255)                            ttAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGcc                                                                                                                       >  2:783441/1‑151 (MQ=255)                                 cGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAgc                                                                                                                  <  1:748240/151‑1 (MQ=255)                                    tGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGcgcg                                                                                                               >  2:711870/1‑151 (MQ=255)                                     ggAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGt                                                                                                              >  2:175089/1‑151 (MQ=255)                                         gCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTc                                                                                                          <  1:731082/151‑1 (MQ=255)                                                tGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAAc                                                                                                   >  2:228867/1‑151 (MQ=255)                                                 gCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACg                                                                                                  <  2:73287/151‑1 (MQ=255)                                                  cTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGt                                                                                                 >  1:301270/1‑151 (MQ=255)                                                    cAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTc                                                                                                <  1:228794/150‑1 (MQ=255)                                                               cGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCg                                                                                    >  1:179771/1‑151 (MQ=255)                                                                 aTTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGc                                                                                  >  2:705625/1‑151 (MQ=255)                                                                     tGGCGTAAATCGCTATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCAttt                                                                              >  2:681768/1‑151 (MQ=255)                                                                        cGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCt                                                                            >  1:384055/1‑150 (MQ=255)                                                                           aaaTCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACgg                                                                        <  1:775587/151‑1 (MQ=255)                                                                                  cATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGcc                                                                   >  1:841034/1‑149 (MQ=255)                                                                                  cATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGcc                                                                   >  1:259926/1‑149 (MQ=255)                                                                                     gTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGAcc                                                              <  2:432218/151‑1 (MQ=255)                                                                                     gTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGAcc                                                              <  2:820189/151‑1 (MQ=255)                                                                                            aTCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATa                                                        >  2:476068/1‑150 (MQ=255)                                                                                            aTCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATa                                                        <  2:148409/150‑1 (MQ=255)                                                                                              cGCCAATAATCGGATGACGCAGATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAAt                                                      <  1:19274/150‑1 (MQ=255)                                                                                              cGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATg                                                     >  1:685185/1‑151 (MQ=255)                                                                                                ccAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCg                                                   >  2:413167/1‑151 (MQ=255)                                                                                                  aataatCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTc                                                 >  2:461223/1‑151 (MQ=255)                                                                                                      atCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCActg                                             >  1:407479/1‑151 (MQ=255)                                                                                                         ggATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCActgctg                                          >  1:132873/1‑151 (MQ=255)                                                                                                          gATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTgctgg                                         <  1:561777/151‑1 (MQ=255)                                                                                                          gATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTgctgg                                         >  2:865694/1‑151 (MQ=255)                                                                                                               cGCAAATGGGCAAGATGTCGGCGAAGCTGATGATTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCg                                    >  1:620569/1‑151 (MQ=255)                                                                                                               agCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCg                                    >  2:430605/2‑151 (MQ=255)                                                                                                                       ggCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTa                            <  1:290501/151‑1 (MQ=255)                                                                                                                          aaGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCt                         >  1:603843/1‑151 (MQ=255)                                                                                                                            gATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTc                       <  2:87382/151‑1 (MQ=255)                                                                                                                               gTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCGCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCCTTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGTGCCTTTATCTTcg                      >  2:396010/1‑149 (MQ=255)                                                                                                                                 cGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCggg                  <  2:135822/151‑1 (MQ=255)                                                                                                                                  ggCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGc                 <  1:23073/151‑1 (MQ=255)                                                                                                                                     gAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCaa               >  2:296884/1‑150 (MQ=255)                                                                                                                                      aaGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAt             <  1:82721/151‑1 (MQ=255)                                                                                                                                      aaGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAt             >  1:361789/1‑151 (MQ=255)                                                                                                                                       aGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCTCCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAt             <  2:244208/150‑1 (MQ=255)                                                                                                                                        gCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAttt           <  2:581969/151‑1 (MQ=255)                                                                                                                                           gATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTAt         <  1:11469/150‑1 (MQ=255)                                                                                                                                            aTGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATc        >  2:324731/1‑150 (MQ=255)                                                                                                                                              gTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGc     >  2:28801/1‑151 (MQ=255)                                                                                                                                                 tGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGCGAt  >  1:530457/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TCAATGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGCGAT  >  NC_000913/2923160‑2923454

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