Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,052,156	C→T	intergenic (+142/+120)	yeeJ → / ← yeeL	putative adhesin/pseudogene, glycosyltransferase homology

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,052,156	0	C	T	100.0%	110.5 / NA	32	intergenic (+142/+120)	yeeJ/yeeL	putative adhesin/pseudogene, glycosyltransferase homology
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/19);  total (13/19)

GGCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATCATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTC  >  NC_000913/2052052‑2052275                                                                                                         |                                                                                                                        ggCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTg                                                                                       <  2:490307/139‑1 (MQ=255)                  tCCCTCTACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTg                                                                      >  1:164567/1‑139 (MQ=255)                       ctACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAg                                                                 <  2:42416/139‑1 (MQ=255)                       ctACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAg                                                                 <  2:422268/139‑1 (MQ=255)                       ctACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAg                                                                 >  1:498629/1‑139 (MQ=255)                             aaCTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAGAGGCTGGTTTTCTGTGATGATTGAc                                                                    <  2:435151/130‑1 (MQ=255)                             aaCTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGAc                                                                    >  1:435151/1‑130 (MQ=255)                               cTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCagag                                                         >  1:122485/1‑139 (MQ=255)                                  gCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTAt                                                      >  2:331743/1‑139 (MQ=255)                                  gCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTAt                                                      >  1:446253/1‑139 (MQ=255)                                             acaGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGAcacca                                           <  2:122485/139‑1 (MQ=255)                                              caGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCaga                                                          <  1:273370/123‑1 (MQ=255)                                              caGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCaga                                                          >  2:273370/1‑123 (MQ=255)                                                 aaaaTGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGt                                       >  2:100339/1‑139 (MQ=255)                                                       aataaataAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATaa                                 <  1:331743/139‑1 (MQ=255)                                                          aaataaaAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAg                              <  2:34971/139‑1 (MQ=255)                                                           aataaaAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGAcacca                                           <  2:12771/125‑1 (MQ=255)                                                           aataaaAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGAcacca                                           >  1:12771/1‑125 (MQ=255)                                                            ataaaAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGAt                            <  1:100339/139‑1 (MQ=255)                                                                  aTGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCAt                      <  2:88824/139‑1 (MQ=255)                                                                    gCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTa                    <  1:495380/139‑1 (MQ=255)                                                                       gCACCGCCAGAATCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGACTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTACTGATGACACCACCGTTCATAACATATCGGCATTACTa                 <  1:343608/139‑1 (MQ=255)                                                                                 aaTCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTa       <  1:428991/139‑1 (MQ=255)                                                                                  aTCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCa                       >  2:22945/1‑122 (MQ=255)                                                                                  aTCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCa                       <  1:22945/122‑1 (MQ=255)                                                                                   tCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACt     <  2:143395/139‑1 (MQ=255)                                                                                   tCGCGTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACt     <  2:498629/139‑1 (MQ=255)                                                                                      cgTTCGATGTGACTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGtc  >  1:68847/1‑139 (MQ=255)                                                                                                 aCTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAg                              >  1:128418/1‑100 (MQ=255)                                                                                                 aCTATGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAg                              <  2:128418/100‑1 (MQ=255)                                                                                                     tGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTa                 >  1:64084/1‑109 (MQ=255)                                                                                                     tGATTATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTa                 <  2:64084/109‑1 (MQ=255)                                                                                                         |                                                                                                                        GGCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTTCTCTCACAGAAAATGAATAAATAAAAATGCGGCACCGCCAGAATCGCGTTCGATGTGACTATGATCATTTTCGGTGAGTACTGCTTTGATCTAAAGGCTGGTTTTCTGTGATGATTGACCAGCCTCAGAGTATTGATGACACCACCGTTCATAACAGATCGGCATTACTATGCATAAGTACTGTC  >  NC_000913/2052052‑2052275

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