Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,341,293	+T	coding (22/1689 nt)	yfaA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,341,290	1	.	T	98.4%	174.6 / ‑3.2	61	coding (25/1689 nt)	yfaA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (21/39);  total (22/39)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑TTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTT  >  W3110S.gb/2341221‑2341347                                                                        |                                                          gCAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑t                                                          >  1:3579465/1‑71 (MQ=255) gCAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑t                                                          >  1:664516/1‑71 (MQ=255)  cAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑tt                                                         <  1:941332/71‑1 (MQ=255)  cAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑tt                                                         <  1:3493005/71‑1 (MQ=255)   aaCGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑tt                                                         <  1:326753/70‑1 (MQ=255)         cAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCt                                                   >  1:3381698/1‑71 (MQ=255)           gACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCttt                                                 >  1:93737/1‑71 (MQ=255)               ccAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCttttttt                                             <  1:1042273/71‑1 (MQ=255)                cAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTc                                            <  1:3423279/71‑1 (MQ=255)                  gCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcac                                          <  1:1001192/71‑1 (MQ=255)                  gCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcac                                          <  1:2182268/71‑1 (MQ=255)                   cGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcacc                                         <  1:4302058/71‑1 (MQ=255)                      aaaGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:1281444/71‑1 (MQ=255)                      aaaGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:2545282/71‑1 (MQ=255)                      aaaGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:170167/71‑1 (MQ=255)                      aaaGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:4681842/71‑1 (MQ=255)                      aaaGCAGTGCCATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:2007811/71‑1 (MQ=255)                          cAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCAtt                                  <  1:279457/71‑1 (MQ=255)                          cAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑tt                                                         >  1:3726161/1‑47 (MQ=255)                            gTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:2015289/65‑1 (MQ=255)                            gTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:1637326/65‑1 (MQ=255)                             tGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCtttt                                                >  1:4345725/1‑54 (MQ=255)                             tGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCa                                    <  1:2134063/66‑1 (MQ=255)                              gCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:4557807/1‑71 (MQ=255)                              gCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:1701842/1‑71 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcac                                          >  1:2108346/1‑58 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcac                                          <  1:1063177/58‑1 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCa                                    >  1:1162138/1‑64 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:3755875/1‑70 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:4189351/1‑70 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:4566024/1‑70 (MQ=255)                               cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              <  1:3210119/70‑1 (MQ=255)                               cTATTGCGCCAAAACCTACAAGAACATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:232769/1‑70 (MQ=255)                                 aTTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcacca                                        >  1:1123515/1‑58 (MQ=255)                                  ttGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCa                                    <  1:3480531/61‑1 (MQ=255)                                    gcgcCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCcaca                            >  1:1814631/1‑67 (MQ=255)                                    gcgcCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:3176043/1‑65 (MQ=255)                                    gcgcCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCGTTTTTTGACCACTCATTGCcac                             >  1:2571656/1‑66 (MQ=255)                                     cgcCAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTc                         <  1:2475086/69‑1 (MQ=255)                                       ccAAAACCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCttttttt                                             <  1:1902891/47‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:326343/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:3215384/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:2126924/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:3035333/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:4344990/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:2319969/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:412667/71‑1 (MQ=255)                                         aaaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:3942333/71‑1 (MQ=255)                                          aaaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:330559/70‑1 (MQ=255)                                           aaCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGta                 <  1:2984284/71‑1 (MQ=255)                                            aCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGtat                <  1:2437647/71‑1 (MQ=255)                                            aCCTACAAGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTAATTGCCACATTCCTTGTGta                 >  1:4640871/1‑70 (MQ=255)                                             ccTACAAGACCATAGAACCGCCAGGCT‑TTTCGCCTTTTTTTCACCACTCATTGCca                              >  1:3244020/1‑56 (MQ=255)                                                   aGACCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATagccagc         <  1:793250/71‑1 (MQ=255)                                                     aCCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTcacc                                         >  1:3815770/1‑37 (MQ=38)                                                     aCCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:1035953/59‑1 (MQ=255)                                                     aCCATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAgccagcca       <  1:3908404/71‑1 (MQ=255)                                                      ccATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAgccagcca       <  1:3106018/70‑1 (MQ=255)                                                       cATAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgtg                   <  1:4478343/57‑1 (MQ=255)                                                         tAGAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACt                                      <  1:3935842/36‑1 (MQ=37)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTc                                     <  1:1828995/35‑1 (MQ=37)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  >  1:364404/1‑70 (MQ=255)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  >  1:3076027/1‑70 (MQ=255)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  <  1:3436188/70‑1 (MQ=255)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  >  1:1728311/1‑70 (MQ=255)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  <  1:1063376/70‑1 (MQ=255)                                                           gAACCGCCA‑GCCTTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt   >  1:1348534/1‑69 (MQ=255)                                                                        |                                                          GCAACGCCCAGACGCCAGCGGAAAGCAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCA‑GCC‑TTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTT  >  W3110S.gb/2341221‑2341347

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