Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,623,057	C→T	15.5%	H49H (CAC→CAT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,623,057	0	C	T	15.5%	220.2 / 26.5	84	H49H (CAC→CAT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (42/29);  new base T (9/4);  total (51/33)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGA  >  W3110S.gb/4622987‑4623122                                                                       |                                                                  tCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCAc                                                                   <  1:3522491/71‑1 (MQ=255) tCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCAc                                                                   <  1:3664535/71‑1 (MQ=255) tCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCAc                                                                   <  1:1484533/71‑1 (MQ=255)   aaCGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACg                                                                  <  1:2495913/70‑1 (MQ=255)        tATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGa                                                            <  1:2623659/71‑1 (MQ=255)          ttCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGat                                                                <  1:1358872/65‑1 (MQ=255)          ttCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGAc                                                           <  1:1884325/70‑1 (MQ=255)           tCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCAc                                                                   <  1:4507104/61‑1 (MQ=255)            cgcgACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACaa                                                         >  1:3682118/1‑70 (MQ=255)                  aacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACg                                                                  <  1:130986/55‑1 (MQ=255)                    cacTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTGAGc                                                >  1:4475555/1‑71 (MQ=255)                    cacTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTACATGATATGACAATGCCTGAGc                                                >  1:781934/1‑71 (MQ=255)                    cacTATCTACGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGc                                                <  1:4240343/71‑1 (MQ=255)                     acTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCg                                               >  1:2335292/1‑71 (MQ=255)                      cTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGa                                                                 <  1:1157556/52‑1 (MQ=255)                      cTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCgt                                              >  1:964773/1‑71 (MQ=255)                      cTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCg                                               >  1:2898904/1‑70 (MQ=255)                       tATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGa                                                                 <  1:2989100/51‑1 (MQ=255)                       tATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGa                                                                 <  1:3188435/51‑1 (MQ=255)                        aTCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTGAGCgt                                              <  1:2513934/69‑1 (MQ=255)                        aTCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCgtgt                                            >  1:490815/1‑71 (MQ=255)                          ctcCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtct                                          >  1:1351406/1‑71 (MQ=255)                           tcCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtctc                                         >  1:4497739/1‑71 (MQ=255)                           tcCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtctc                                         >  1:4465655/1‑71 (MQ=255)                            ccGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGatat                                                              >  1:121850/1‑49 (MQ=255)                             cGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCg                                        <  1:1879684/70‑1 (MQ=255)                             cGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGAc                                                           >  1:1977768/1‑51 (MQ=255)                             cGAAGGGCAGATTGCCGCCCTCGCGATGACCAATTTCTTCCACGATATGAc                                                           >  1:4086172/1‑51 (MQ=255)                              taaGGGCAGATTGCCTCCCTCGCGATGACCATTTTCTTCCACg                                                                  <  1:2076664/42‑1 (MQ=39)                                    cAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACaa                                                         >  1:1421952/1‑46 (MQ=255)                                     aGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGa                                                                 >  1:378534/1‑37 (MQ=255)                                      gATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATgg                              >  1:2922521/1‑71 (MQ=255)                                       aTTGCCGCCCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                             >  1:3576177/1‑71 (MQ=255)                                       aTTGCCGCCCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                             <  1:1089233/71‑1 (MQ=255)                                       aTTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtc                                           <  1:1160701/57‑1 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACaa                                                         >  1:2341858/1‑41 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACaa                                                         >  1:1567275/1‑41 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:1697582/1‑71 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:3387063/1‑71 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:4554132/1‑71 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:3921208/1‑71 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:2499816/1‑71 (MQ=255)                                         tGCCGCCCTCGCGATGACCATTTTCTTCCACCATATGACAATGCCTGAGCgtgt                                            >  1:4243126/1‑54 (MQ=255)                                          gccgccCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTTAGCGTGTCTCGCTGACCATGGCGa                           <  1:2656223/70‑1 (MQ=255)                                            cgccCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGAcc                                  <  1:3466264/61‑1 (MQ=255)                                             gccCTCGCGATGACCATTTTCTTCCATGATATGACAATGCCTg                                                   >  1:4540500/1‑43 (MQ=255)                                              cccTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCa                                 >  1:444665/1‑60 (MQ=255)                                              cccTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCga                      <  1:3910543/71‑1 (MQ=255)                                                cTCGCGATGACCATTTTCTTCCACGATATGACAATGc                                                      <  1:1702894/37‑1 (MQ=255)                                                cTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGc                        >  1:3482249/1‑67 (MQ=255)                                                cTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCgaga                    <  1:2862870/71‑1 (MQ=255)                                                cTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCgaga                    <  1:951881/71‑1 (MQ=255)                                                  cgcgATGACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAt                   >  1:1113814/1‑70 (MQ=255)                                                  cgcgATGACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAt                   >  1:3385335/1‑70 (MQ=255)                                                  cgcgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAtt                  >  1:1551289/1‑71 (MQ=255)                                                  cgcgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAtt                  >  1:1640585/1‑71 (MQ=255)                                                   gcgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAtt                  <  1:2547398/70‑1 (MQ=255)                                                   gcgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAtt                  >  1:2627390/1‑70 (MQ=255)                                                   gcgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGAtt                  >  1:2765317/1‑70 (MQ=255)                                                    cgATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCgaga                    <  1:151364/67‑1 (MQ=255)                                                      aTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCa                                 >  1:3670620/1‑52 (MQ=255)                                                      aTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAgg              <  1:959047/71‑1 (MQ=255)                                                        gACCATTTTCTTCCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:3052413/1‑56 (MQ=255)                                                         aCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATgg                              >  1:4357497/1‑52 (MQ=255)                                                         aCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATgg                              >  1:4228218/1‑52 (MQ=255)                                                         aCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGaa            <  1:1595637/70‑1 (MQ=255)                                                          ccATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAAcc          >  1:1180808/1‑71 (MQ=255)                                                          ccATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAAcc          >  1:2617947/1‑71 (MQ=255)                                                            aTTTTCTTCCACGATATGACAATGCCTGAGCgtgt                                            <  1:2748156/35‑1 (MQ=255)                                                            aTTTTCTTCCACGATATGACAATGCCTGAGCGTGtctc                                         <  1:4503292/38‑1 (MQ=255)                                                             ttttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                             <  1:3030555/49‑1 (MQ=255)                                                             ttttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGtt       >  1:2662595/1‑71 (MQ=255)                                                             ttttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGt        <  1:3228421/70‑1 (MQ=255)                                                              tttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCAt                                >  1:2648271/1‑45 (MQ=255)                                                              tttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCAt                                >  1:1524533/1‑45 (MQ=255)                                                              tttCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTc                 <  1:2091414/60‑1 (MQ=255)                                                               ttcttcCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtct     >  1:3213801/1‑71 (MQ=255)                                                               ttcttcCATGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtct     >  1:3111865/1‑71 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCa                                 >  1:3943141/1‑43 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:2077932/1‑49 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGa                           >  1:137570/1‑49 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAgg              >  1:1811751/1‑62 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAg               >  1:1618243/1‑61 (MQ=255)                                                               ttcttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtc      <  1:4030112/70‑1 (MQ=255)                                                                 cttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTTATGAACCGTtctc    <  1:823539/70‑1 (MQ=255)                                                                 cttcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtctc    <  1:4255578/70‑1 (MQ=255)                                                                   tcCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGa             >  1:3328191/1‑59 (MQ=255)                                                                    ccACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGa  <  1:980005/69‑1 (MQ=255)                                                                       |                                                                  TCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGA  >  W3110S.gb/4622987‑4623122

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