Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,093,680	G→A	10.2%	H19H (CAC→CAT)	yggR ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,093,680	0	G	A	10.2%	141.1 / 5.5	49	H19H (CAC→CAT)	yggR	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (31/13);  new base A (2/3);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACT  >  W3110S.gb/3093611‑3093739                                                                      |                                                            cGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAgg                                                             >  1:2933524/1‑70 (MQ=255)    cGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                                                            >  1:510341/1‑68 (MQ=255)        tCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATc                                                     >  1:2158085/1‑71 (MQ=255)         cAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTa                                                         >  1:1785660/1‑66 (MQ=255)          aaccGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATcc                                                    <  1:882912/67‑1 (MQ=255)          aaaCGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATcc                                                    <  1:3266790/70‑1 (MQ=255)           aaCGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCga                                                  <  1:317581/71‑1 (MQ=255)           aaCGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGATGAAGATCCga                                                  <  1:3246023/71‑1 (MQ=255)             cGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCgaga                                                >  1:244933/1‑71 (MQ=255)               tcgcAGCTTCCATTCTGCCGCGAATGCGCCATCGTGC‑GGCCAGGCGCTGCACAGGTGTAGATCCGAGACGt                                             <  1:3671884/70‑1 (MQ=255)                 gcAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgtgt                                                          <  1:2061029/57‑1 (MQ=255)                   aGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATcc                                                    <  1:1632873/61‑1 (MQ=255)                   aGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTAt                                          >  1:4010114/1‑71 (MQ=255)                       tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTTTAGATCCGAGACGTTATGCtt                                      >  1:17865/1‑71 (MQ=255)                       tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      >  1:4338065/1‑71 (MQ=255)                       tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      >  1:3682369/1‑71 (MQ=255)                       tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      >  1:4457326/1‑71 (MQ=255)                       tCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      >  1:3744830/1‑71 (MQ=255)                        ccATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      <  1:227036/70‑1 (MQ=255)                         cATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCttt                                     >  1:3752535/1‑70 (MQ=255)                         cATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa                                    >  1:409766/1‑71 (MQ=255)                         cATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa                                    >  1:1180558/1‑71 (MQ=255)                           ttCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATc                                                     >  1:2305875/1‑52 (MQ=255)                            tCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAgg                                                             <  1:928426/43‑1 (MQ=255)                             cTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTacac                                 >  1:4279386/1‑70 (MQ=255)                             cTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTacac                                 >  1:3933027/1‑70 (MQ=255)                             cTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTacac                                 >  1:963427/1‑70 (MQ=255)                                  gcgAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCt                                       >  1:1100260/1‑59 (MQ=255)                                   cgAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGATGTAGATCCGAGACGTTATGCTTTACACTAAGGGc                          >  1:1831260/1‑71 (MQ=255)                                      aTGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCcac                       <  1:1105943/71‑1 (MQ=255)                                         cgcCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCtt                                      <  1:2810291/53‑1 (MQ=255)                                         cgcCATCGTGCGGGCCAGGCGCTGCACAGATGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACaa                     <  1:734174/70‑1 (MQ=255)                                              tCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttctt               >  1:550746/1‑71 (MQ=255)                                               cGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttcttc              >  1:3361489/1‑71 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCTAGACGTTATGCTTTACACTAAGGGCCACAAttt                  >  1:811786/1‑66 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTaca                                  >  1:206176/1‑50 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa                                    >  1:4574095/1‑48 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAg                             >  1:3336600/1‑55 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAg                             >  1:2051559/1‑55 (MQ=255)                                                gTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttcttc              >  1:334085/1‑70 (MQ=255)                                                gTGCGGGCCAGGCGCTGAACAGGTGTAGATCCGAGACGt                                             >  1:4119638/1‑39 (MQ=255)                                                   cGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCAAAATTTCTTCCat           <  1:2512496/70‑1 (MQ=255)                                                     ggCCAGGCGCTGCACAGATGTAGATCCGAGACGTTATGCTTTACACTAAGGGCcaca                      <  1:4025908/57‑1 (MQ=255)                                                      gCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTc       <  1:1617072/71‑1 (MQ=255)                                                        cAGGCGCTGCACAGATGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCa      >  1:1126427/1‑70 (MQ=255)                                                         aGGCGCTGCACAGGTGTAGATCCGAGACGTTATGc                                        <  1:1216796/35‑1 (MQ=255)                                                             gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCa      >  1:1390802/1‑65 (MQ=255)                                                             gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCa      >  1:4063903/1‑65 (MQ=255)                                                             gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACt  >  1:1476281/1‑69 (MQ=255)                                                             gcTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACt  >  1:967439/1‑69 (MQ=255)                                                               tGCACATGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATttc                 <  1:279875/52‑1 (MQ=255)                                                               tGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACt  >  1:1201433/1‑67 (MQ=255)                                                                 cacaGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCatat         <  1:1191700/58‑1 (MQ=255)                                                                      gtgtAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCa            <  1:3744335/50‑1 (MQ=255)                                                                      |                                                            CGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACACTAAGGGCCACAATTTCTTCCATATTCATACT  >  W3110S.gb/3093611‑3093739

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