Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,924	C→T	14.3%	A95A (GCG→GCA)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,924	0	C	T	14.3%	83.5 / 5.7	35	A95A (GCG→GCA)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (3/2);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGC  >  W3110S.gb/4563862‑4563985                                                               |                                                              aGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGt                                                       <  1:853778/71‑1 (MQ=255)        ccGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGa                                                >  1:2165358/1‑71 (MQ=255)           ggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCTGCTTCCGTCAGACGACtt                                             <  1:1427697/71‑1 (MQ=255)           ggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCTGCTTCCGTCAGACGACtt                                             <  1:3836108/71‑1 (MQ=255)           ggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACt                                              >  1:3693054/1‑70 (MQ=255)             tGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAg                                           <  1:3075428/71‑1 (MQ=255)                cGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCg                                                 <  1:2076976/62‑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:1837504/1‑71 (MQ=255)                  agagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCa                                      >  1:197896/1‑71 (MQ=255)                      aTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTa                                            <  1:91332/61‑1 (MQ=255)                                tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                        <  1:1800280/71‑1 (MQ=255)                                tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                        <  1:3294010/71‑1 (MQ=255)                                  cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                      >  1:2155149/1‑71 (MQ=255)                                  cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                      >  1:2039097/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:1305982/1‑55 (MQ=255)                                      gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTc                                  >  1:3261133/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:2740521/1‑71 (MQ=255)                                                     acgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccg   >  1:1980605/1‑71 (MQ=255)                                                      cgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccgc  <  1:1228987/71‑1 (MQ=255)                                                       gtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccgc  >  1:3826811/1‑70 (MQ=255)                                                               |                                                              AGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGC  >  W3110S.gb/4563862‑4563985

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