Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,944	G→A	12.8%	L89L (CTA→TTA)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,944	0	G	A	12.8%	95.8 / 6.0	39	L89L (CTA→TTA)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (23/11);  new base A (3/2);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGAT  >  W3110S.gb/4563874‑4564012                                                                       |                                                                     tGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAg                                                                      <  1:3075428/71‑1 (MQ=255)     gagagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGcc                                                                  <  1:2355513/71‑1 (MQ=255)     gagagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCCCTTCCGTCAGGCGACTTAGCGcc                                                                  <  1:3534175/71‑1 (MQ=255)      agagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCa                                                                 >  1:847508/1‑71 (MQ=255)      agagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCa                                                                 >  1:197896/1‑71 (MQ=255)      agagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCa                                                                 >  1:1837504/1‑71 (MQ=255)                    tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                                                   <  1:3294010/71‑1 (MQ=255)                    tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                                                   <  1:1800280/71‑1 (MQ=255)                      cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                                                 >  1:2039097/1‑71 (MQ=255)                      cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                                                 >  1:2155149/1‑71 (MQ=255)                      cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                                                 >  1:3739023/1‑71 (MQ=255)                        cagcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgtcg                                               <  1:1295124/71‑1 (MQ=255)                         agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgtcgt                                              <  1:777620/71‑1 (MQ=255)                          gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTc                                                             >  1:3261133/1‑55 (MQ=255)                          gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTc                                                             >  1:1305982/1‑55 (MQ=255)                              accaaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAc                                         <  1:2627025/71‑1 (MQ=255)                               ccaaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAcc                                        >  1:2193455/1‑71 (MQ=255)                                 aaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                     >  1:2582159/1‑40 (MQ=255)                                  accaCTGACGTGACGCCTGCTTCCGTCAGACGACTTAACGCCACTTCCGGCGTGcgcg                                                  >  1:1320913/1‑58 (MQ=255)                                  accaCTGACGTGACGCCTGCTTCCGTCAGACGACTTAACGCCACTTCCGGCGTGCGCGTCGTGGGACGTGc                                     >  1:1177953/1‑71 (MQ=255)                                   ccaCTGACGTGACGCCTGCTTCCGTCAGACGACTTAACGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCt                                    >  1:1590027/1‑71 (MQ=255)                                      cTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcg                                 >  1:3600727/1‑71 (MQ=255)                                       tgacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgc                                >  1:3206174/1‑71 (MQ=255)                                       tgacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgc                                >  1:3201928/1‑71 (MQ=255)                                         acgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccg                              >  1:1980605/1‑71 (MQ=255)                                         acgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccg                              >  1:2740521/1‑71 (MQ=255)                                          cgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccgc                             <  1:1228987/71‑1 (MQ=255)                                           gtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccgc                             >  1:3826811/1‑70 (MQ=255)                                                        ccGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCaaa               >  1:2551404/1‑71 (MQ=255)                                                        ccGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCa                 <  1:939900/69‑1 (MQ=255)                                                        ccGTCAGACGACTTAACGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCa                 <  1:1602673/69‑1 (MQ=255)                                                        ccGTCAGACGACTTAACGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCa                 <  1:975642/69‑1 (MQ=255)                                                          gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATg             >  1:3573864/1‑71 (MQ=255)                                                            cAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCACCAATCAAATGGa           >  1:2840383/1‑71 (MQ=255)                                                             aGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAc          >  1:2666550/1‑71 (MQ=255)                                                             aGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAc          >  1:1814741/1‑71 (MQ=255)                                                                   cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg    <  1:259384/71‑1 (MQ=255)                                                                    ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa   >  1:2656434/1‑71 (MQ=255)                                                                    ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa   >  1:408418/1‑71 (MQ=255)                                                                      aGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa   <  1:2968858/69‑1 (MQ=255)                                                                      aGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGAt  >  1:1950444/1‑70 (MQ=255)                                                                       |                                                                     TGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGAT  >  W3110S.gb/4563874‑4564012

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: ‑