Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	405,632	A→G	9.6%	T2A (ACA→GCA)	aroL →	shikimate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	405,632	0	A	G	9.6%	166.0 / 8.2	62	T2A (ACA→GCA)	aroL	shikimate kinase II
Reads supporting (aligned to +/- strand):  ref base A (20/36);  new base G (3/3);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGATCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGA  >  W3110S.gb/405566‑405699                                                                   |                                                                    tATGATCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGacac                                                                  >  1:4340956/1‑70 (MQ=255)  aTGATCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACaa                                                                >  1:2958340/1‑71 (MQ=255)  aTGATCGCTATTCTCATGACACAGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACaa                                                                >  1:819219/1‑71 (MQ=255)     aTCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACaa                                                                >  1:1265574/1‑68 (MQ=255)     aTCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACct                                                             >  1:2381412/1‑71 (MQ=255)     aTCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACct                                                             <  1:47424/71‑1 (MQ=255)          agatctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCtttt                                                        <  1:4168324/68‑1 (MQ=255)            ttCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGGCACAACCTCtttt                                                        <  1:2095934/69‑1 (MQ=255)            ttCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCttttt                                                       <  1:4575209/70‑1 (MQ=255)            ttCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCttttt                                                       <  1:2062798/70‑1 (MQ=255)              ctcATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAAc                                                               >  1:932686/1‑60 (MQ=255)                 aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:1830695/68‑1 (MQ=255)                 aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:2068697/68‑1 (MQ=255)                 aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:2423424/68‑1 (MQ=255)                 aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                  <  1:200157/70‑1 (MQ=255)                 aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATCACACAACCTCTTTTTCTGAt                                                  <  1:3315491/70‑1 (MQ=255)                   gACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGa                                                                     >  1:2636403/1‑49 (MQ=255)                   gACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                  >  1:3054415/1‑68 (MQ=255)                    acacCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCggg                                              >  1:2347488/1‑71 (MQ=255)                       ccGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:1591881/62‑1 (MQ=255)                       ccGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCt                                           <  1:3643485/71‑1 (MQ=255)                         ggCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCg                                         >  1:865217/1‑71 (MQ=255)                           cTTTCGCCGCATTGCGACCTATTGGCGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCggg                                       <  1:4359745/71‑1 (MQ=255)                                  cgcATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCt                                                     <  1:3943698/50‑1 (MQ=255)                                      ttGCGACCTATTGGGGAAAACCCACGATGACACAACCTCttt                                                         <  1:4544737/42‑1 (MQ=255)                                      ttGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:488105/47‑1 (MQ=255)                                      ttGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                             >  1:2226760/1‑70 (MQ=255)                                      ttGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                             <  1:3867929/70‑1 (MQ=255)                                        gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa                           >  1:262887/1‑70 (MQ=255)                                        gCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaaaa                           >  1:345739/1‑70 (MQ=255)                                         cGACCTATTGTGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaac                          <  1:3615498/70‑1 (MQ=255)                                         cGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaaca                         <  1:2159522/71‑1 (MQ=255)                                         cGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaac                          <  1:990249/70‑1 (MQ=255)                                         cGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaac                          <  1:1115369/70‑1 (MQ=255)                                         cGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaac                          <  1:2825101/70‑1 (MQ=255)                                           aCCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGt                               >  1:3108157/1‑63 (MQ=255)                                             cTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTg                                                    <  1:1725928/40‑1 (MQ=255)                                               aTTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACg                      >  1:2742572/1‑68 (MQ=255)                                                ttGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCg                                         <  1:507531/48‑1 (MQ=255)                                                ttGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaaca                         <  1:3285725/64‑1 (MQ=255)                                                 tGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                  >  1:2712211/1‑38 (MQ=255)                                                 tGGGGAAAACCCACGATGACACAACCTCTTTTTCTGAt                                                  >  1:42628/1‑38 (MQ=255)                                                  ggggAAAACCCACGATGGCACAACCTCTTTTTCTGAt                                                  >  1:1513819/1‑37 (MQ=255)                                                  ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:3690590/70‑1 (MQ=255)                                                  ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:1448199/70‑1 (MQ=255)                                                  ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:2631203/70‑1 (MQ=255)                                                  ggggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:3125126/70‑1 (MQ=255)                                                   gggAAAACCCACGATGGCACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:4310049/69‑1 (MQ=255)                                                   gggAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTgg                                <  1:3982913/54‑1 (MQ=255)                                                     gAAAACCCACGATGGCACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaaca                         <  1:1417492/59‑1 (MQ=255)                                                     gAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATg             >  1:3112378/1‑71 (MQ=255)                                                     gAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATg             >  1:1095487/1‑71 (MQ=255)                                                       aaaCCCACGATGGCACAACCTCTTATTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACCGTCGGAATGGc           >  1:2634412/1‑71 (MQ=255)                                                       aaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGc           >  1:3347985/1‑71 (MQ=255)                                                       aaaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGc           <  1:2458823/71‑1 (MQ=255)                                                        aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                             <  1:4113095/52‑1 (MQ=255)                                                        aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTaa                             <  1:2579186/52‑1 (MQ=255)                                                        aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGcc          >  1:3015390/1‑71 (MQ=255)                                                        aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGcc          <  1:3109489/71‑1 (MQ=255)                                                         aCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCgg                 <  1:1495675/63‑1 (MQ=255)                                                          cccACGATGACACAACATCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGccc         <  1:749855/70‑1 (MQ=255)                                                           ccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAaacaa                        <  1:2758496/54‑1 (MQ=255)                                                               gATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGcc    >  1:3964059/1‑70 (MQ=255)                                                                aTGGCACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGa  >  1:1960904/1‑71 (MQ=255)                                                                   |                                                                    TATGATCGCTATTCTCATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGA  >  W3110S.gb/405566‑405699

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