Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,156,064	A→C	L46R (CTA→CGA)	yheU ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,156,064	0	A	C	100.0%	210.4 / NA	62	L46R (CTA→CGA)	yheU	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (33/29);  total (33/29)

GACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGTGCC  >  W3110S.gb/4155996‑4156128                                                                     |                                                                 gactgCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGc                                                                >  1:1957995/1‑71 (MQ=255)         ggCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGttt                                                        <  1:1003792/71‑1 (MQ=255)            ataataTTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGa                                                      >  1:1949436/1‑70 (MQ=255)             taataTTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGttt                                                        <  1:812194/67‑1 (MQ=255)             taataTTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGAc                                                     >  1:345851/1‑70 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGa                                                      >  1:2477933/1‑67 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACg                                                    >  1:3174511/1‑69 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:3878113/1‑71 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:2556073/1‑71 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:1653050/1‑71 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:1992445/1‑71 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:3972227/1‑71 (MQ=255)               atatTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:1875360/1‑71 (MQ=255)               atatTGAACGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTc                                                  >  1:842111/1‑71 (MQ=255)                  ttGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCgg                                                >  1:3260857/1‑70 (MQ=255)                    gACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCg                                              <  1:2571704/70‑1 (MQ=255)                    gACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGa                                             >  1:899891/1‑71 (MQ=255)                       cGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACtt                                          >  1:884287/1‑71 (MQ=255)                        gTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACttt                                         >  1:2281097/1‑71 (MQ=255)                         tttCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGa                                             <  1:1429195/66‑1 (MQ=255)                          ttCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTc                                        >  1:4093525/1‑70 (MQ=255)                          ttCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCt                                       >  1:2309777/1‑71 (MQ=255)                          ttCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCt                                       >  1:853383/1‑71 (MQ=255)                           tCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTg                                      >  1:2012710/1‑71 (MQ=255)                           tCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTg                                      >  1:971925/1‑71 (MQ=255)                           tCGTGCAGCTCCGACCAAACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTg                                      >  1:1574154/1‑71 (MQ=255)                                cAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCa                                  <  1:4259224/70‑1 (MQ=255)                                cAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCa                                  <  1:3545234/70‑1 (MQ=255)                                cAGCTCCGACCATACCAGCACAGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTc                                   >  1:4011391/1‑69 (MQ=255)                                   cTCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAgtg                              <  1:4479349/71‑1 (MQ=255)                                    tCCGACCATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAgtg                              >  1:4660030/1‑70 (MQ=255)                                         ccATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCt                        <  1:258005/71‑1 (MQ=255)                                          cATACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCt                        >  1:3173319/1‑70 (MQ=255)                                            tACCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCt                        <  1:1537418/68‑1 (MQ=255)                                             aCCAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATg                    >  1:375093/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCa                      >  1:3493333/1‑68 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:1155510/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:327294/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:872593/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:4352023/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:4295809/1‑71 (MQ=255)                                              ccAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGt                   >  1:3906192/1‑71 (MQ=255)                                               cAGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGtt                  <  1:3663041/71‑1 (MQ=255)                                                aGCACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTc                 <  1:4391872/71‑1 (MQ=255)                                                  cACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCAc               <  1:958206/71‑1 (MQ=255)                                                  cACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCAc               <  1:1756244/71‑1 (MQ=255)                                                  cACCGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCAc               <  1:45555/71‑1 (MQ=255)                                                     cGCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCAt            <  1:3880159/71‑1 (MQ=255)                                                      gCTTCTCCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCAt            <  1:406257/70‑1 (MQ=255)                                                          ctcCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATc        <  1:4639080/70‑1 (MQ=255)                                                          ctcCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATc        <  1:1563258/70‑1 (MQ=255)                                                          ctcCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATc        <  1:1664665/70‑1 (MQ=255)                                                           tcCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATg                    <  1:3768876/57‑1 (MQ=255)                                                           tcCGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATc        <  1:4053578/69‑1 (MQ=255)                                                             cGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCgg      <  1:773082/69‑1 (MQ=255)                                                             cGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGt     <  1:1437778/70‑1 (MQ=255)                                                             cGCACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGt     <  1:2942277/70‑1 (MQ=255)                                                               cACTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGTGcc  <  1:1215462/71‑1 (MQ=255)                                                                 cTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCAc               <  1:4020748/56‑1 (MQ=255)                                                                 cTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCAc               <  1:1536789/56‑1 (MQ=255)                                                                 cTGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGt     <  1:1667701/66‑1 (MQ=255)                                                                  tGTCGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCgg      <  1:4101868/64‑1 (MQ=255)                                                                     |                                                                 GACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGTTCACCATAATCGGTGCC  >  W3110S.gb/4155996‑4156128

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