Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	542,604	(T)6→5	coding (667/786 nt)	ylbA ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	542,599	0	T	.	100.0%	183.9 / NA	46	coding (672/786 nt)	ylbA	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (29/17);  total (29/17)

CACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGG  >  W3110S.gb/542550‑542660                                                  |                                                              cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGtt                                           <  1:4164362/69‑1 (MQ=255) cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGtt                                           <  1:2099766/69‑1 (MQ=255) cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTAt                                         >  1:1490110/1‑71 (MQ=255) cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTAt                                         >  1:1495688/1‑71 (MQ=255) cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTAt                                         >  1:1660977/1‑71 (MQ=255) cacCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTAt                                         >  1:1873833/1‑71 (MQ=255)        aCCAGCCTGTAAAGAATAAGCGCTCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAg                                  >  1:930783/1‑71 (MQ=255)        aCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTAt                                         >  1:2664261/1‑64 (MQ=255)        aCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAg                                  >  1:89246/1‑71 (MQ=255)        aCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAg                                  >  1:3531322/1‑71 (MQ=255)        aCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAg                                  >  1:2223235/1‑71 (MQ=255)          cAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGt                                <  1:2295868/71‑1 (MQ=255)             ccTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTta                              >  1:1201824/1‑70 (MQ=255)             ccTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTta                              >  1:33478/1‑70 (MQ=255)             ccTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTta                              >  1:1124843/1‑70 (MQ=255)              cTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCACTTATTGTCGAGGTtat                             >  1:2104183/1‑70 (MQ=255)                   aaGAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAAccc                       >  1:2556458/1‑71 (MQ=255)                     gAATAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAAccccc                     >  1:3092062/1‑71 (MQ=255)                      aaTAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGt                                <  1:3155832/59‑1 (MQ=255)                        tAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGa                  >  1:3479395/1‑71 (MQ=255)                        tAAGCGCCCATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGa                  >  1:2930372/1‑71 (MQ=255)                                cATAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAAcccc                      >  1:1960744/1‑59 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGt                                <  1:2751178/47‑1 (MQ=39)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGt                                <  1:1951991/47‑1 (MQ=39)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:2696885/71‑1 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:4219684/71‑1 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:276014/71‑1 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:88246/71‑1 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:1368826/71‑1 (MQ=255)                                  tAAAGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:39786/71‑1 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      <  1:502332/71‑1 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:1284634/1‑71 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:39640/1‑71 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:3283577/1‑71 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      <  1:1900386/71‑1 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:3343561/1‑71 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      <  1:2678320/71‑1 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:1969055/1‑71 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAat       >  1:1513521/1‑70 (MQ=255)                                    aaGATGTAATCGCC‑TTTTTCACAGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:2708443/1‑71 (MQ=255)                                        tGTAATCGCC‑TTTTTCACCGGGATCCATTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAgg  <  1:3107596/71‑1 (MQ=255)                                        tGTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGaa        <  1:4326766/65‑1 (MQ=255)                                         gTAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAgg  >  1:3733582/1‑70 (MQ=255)                                          tAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAAc                         >  1:2507327/1‑46 (MQ=39)                                          tAATCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata      >  1:4149062/1‑65 (MQ=255)                                            aTCGCC‑TTTTTCACCGGGATCCAGTTATTGTCGAGGTTATaaa                          >  1:2614730/1‑43 (MQ=38)                                                  |                                                              CACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGG  >  W3110S.gb/542550‑542660

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