Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	798,224	C→T	6.1%	A210A (GCG→GCA)	ybhA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	798,224	0	C	T	6.1%	261.4 / 4.6	82	A210A (GCG→GCA)	ybhA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base T (3/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGATCGT  >  W3110S.gb/798156‑798291                                                                     |                                                                    gCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTgcgc                                                                   >  1:2168728/1‑71 (MQ=255)    tCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCt                                                                 >  1:201393/1‑70 (MQ=255)     ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                                >  1:3393097/1‑70 (MQ=255)     ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                                >  1:868460/1‑70 (MQ=255)     ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                                >  1:2910058/1‑70 (MQ=255)     ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                                >  1:878473/1‑70 (MQ=255)        tcATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAAcc                                                            <  1:1155338/70‑1 (MQ=255)           tACTGATATCATTAAGGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCAt                                                         <  1:955069/71‑1 (MQ=255)             cTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCAtt                                                        <  1:1849945/70‑1 (MQ=255)             cTGATATCATTAAAGTTATAGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCAtt                                                        >  1:4288724/1‑70 (MQ=255)              tGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTc                                                      >  1:836290/1‑71 (MQ=255)              tGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTc                                                      >  1:1627517/1‑71 (MQ=255)                atatCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTc                                                      <  1:825417/69‑1 (MQ=255)                   tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                                >  1:3024482/1‑56 (MQ=255)                   tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                                  <  1:2636200/70‑1 (MQ=255)                   tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                                  <  1:663640/70‑1 (MQ=255)                      ttAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCaa                                                              <  1:3185034/55‑1 (MQ=255)                      ttAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                                  <  1:3146047/67‑1 (MQ=255)                      ttAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACg                                              <  1:3288730/71‑1 (MQ=255)                       tAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTgcg                                                                    >  1:1965209/1‑48 (MQ=255)                       tAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTgcg                                                                    >  1:2251775/1‑48 (MQ=255)                       tAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGt                                             >  1:1018639/1‑71 (MQ=255)                       tAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGt                                             >  1:3859420/1‑71 (MQ=255)                        aaaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAAc                                               >  1:2208905/1‑68 (MQ=255)                        aaaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGt                                             <  1:2929733/70‑1 (MQ=255)                         aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                            <  1:3858246/70‑1 (MQ=255)                         aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                            <  1:3853013/70‑1 (MQ=255)                         aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                            >  1:1441232/1‑70 (MQ=255)                         aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                            >  1:3157147/1‑70 (MQ=255)                           gTTATCGCCGAATGCCCCGACGTTTTCCATCGACCAACCTTGCGCCTCa                                                               >  1:460953/1‑49 (MQ=255)                             tATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCaa                                                 <  1:922558/61‑1 (MQ=255)                             tATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCaa                                                 <  1:595861/61‑1 (MQ=255)                              aTCGCCGAACGCCACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCaaa                                                >  1:4308478/1‑61 (MQ=255)                              aTCGCCGAACGCCACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCaaa                                                >  1:4297423/1‑61 (MQ=255)                              aTCGCCGAACGCCACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCaaa                                                >  1:562976/1‑61 (MQ=255)                               tCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCttt                                     <  1:2507155/71‑1 (MQ=255)                               tCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCttt                                     <  1:519702/71‑1 (MQ=255)                               tCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCtt                                      <  1:2231185/70‑1 (MQ=255)                                 gCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGc                                   <  1:4562591/71‑1 (MQ=255)                                   cGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCt                                  <  1:3981835/70‑1 (MQ=255)                                    gAATGCCACGACGTTTTCCATCGACCAACCTTGCGcc                                                                  >  1:937163/1‑37 (MQ=255)                                    gAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGt                                >  1:3899107/1‑71 (MQ=255)                                      aTGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                                  >  1:934284/1‑51 (MQ=255)                                       tGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgc                             >  1:2477087/1‑71 (MQ=255)                                       tGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgc                             >  1:3562860/1‑71 (MQ=255)                                         ccACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCAAACg                                              <  1:2928949/52‑1 (MQ=255)                                         ccACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCAAACg                                              <  1:3526026/52‑1 (MQ=255)                                         ccACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgccg                           <  1:772919/71‑1 (MQ=255)                                           acgacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgccgcc                         <  1:3768822/71‑1 (MQ=255)                                            cgacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGccgccg                        <  1:576022/71‑1 (MQ=255)                                            cgacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGccgccg                        <  1:775800/71‑1 (MQ=255)                                            cgacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGccgccg                        <  1:665252/71‑1 (MQ=255)                                             gacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCcgccgc                       <  1:2560191/71‑1 (MQ=255)                                               cgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGt                     >  1:1768401/1‑71 (MQ=255)                                                gTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTAc                                         >  1:1173069/1‑50 (MQ=255)                                                gTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTg                    <  1:4476455/71‑1 (MQ=255)                                                 ttttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                   <  1:3757720/71‑1 (MQ=255)                                                 ttttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                   <  1:1306414/71‑1 (MQ=255)                                                 ttttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                   <  1:3031884/71‑1 (MQ=255)                                                  tttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                                  <  1:3220205/39‑1 (MQ=255)                                                   ttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCa                  >  1:759774/1‑70 (MQ=255)                                                    tCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCaa                 <  1:194501/70‑1 (MQ=255)                                                    tCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCaa                 <  1:2648057/70‑1 (MQ=255)                                                      cATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAata               <  1:2207775/70‑1 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                            >  1:2288898/1‑39 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgccgcc                         <  1:209991/58‑1 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCTGCGTGCAATATc             >  1:111663/1‑70 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                   <  1:4483246/64‑1 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATc             >  1:2193925/1‑70 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATc             >  1:364518/1‑70 (MQ=255)                                                        tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCa            <  1:1491966/71‑1 (MQ=255)                                                         cGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCaa           >  1:131516/1‑71 (MQ=255)                                                          gACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCttt                                     <  1:4503209/44‑1 (MQ=255)                                                          gACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGtt                               >  1:2396084/1‑50 (MQ=255)                                                          gACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCaa           <  1:1387245/70‑1 (MQ=255)                                                          gACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAAc          <  1:4564150/71‑1 (MQ=255)                                                          gACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAAc          >  1:4276673/1‑71 (MQ=255)                                                             cAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCt        >  1:721419/1‑70 (MQ=255)                                                             cAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCt        >  1:2199104/1‑70 (MQ=255)                                                             cAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTg       <  1:4489422/71‑1 (MQ=255)                                                               aCCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGa      >  1:4074254/1‑70 (MQ=255)                                                               aCCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGAt     >  1:800219/1‑71 (MQ=255)                                                                ccTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGt                                <  1:1830219/43‑1 (MQ=255)                                                                 cTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGAt     >  1:4085998/1‑69 (MQ=255)                                                                 cTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGAt     >  1:3384143/1‑69 (MQ=255)                                                                    gcgcCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGATCGt  <  1:267242/69‑1 (MQ=255)                                                                     cgcCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGACGCCGCGTGCAATATCAACCTGATCGt  <  1:1175343/68‑1 (MQ=255)                                                                     |                                                                    GCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGATCGT  >  W3110S.gb/798156‑798291

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