Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,936	G→A	16.6%	R91R (CGC→CGT)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,936	0	G	A	16.6%	93.5 / 10.9	42	R91R (CGC→CGT)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (23/12);  new base A (3/4);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAC  >  W3110S.gb/4563869‑4564004                                                                    |                                                                     ccGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGa                                                                   >  1:2165358/1‑71 (MQ=255)    ggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCTGCTTCCGTCAGACGACtt                                                                <  1:3836108/71‑1 (MQ=255)    ggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCTGCTTCCGTCAGACGACtt                                                                <  1:1427697/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:3739023/1‑71 (MQ=255)                           cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                                         >  1:2039097/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:3201928/1‑71 (MQ=255)                                            tgacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgc                        >  1:3206174/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)                                                             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)                                                                    |                                                                     CCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAC  >  W3110S.gb/4563869‑4564004

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