Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	494,006	T→A	9.1%	D126E (GAT→GAA)	recR →	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	494,006	0	T	A	9.1%	182.9 / 9.9	66	D126E (GAT→GAA)	recR	gap repair protein
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base A (4/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGA  >  W3110S.gb/493939‑494074                                                                    |                                                                     gTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCg                                                                    <  1:1449359/70‑1 (MQ=255)                  gggCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTgg                                                   >  1:3318594/1‑70 (MQ=255)                  gggCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTgg                                                   >  1:3257463/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:1786799/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:127174/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:351664/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:1452116/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:3248187/1‑70 (MQ=255)                         tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                            >  1:1768366/1‑70 (MQ=255)                          gTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGa                                                              >  1:1211589/1‑51 (MQ=255)                          gTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGa                                                              >  1:1356677/1‑51 (MQ=255)                          gTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGa                                                              >  1:850439/1‑51 (MQ=255)                           tCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTgg                                                   <  1:397295/61‑1 (MQ=255)                             aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                       >  1:1363854/1‑71 (MQ=255)                             aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                       >  1:3328931/1‑71 (MQ=255)                             aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                       >  1:2469520/1‑71 (MQ=255)                             aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                        <  1:3711806/70‑1 (MQ=255)                                    gACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                                 <  1:28070/70‑1 (MQ=255)                                    gACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                                 <  1:1239258/70‑1 (MQ=255)                                    gACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                                 <  1:2574618/70‑1 (MQ=255)                                    gACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGa                                <  1:2390814/71‑1 (MQ=255)                                    gACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGa                                <  1:759997/71‑1 (MQ=255)                                           tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc                         >  1:807495/1‑71 (MQ=255)                                           tCGGTCCGGATGATATCGGGCTTGATCGTCTCGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc                         >  1:1402260/1‑71 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:348357/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:905415/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:318438/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:3171685/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:284412/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:3075369/1‑68 (MQ=255)                                             ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                       >  1:1091643/1‑71 (MQ=255)                                              gTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATcc                        <  1:2917863/69‑1 (MQ=255)                                              gTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATcc                        <  1:2508339/69‑1 (MQ=255)                                                ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc                      >  1:2515860/1‑69 (MQ=255)                                                       atatCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAc             >  1:1164905/1‑71 (MQ=255)                                                       atatCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAc             >  1:376000/1‑71 (MQ=255)                                                         atCGGGCTTGATCGTCTGGAACAGCGTCTGGGAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAccc           >  1:1019325/1‑71 (MQ=255)                                                         atCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAccc           >  1:728725/1‑71 (MQ=255)                                                         atCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAcc            >  1:276338/1‑70 (MQ=255)                                                          tCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGa                                            <  1:1754248/37‑1 (MQ=255)                                                            gggCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCAc        <  1:973993/71‑1 (MQ=255)                                                            gggCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCAc        >  1:744434/1‑71 (MQ=255)                                                            gggCTTGATCGTCTGGAACAGCGTCGCGCAGAGGAAAAAAt                                      >  1:3280744/1‑41 (MQ=38)                                                             ggCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                            <  1:3132066/50‑1 (MQ=255)                                                             ggCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                            <  1:2465826/50‑1 (MQ=255)                                                               cTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:87603/1‑50 (MQ=255)                                                               cTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                          >  1:1544415/1‑50 (MQ=255)                                                               cTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                            <  1:3199491/48‑1 (MQ=255)                                                               cTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                            <  1:830052/48‑1 (MQ=255)                                                                ttGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACg       >  1:107451/1‑68 (MQ=255)                                                                ttGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGtt    >  1:3277554/1‑71 (MQ=255)                                                                ttGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGtt    >  1:846040/1‑71 (MQ=255)                                                                ttGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGtt    >  1:2303295/1‑71 (MQ=255)                                                                  gATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTg   <  1:2179485/70‑1 (MQ=255)                                                                  gATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTg   <  1:832681/70‑1 (MQ=255)                                                                  gATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTg   <  1:87497/70‑1 (MQ=255)                                                                  gATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTg   <  1:1133823/70‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                            <  1:600928/44‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                           <  1:405382/45‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                           <  1:2195335/45‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                           <  1:1799248/45‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                           <  1:400906/45‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGa  <  1:3578219/70‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGa  <  1:152852/70‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGa  <  1:3227210/70‑1 (MQ=255)                                                                   aTCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGa  <  1:3414781/70‑1 (MQ=255)                                                                    |                                                                     GTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGA  >  W3110S.gb/493939‑494074

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