Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,223,904	C→T	Y424Y (TAC→TAT)	ebgA →	evolved beta‑D‑galactosidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,223,904	0	C	T	100.0%	133.1 / NA	37	Y424Y (TAC→TAT)	ebgA	evolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/22);  total (15/22)

GACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATG  >  NC_000913/3223827‑3223977                                                                              |                                                                          gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTc                                                                        >  1:606787/1‑81 (MQ=255) gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTc                                                                        >  1:606775/1‑81 (MQ=255)  aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCg                                                                       <  1:606723/81‑1 (MQ=255)    gTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAg                                                                     >  1:606796/1‑81 (MQ=255)     tCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAgc                                                                    <  2:606724/81‑1 (MQ=255)      cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAgcg                                                                   >  2:606770/1‑81 (MQ=255)       gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAgcgc                                                                  <  1:641232/81‑1 (MQ=255)        aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCa                                                                 <  1:606725/81‑1 (MQ=255)          tCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCAtt                                                               <  1:606726/81‑1 (MQ=255)           cGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCAtcg                                                              >  2:606793/1‑79 (MQ=255)            gCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGt                                                             <  1:606727/81‑1 (MQ=255)              aCGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTc                                                           >  2:606777/1‑81 (MQ=255)              aCGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTc                                                           <  1:606728/81‑1 (MQ=255)                ggCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGc                                                         >  2:606772/1‑81 (MQ=255)                 gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGcc                                                        <  1:606729/81‑1 (MQ=255)                     tGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCat                                                      <  2:606730/79‑1 (MQ=255)                       cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATc                                                   <  2:606731/80‑1 (MQ=255)                        tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCa                                                 >  2:606792/1‑81 (MQ=255)                         aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCAc                                                >  1:606804/1‑81 (MQ=255)                          aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCAc                                                <  1:606732/80‑1 (MQ=255)                           tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCAcgc                                              <  2:606733/81‑1 (MQ=255)                             tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCAcgcg                                             >  1:606788/1‑80 (MQ=255)                                 cGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGa                                         >  1:606774/1‑80 (MQ=255)                                    tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGaaaaa                                     <  1:606734/81‑1 (MQ=255)                                      ttAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAAc                                    <  1:606735/80‑1 (MQ=255)                                         gCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCAt                                 <  2:606736/80‑1 (MQ=255)                                          ccGTATTACCGACGATCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATcc                               <  1:606737/81‑1 (MQ=255)                                                    gacgaTCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGatcatc                     <  1:606738/81‑1 (MQ=255)                                                    gacgaTCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGatcat                      >  1:606799/1‑80 (MQ=255)                                                    gacgaTCCGCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGatcat                      >  1:606791/1‑80 (MQ=255)                                                            gCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTc             <  2:606739/81‑1 (MQ=255)                                                            gCAGTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTc             >  2:606790/1‑81 (MQ=255)                                                               gTGGGAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt          >  1:606789/1‑81 (MQ=255)                                                                   gAAAAAGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGc      <  2:606740/81‑1 (MQ=255)                                                                      aaaGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAAt   <  1:606742/81‑1 (MQ=255)                                                                      aaaGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAAt   <  1:606741/81‑1 (MQ=255)                                                                       aaGTCTATGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATg  <  2:606743/81‑1 (MQ=255)                                                                              |                                                                          GACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATG  >  NC_000913/3223827‑3223977

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