Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	363,547	A→G	17.8%	G661G (GGT→GGC)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	363,547	0	A	G	17.8%	55.9 / 6.9	28	G661G (GGT→GGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (13/10);  new base G (3/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTG  >  W3110S.gb/363478‑363610                                                                      |                                                                aGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTa                                                                 <  1:3559586/70‑1 (MQ=255)              aaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGccac                                                   >  1:3598050/1‑71 (MQ=255)                   tttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATc                                               <  1:924493/70‑1 (MQ=255)                     tACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTAc                                                                >  1:3310080/1‑51 (MQ=255)                         ttGTGGAGCGACATCGAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGc                                         <  1:218337/70‑1 (MQ=255)                         ttGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTg                                          >  1:3592426/1‑69 (MQ=255)                         ttGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGc                                         >  1:1069537/1‑70 (MQ=255)                                 cGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCg                                >  1:3700842/1‑71 (MQ=255)                                  gACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCa                                             >  1:3029142/1‑57 (MQ=255)                                    cATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTa                             <  1:2998478/71‑1 (MQ=255)                                    cATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTa                             <  1:3613852/71‑1 (MQ=255)                                        cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATc                                               >  1:1790683/1‑49 (MQ=255)                                        cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTAt                            >  1:3803948/1‑68 (MQ=255)                                                 ttCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGCAGGAACTCGTTATCGCTATGACgg                 >  1:1419361/1‑70 (MQ=255)                                                 ttCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGCAGGAACTCGTTATCGCTATGACgg                 >  1:1344680/1‑70 (MQ=255)                                                 ttCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGCAGGAACTCGTTATCGCTATGACgg                 >  1:3628319/1‑70 (MQ=255)                                                 ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCAt                                                <  1:2587699/39‑1 (MQ=255)                                                   cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa               >  1:2337708/1‑70 (MQ=255)                                                   cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCAACCAGTGCAGGAGCTCGTTATCGCTATGACGGaa               >  1:2798856/1‑70 (MQ=255)                                                    aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACa             <  1:3839739/71‑1 (MQ=255)                                                     ccGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGCAGGAACTCGTTATCGCTATGACGGaa               <  1:491466/68‑1 (MQ=255)                                                     ccGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGCAGGAACTCGTTATCGCTATGACGGaa               <  1:1603659/68‑1 (MQ=255)                                                     ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCg                          <  1:1924078/57‑1 (MQ=255)                                                      cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAgg           <  1:472801/71‑1 (MQ=255)                                                      cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAgg           <  1:1035945/71‑1 (MQ=255)                                                       gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTa                             >  1:2463070/1‑52 (MQ=255)                                                          tGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAAc              <  1:1877245/64‑1 (MQ=255)                                                           gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTc      >  1:1189472/1‑71 (MQ=255)                                                               gCGGCTTACCATCCAGCTCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTg  >  1:1989846/1‑71 (MQ=255)                                                                      |                                                                AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTG  >  W3110S.gb/363478‑363610

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