Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	798,218	A→G	5.9%	G212G (GGT→GGC)	ybhA ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	798,218	0	A	G	5.9%	229.2 / 5.3	85	G212G (GGT→GGC)	ybhA	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base G (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

CCAGCGGCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGAT  >  W3110S.gb/798150‑798288                                                                     |                                                                       ccAGCGGCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAAcc                                                                      >  1:745454/1‑71 (MQ=255)   aGCGGCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCtt                                                                    >  1:2739220/1‑71 (MQ=255)    gCGGCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTg                                                                   >  1:2975929/1‑71 (MQ=255)       gCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTgcgc                                                                >  1:2168728/1‑71 (MQ=255)          tCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCt                                                              >  1:201393/1‑70 (MQ=255)           ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                             >  1:878473/1‑70 (MQ=255)           ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                             >  1:868460/1‑70 (MQ=255)           ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                             >  1:2910058/1‑70 (MQ=255)           ccAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                             >  1:3393097/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:1627517/1‑71 (MQ=255)                    tGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTc                                                   >  1:836290/1‑71 (MQ=255)                      atatCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTc                                                   <  1:825417/69‑1 (MQ=255)                         tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTc                                                             >  1:3024482/1‑56 (MQ=255)                         tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                               <  1:663640/70‑1 (MQ=255)                         tCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCa                                               <  1:2636200/70‑1 (MQ=255)                           aTTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAAc                                                                       >  1:2136568/1‑44 (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:2251775/1‑48 (MQ=255)                             tAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTgcg                                                                 >  1:1965209/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)                               aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAAc                                                                       >  1:2701127/1‑40 (MQ=255)                               aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                         <  1:3853013/70‑1 (MQ=255)                               aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                         <  1:3858246/70‑1 (MQ=255)                               aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                         >  1:3157147/1‑70 (MQ=255)                               aaGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGtt                                         >  1:1441232/1‑70 (MQ=255)                                 gTTATCGCCGAATGCCCCGACGTTTTCCATCGACCAACCTTGCGCCTCa                                                            >  1:460953/1‑49 (MQ=255)                                 gTTATCGCCGAATGCCACGACGTTTTCCATCGACCaa                                                                        <  1:2056132/37‑1 (MQ=255)                                   tATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCaa                                              <  1:595861/61‑1 (MQ=255)                                   tATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCaa                                              <  1:922558/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:519702/71‑1 (MQ=255)                                     tCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCttt                                  <  1:2507155/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:3562860/1‑71 (MQ=255)                                             tGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTgc                          >  1:2477087/1‑71 (MQ=255)                                               ccACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCAAACg                                           <  1:3526026/52‑1 (MQ=255)                                               ccACGACGTTTTCCATCGACCAGCCTTGTGCCTCAACCCATTTCGTCAAACg                                           <  1:2928949/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:665252/71‑1 (MQ=255)                                                  cgacgTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGccgccg                     <  1:775800/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:1306414/71‑1 (MQ=255)                                                       ttttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                <  1:3031884/71‑1 (MQ=255)                                                       ttttCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGc                <  1:3757720/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:364518/1‑70 (MQ=255)                                                              tCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATc          >  1:2193925/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:2199104/1‑70 (MQ=255)                                                                   cAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCt     >  1:721419/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)                                                                     |                                                                       CCAGCGGCTTCCAGCATACTGATATCATTAAAGTTATCGCCGAATGCCACGACGTTTTCCATCGACCAACCTTGCGCCTCAACCCATTTCGTCAAACGTTTACCTTTGCTGTTGCCGCCGCGTGCAATATCAACCTGAT  >  W3110S.gb/798150‑798288

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