Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,588,402	C→T	15.3%	L1032L (CTG→CTA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,588,402	0	C	T	15.3%	191.8 / 22.0	72	L1032L (CTG→CTA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (6/5);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCAC  >  W3110S.gb/4588339‑4588470                                                                |                                                                     gagaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGt                                                               >  1:309341/1‑71 (MQ=255) gagaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGt                                                               >  1:4008729/1‑71 (MQ=255)  agaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTc                                                              <  1:4951357/71‑1 (MQ=255)   gaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCa                                                             <  1:2587935/71‑1 (MQ=255)    aTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTc                                                              <  1:2530333/69‑1 (MQ=255)     tCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCa                                                             >  1:1111600/1‑69 (MQ=255)     tCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCa                                                             >  1:148136/1‑69 (MQ=255)       gcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAg                                                          <  1:1346026/70‑1 (MQ=255)        cgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGc                                                        <  1:3916547/71‑1 (MQ=255)             gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGACCGTTGCACCAGCGAGTCAAAGGCTTc                                                     >  1:1766984/1‑69 (MQ=255)             gCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGa                                                   <  1:2268006/71‑1 (MQ=255)                aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAgg                                                 <  1:317862/70‑1 (MQ=255)                aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAgg                                                 >  1:837434/1‑70 (MQ=255)                   aataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGa                                                   >  1:1330690/1‑65 (MQ=255)                   aataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAAt                                             >  1:464771/1‑71 (MQ=255)                    ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                                         >  1:1980342/1‑58 (MQ=255)                    ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAgg                                                         >  1:664147/1‑58 (MQ=255)                    ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGaaa                                              <  1:3199109/69‑1 (MQ=255)                     taaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCaa                                                            <  1:3823066/54‑1 (MQ=255)                       aCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCt                                          <  1:1566190/70‑1 (MQ=255)                         tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTgcg                                       >  1:2337787/1‑71 (MQ=255)                          gcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGa                                                                 >  1:3734576/1‑44 (MQ=255)                          gcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGa                                                                 >  1:2366881/1‑44 (MQ=255)                           cctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTgcgg                                      <  1:3562083/70‑1 (MQ=255)                               cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtgt                                 >  1:1450462/1‑71 (MQ=255)                               cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtgt                                 >  1:225311/1‑71 (MQ=255)                               cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtgt                                 >  1:1456319/1‑71 (MQ=255)                                 aTGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTG‑GGCGTGTCGt                              <  1:3138190/71‑1 (MQ=255)                                  tGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                              <  1:3879234/71‑1 (MQ=255)                                   gCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                              >  1:3481759/1‑70 (MQ=255)                                   gCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                              >  1:328881/1‑70 (MQ=255)                                   gCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                              >  1:4963810/1‑70 (MQ=255)                                    ccGGTTGCGCGTTCGGTGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:3101021/70‑1 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtg                                  >  1:2783308/1‑65 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAAt                                             >  1:2993377/1‑54 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTaa                            >  1:193041/1‑71 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:718920/70‑1 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:4184242/70‑1 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:202644/70‑1 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:3431988/70‑1 (MQ=255)                                    ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                             <  1:2364317/70‑1 (MQ=255)                                      ggTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATc                          <  1:1659172/71‑1 (MQ=255)                                      ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtgt                                 <  1:2804200/64‑1 (MQ=255)                                      ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCgtgt                                 <  1:3560865/64‑1 (MQ=255)                                      ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATc                          <  1:3030764/71‑1 (MQ=255)                                      ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATc                          <  1:119765/71‑1 (MQ=255)                                           cgcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAgg                                                 <  1:4358861/43‑1 (MQ=255)                                             cgTTCGGGGAATGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                    >  1:2071341/1‑70 (MQ=255)                                             cgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                    >  1:2349603/1‑70 (MQ=255)                                               ttCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAgg                                                 <  1:3637620/39‑1 (MQ=255)                                                tCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATcc                                           >  1:3469220/1‑44 (MQ=255)                                                tCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATcc                                           >  1:3297349/1‑44 (MQ=255)                                                tCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATcc                                           >  1:1012320/1‑44 (MQ=255)                                                tCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGc                >  1:2179823/1‑71 (MQ=255)                                                 cGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATc                          >  1:4113825/1‑60 (MQ=255)                                                  ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGa              >  1:3016752/1‑71 (MQ=255)                                                  ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGa              >  1:2565122/1‑71 (MQ=255)                                                     gAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAttt           <  1:4221113/71‑1 (MQ=255)                                                     gAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAttt           <  1:3248735/71‑1 (MQ=255)                                                      aaCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAtttt          <  1:2010787/71‑1 (MQ=255)                                                      aaCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAtttt          <  1:65243/71‑1 (MQ=255)                                                       aCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTa         <  1:2484989/71‑1 (MQ=255)                                                       aCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTa         <  1:973817/71‑1 (MQ=255)                                                        cGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTac        <  1:2910972/71‑1 (MQ=255)                                                        cGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTac        <  1:576144/71‑1 (MQ=255)                                                           tGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTg                                         >  1:4603832/1‑35 (MQ=255)                                                            gCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgc                                        <  1:5076291/35‑1 (MQ=255)                                                            gCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                              <  1:561411/45‑1 (MQ=255)                                                            gCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTaa                            >  1:4805564/1‑47 (MQ=255)                                                            gCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTaa                            >  1:823795/1‑47 (MQ=255)                                                            gCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTaa                            >  1:2045689/1‑47 (MQ=255)                                                              aCCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTcac  >  1:1420765/1‑71 (MQ=255)                                                                cAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                    <  1:729822/51‑1 (MQ=255)                                                                |                                                                     GAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCAC  >  W3110S.gb/4588339‑4588470

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