Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,923,326	T→C	E153E (GAA→GAG)	truC ←	tRNA pseudouridine(65) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,326	0	T	C	100.0%	142.4 / NA	40	E153E (GAA→GAG)	truC	tRNA pseudouridine(65) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/23);  total (17/23)

CCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGC  >  NC_000913/2923191‑2923451                                                                                                                                        |                                                                                                                              ccGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGgctccagctcca                                                                                                                           <  1:864762/140‑1 (MQ=255)   gCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTccacc                                                                                                                         <  2:876852/140‑1 (MQ=255)    cTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTccacca                                                                                                                        <  2:112836/140‑1 (MQ=255)      ggAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAgg                                                                                                                      <  2:202887/140‑1 (MQ=255)       gAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGc                                                                                                                     <  2:75036/140‑1 (MQ=255)          gCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTccacca                                                                                                                        >  2:98476/1‑134 (MQ=255)          gCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTccacca                                                                                                                        <  1:98476/134‑1 (MQ=255)          gCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGt                                                                                                                  <  2:647511/140‑1 (MQ=255)                aTGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGcgcg                                                                                                            <  1:209702/140‑1 (MQ=255)                aTGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGcgcg                                                                                                            <  2:500432/140‑1 (MQ=255)                    tCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGt                                                                                                        >  2:234447/1‑140 (MQ=255)                    tCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGt                                                                                                        <  1:85326/140‑1 (MQ=255)                        cACCGCTGCGATTCTGGCGTAAATCACCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCggg                                                                                                    >  2:572257/1‑140 (MQ=255)                         aCCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTc                                                                                                       <  2:766279/136‑1 (MQ=255)                         aCCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTc                                                                                                       >  1:766279/1‑136 (MQ=255)                                cGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTcc                                                                                            >  2:48903/1‑140 (MQ=255)                                       ggCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCt                                                                                     <  1:671588/140‑1 (MQ=255)                                             aaTCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGc                                                                               >  1:548400/1‑140 (MQ=255)                                             aaTCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGc                                                                               <  1:48903/140‑1 (MQ=255)                                                  ccATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTc                                                                          >  1:3921/1‑140 (MQ=255)                                                     tGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTAc                                                                       >  2:696711/1‑140 (MQ=255)                                                       tttGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACgg                                                                     >  1:658821/1‑140 (MQ=255)                                                               cGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGa                                                             >  1:84423/1‑140 (MQ=255)                                                                 ccAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGAcc                                                           >  1:622270/1‑140 (MQ=255)                                                                 ccAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGAcc                                                           <  1:234447/140‑1 (MQ=255)                                                                     taatCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGa                                                       <  2:874675/140‑1 (MQ=255)                                                                      aatCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGAt                                                      <  2:460273/140‑1 (MQ=255)                                                                          ggATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATg                                                  <  2:3921/140‑1 (MQ=255)                                                                            aTGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCg                                                <  2:548400/140‑1 (MQ=255)                                                                            aTGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCg                                                >  2:877420/1‑140 (MQ=255)                                                                            aTGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCg                                                >  1:608734/1‑140 (MQ=255)                                                                                          cAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGc                                  <  2:608734/140‑1 (MQ=255)                                                                                          cAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGc                                  <  1:877420/140‑1 (MQ=255)                                                                                            aGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCgg                                <  2:622270/140‑1 (MQ=255)                                                                                            aGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCgg                                <  1:696711/140‑1 (MQ=255)                                                                                                 tCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCtt                           >  2:806228/1‑140 (MQ=255)                                                                                                  cGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCttt                          >  1:311836/1‑140 (MQ=255)                                                                                                                    cGTCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAttt        >  2:355052/1‑140 (MQ=255)                                                                                                                      tCCGGTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTAt      <  2:311836/140‑1 (MQ=255)                                                                                                                         ggTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAg   <  1:355052/140‑1 (MQ=255)                                                                                                                          gTTTTCGGCTCCAGCTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGc  >  2:17228/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                              CCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGC  >  NC_000913/2923191‑2923451

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