Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,569,965	T→C	62.8%	Q191Q (CAA→CAG)	yddW ←	predicted liprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,569,965	0	T	C	62.8%	37.2 / 42.7	43	Q191Q (CAA→CAG)	yddW	predicted liprotein
Reads supporting (aligned to +/- strand):  ref base T (5/11);  new base C (10/17);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGT  >  W3110S.gb/1569899‑1570021                                                                   |                                                         cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTg                                                        >  1:2023922/1‑69 (MQ=255) cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTg                                                        >  1:1980638/1‑69 (MQ=255)   gAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGTTCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca                                                      <  1:3216822/69‑1 (MQ=255)   gAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                      <  1:209183/69‑1 (MQ=255)       cTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATa                                                  <  1:1615137/69‑1 (MQ=255)       cTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATa                                                  <  1:3541530/69‑1 (MQ=255)            ggATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                      <  1:179126/60‑1 (MQ=255)               tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCg                                          >  1:3410586/1‑69 (MQ=255)               tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCg                                          >  1:3540939/1‑69 (MQ=255)                 cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGcc                                        >  1:1131032/1‑69 (MQ=255)                 cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGcc                                        >  1:1196600/1‑69 (MQ=255)                 cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGcc                                        >  1:2964521/1‑69 (MQ=255)                       tCGAGGACAAAGCGATCGCCATACGTTCTGATCCAGTCGCGGTGTTGca                                                      <  1:2475415/49‑1 (MQ=255)                       tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCTGTGTTGca                                                      <  1:2715947/49‑1 (MQ=255)                       tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                      <  1:3492190/49‑1 (MQ=255)                       tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                      >  1:3629121/1‑49 (MQ=255)                       tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTCTTGca                                                      <  1:628301/49‑1 (MQ=255)                       tCGAGGACAAAGCGATCGCCAGAAGTTCTGATCCAGTCGCGGTGTTGc                                                       >  1:1775710/1‑48 (MQ=255)                         gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttg                                   >  1:2230124/1‑66 (MQ=255)                         gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttg                                   >  1:2408469/1‑66 (MQ=255)                              cAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttgt                                  >  1:3204391/1‑62 (MQ=255)                              cAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTgaga                             >  1:1522930/1‑67 (MQ=255)                              cAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGAACATAGCCGCTCGCCGGTTGTTgaga                             >  1:2248900/1‑67 (MQ=255)                                  gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGcaca                                                    <  1:3367525/40‑1 (MQ=255)                                  gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGcaca                                                    <  1:703791/40‑1 (MQ=255)                                  gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGcaca                                                    <  1:2875301/40‑1 (MQ=255)                                  gCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                          >  1:2762955/1‑66 (MQ=255)                                     aTCGCCAGACGTTCTGATCCAGTCGCGGTTCTGCACATa                                                  <  1:2957527/39‑1 (MQ=37)                                      tCGCCAGACGTTCTGATCCAGTCGCGGTGCTGcaca                                                    <  1:2073153/36‑1 (MQ=255)                                       cGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTAt                  <  1:3429184/69‑1 (MQ=255)                                       cGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTAt                  <  1:702557/69‑1 (MQ=255)                                           aGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAg              <  1:2517749/69‑1 (MQ=255)                                           aGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAg              <  1:1952752/69‑1 (MQ=255)                                                ttCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGAGCTATTCa               <  1:2567538/63‑1 (MQ=255)                                                 tCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt        <  1:2677822/69‑1 (MQ=255)                                                    gATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa      <  1:3329525/68‑1 (MQ=255)                                                       ccAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt  <  1:1861026/69‑1 (MQ=255)                                                       ccAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGt  <  1:1470154/69‑1 (MQ=255)                                                            cgcgGTGCTGCACATAGACGCTCGCCGGTTGTTga                               <  1:1031017/35‑1 (MQ=255)                                                             gcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTcc          <  1:1744496/55‑1 (MQ=255)                                                             gcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa      <  1:34360/59‑1 (MQ=255)                                                             gcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa      <  1:1505892/59‑1 (MQ=255)                                                             gcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa      <  1:1498203/59‑1 (MQ=255)                                                                   |                                                         CTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGT  >  W3110S.gb/1569899‑1570021

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