Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,394,072	A→G	23.5%	D97G (GAT→GGT)	yfbP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,394,072	0	A	G	23.5%	150.1 / 49.6	81	D97G (GAT→GGT)	yfbP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base G (7/12);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACG  >  W3110S.gb/2394005‑2394140                                                                    |                                                                     aagaagTTCGAAGAATGGCTGGATGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    >  1:3620103/1‑70 (MQ=255) aagaagTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    <  1:1077922/70‑1 (MQ=255) aagaagTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    <  1:878512/70‑1 (MQ=255) aagaagTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    >  1:871009/1‑70 (MQ=255) aagaagTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    >  1:1856810/1‑70 (MQ=255)   gaagTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATc                                                                    >  1:944978/1‑68 (MQ=255)              aaTGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTgcagca                                                      <  1:1634470/71‑1 (MQ=255)               aTGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTgcagca                                                      <  1:2101216/70‑1 (MQ=255)               aTGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCaa                                                     <  1:3126511/71‑1 (MQ=255)               aTGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCaa                                                     <  1:2542809/71‑1 (MQ=255)                tGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTgcagca                                                      <  1:3904838/69‑1 (MQ=255)                  gCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAAtc                                                  >  1:1810977/1‑71 (MQ=255)                  gCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAAt                                                   >  1:2878468/1‑70 (MQ=255)                         gTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCaa                                                                  <  1:3958891/48‑1 (MQ=255)                         gTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCaaa                                                    <  1:3348583/62‑1 (MQ=255)                            aaTGATGATGATCCCCATTAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCg                                           <  1:1677337/68‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCgtg                                         <  1:3095965/70‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCgtg                                         <  1:3712758/70‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCgtg                                         <  1:2877594/70‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCg                                           <  1:3675040/68‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCg                                           <  1:2211451/68‑1 (MQ=255)                            aaTGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCg                                           <  1:3108271/68‑1 (MQ=255)                              tgatgatgaTCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTg                                      >  1:3792272/1‑71 (MQ=255)                                atgatgatCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCt                                    >  1:3997812/1‑71 (MQ=255)                                atgatgatCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCt                                    >  1:3984751/1‑71 (MQ=255)                                 tgatgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTc                                   >  1:4067662/1‑71 (MQ=255)                                  gatgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTgcagca                                                      <  1:1529170/51‑1 (MQ=255)                                  gatgatCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCg                                  >  1:3959993/1‑71 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGGaaa                                                    >  1:219779/1‑51 (MQ=38)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaaa                                                    >  1:1038711/1‑51 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaaa                                                    >  1:42056/1‑51 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaaa                                                    >  1:1461520/1‑51 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaaa                                                    >  1:1112464/1‑51 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGa                                <  1:165001/71‑1 (MQ=255)                                    tgatCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGa                                <  1:662177/71‑1 (MQ=255)                                      atCCCCATGAAAATGGCATCAGCATCGAAGATCaa                                                                  <  1:1331366/35‑1 (MQ=255)                                       tCCCCATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGAAAATCTCTTCCGTGTTGCTCGGAAAt                             <  1:603127/71‑1 (MQ=255)                                         cccATGAAAATGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaa                                                     <  1:1462465/45‑1 (MQ=255)                                         cccATGAAAATGGCATCAGCATCGAAGGGCAAGACTTTGCAGCaa                                                     <  1:530144/45‑1 (MQ=38)                                          ccATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTc                          <  1:1846233/71‑1 (MQ=255)                                          ccATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTc                          <  1:1547798/71‑1 (MQ=255)                                           cATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCaa                                                     <  1:291993/43‑1 (MQ=255)                                             tGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCttt                       >  1:3304653/1‑71 (MQ=255)                                             tGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCttt                       >  1:3363513/1‑71 (MQ=255)                                             tGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCttt                       >  1:3581690/1‑71 (MQ=255)                                             tGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCttt                       >  1:2982248/1‑71 (MQ=255)                                               aaaaTGGCATCAGCATCGAAGGTCAAGACTTTGCAGCaaa                                                    >  1:3488888/1‑40 (MQ=255)                                               aaaaTGGCATCAGCATCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGa                                <  1:968245/60‑1 (MQ=255)                                                     gCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACggg                <  1:452533/70‑1 (MQ=255)                                                     gCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACggg                <  1:2294726/70‑1 (MQ=255)                                                      cATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACggg                <  1:2244165/69‑1 (MQ=255)                                                        tCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGaaa                              >  1:1051093/1‑53 (MQ=255)                                                        tCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGaaa                              >  1:1118794/1‑53 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCgt                                          >  1:3830570/1‑39 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGc           >  1:701456/1‑70 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCt          >  1:2145127/1‑71 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCt          >  1:1754929/1‑71 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCt          >  1:4034684/1‑71 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCt          >  1:1032735/1‑71 (MQ=255)                                                          aGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCt          >  1:2571675/1‑71 (MQ=255)                                                           gCATCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGGAGGCtt         <  1:1826246/71‑1 (MQ=255)                                                           gCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCtt         <  1:2553678/71‑1 (MQ=255)                                                            cATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtg        <  1:951647/71‑1 (MQ=255)                                                            cATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtg        <  1:1248805/71‑1 (MQ=255)                                                             aTCGAATATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGTAAATGTCTTTCAACGTGAAGGCtt         <  1:3064918/69‑1 (MQ=255)                                                              tCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTc                      <  1:2478671/55‑1 (MQ=255)                                                              tCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTc                      <  1:1862507/55‑1 (MQ=255)                                                              tCGAAGGTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTc                      <  1:482848/55‑1 (MQ=255)                                                              tCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCtt         <  1:2185919/68‑1 (MQ=255)                                                              tCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCtt         <  1:4026308/68‑1 (MQ=255)                                                              tCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgat      <  1:3004897/71‑1 (MQ=255)                                                              tCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgat      <  1:2027573/71‑1 (MQ=255)                                                              tCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgat      <  1:364625/71‑1 (MQ=255)                                                                gAAGATCAAGACTTTGCAGCAAATCTCTTCCgtgt                                        >  1:3752216/1‑35 (MQ=255)                                                                gAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCtt                        >  1:463270/1‑51 (MQ=255)                                                                gAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACggg                <  1:1271319/59‑1 (MQ=255)                                                                gAAGATCAAGACTTT‑‑AGCAAATCTCTTCCGTGTTGCTCGGAAATGTCtt                        >  1:1570266/1‑49 (MQ=255)                                                                  aGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACg  <  1:1767076/71‑1 (MQ=255)                                                                  aGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACg  >  1:287760/1‑71 (MQ=255)                                                                  aGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACg  >  1:2053992/1‑71 (MQ=255)                                                                  aGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACg  <  1:950122/71‑1 (MQ=255)                                                                    aTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgat      <  1:2706200/65‑1 (MQ=255)                                                                    aTCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTtgat      <  1:3740048/65‑1 (MQ=255)                                                                    |                                                                     AAGAAGTTCGAAGAATGGCTGGAAGTGAATGATGATGATCCCCATGAAAATGGCATCAGCATCGAAGATCAAGACTTTGCAGCAAATCTCTTCCGTGTTGCTCGGAAATGTCTTTCAACGGGAAGGCTTGATGACG  >  W3110S.gb/2394005‑2394140

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