Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,258,212	T→C	E52E (GAA→GAG)	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,212	0	T	C	100.0%	126.7 / NA	36	E52E (GAA→GAG)	psuG	pseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/18);  total (18/18)

AACTTTGGTCACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTT  >  NC_000913/2258086‑2258355                                                                                                                               |                                                                                                                                                aaCTTTGGTCACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTggg                                                                                                                           <  1:276031/149‑1 (MQ=255) aaCTTTGGTCACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGc                                                                                                                          <  1:153198/150‑1 (MQ=255)    tttGGTCACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGGCTGGGCAttt                                                                                                                      >  2:372370/1‑151 (MQ=255)          cACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTggg                                                                                                                <  2:379426/151‑1 (MQ=255)            cGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGa                                                                                                               <  1:577552/150‑1 (MQ=255)                          cGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAACCCCCTGTTTACGAATAGGTTCCTCAACTTCAATTGCGGTCTGGGCCTTTTTTGGGGACGGCATccc                                                                                                     >  2:2421/1‑146 (MQ=255)                                       cAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCg                                                                                    <  2:65640/150‑1 (MQ=255)                                         aTTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCCTTTTGTGGGAACGGCATCCCGTGGGAAATAATGGTCGAtt                                                                                 >  2:549802/1‑151 (MQ=255)                                            tcttctTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGGTTACAAAGAGATTCCTCAACTTCACTTGCGGTCTGGGCATTTTGTGGGAACGGGATCCCGTGAGAAATAATGGGCGATTcc                                                                               >  1:269450/1‑150 (MQ=255)                                            tcttctTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAATGCCCCGTTTACGAATAGGTTTCTCAACTTCCATTGCGGTCTGGGCATTTTTTGGGAAAGGCATCCCGTGAGAAATAATGGTCGCtt                                                                                 >  2:390647/1‑148 (MQ=255)                                                  ttGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAATTTCCTCAACTTCAATTGGGGTCTGGGGATTTTGGGGGAACGGCATCCCGTGAGAAATAATCGGCGATTCCAGCGcc                                                                         >  1:336727/1‑150 (MQ=255)                                                    gCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCcaca                                                                      <  2:215968/151‑1 (MQ=255)                                                          ccccACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCCGTTTACGAATAGTTTCCTCAAATTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCgg                                                                 >  1:491190/1‑150 (MQ=255)                                                           cccACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGGCTGGGCCTTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCgg                                                                 >  2:409837/1‑149 (MQ=255)                                                            ccACTTTCATCACGCCGCCGATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAAGTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgtt                                                               <  2:97104/150‑1 (MQ=255)                                                               cTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCCCAACCGggttttt                                                           >  1:140209/1‑151 (MQ=255)                                                                tttCATCACGCCGCCAATAATGGCGATCGTGGCAGGGACAGCGCCCTGTTTACGAATAGTTTTCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgtttttt                                                           >  1:662919/1‑150 (MQ=255)                                                                 ttCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttg                                                          >  1:97129/1‑150 (MQ=255)                                                                 ttCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttg                                                          <  2:164371/150‑1 (MQ=255)                                                                      caCGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgttttt                                                     <  2:3737/150‑1 (MQ=255)                                                                        cgccgcCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTaa                                                  <  1:5576/151‑1 (MQ=255)                                                                         gccgccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGgttttttgtttttt                                                    >  2:560222/1‑148 (MQ=255)                                                                             ccAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGt                                             <  1:19045/151‑1 (MQ=255)                                                                                       cGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTc                                   <  2:618661/151‑1 (MQ=255)                                                                                             tGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCgggg                              >  2:116813/1‑150 (MQ=255)                                                                                               gCAGGTACAGCGCCCTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGaaa                           <  2:615758/151‑1 (MQ=255)                                                                                                           cccTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAAtt               <  2:354584/151‑1 (MQ=255)                                                                                                            ccTGTTTACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTc              <  1:559116/151‑1 (MQ=255)                                                                                                                tttACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAggg          >  1:267999/1‑151 (MQ=255)                                                                                                                tttACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAggg          >  1:601541/1‑151 (MQ=255)                                                                                                                 ttACGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAgggg         >  2:558983/1‑151 (MQ=255)                                                                                                                   aCGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGaa       <  1:409848/151‑1 (MQ=255)                                                                                                                   aCGAATAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGa        >  1:595671/1‑150 (MQ=255)                                                                                                                      aaTAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAtt    <  1:525539/151‑1 (MQ=255)                                                                                                                        tAGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAtttt  >  2:566294/1‑151 (MQ=255)                                                                                                                         aGTTTCCTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAAtttt  <  2:618122/150‑1 (MQ=255)                                                                                                                               |                                                                                                                                                AACTTTGGTCACGTTATGCCCTTCACGACCCAGTAATTCAATTTCTTCTTTGCTTAACCCCACTTTCATCACGCCGCCAATAATGGCGATCGTGGCAGGTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTT  >  NC_000913/2258086‑2258355

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