Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,123,875	G→A	43.9%	R255H (CGC→CAC)	damX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,123,875	0	G	A	43.9%	30.3 / 82.7	66	R255H (CGC→CAC)	damX	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (22/7);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGA  >  W3110S.gb/4123809‑4123940                                                                   |                                                                  ttCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                <  1:4289794/70‑1 (MQ=255) ttCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                <  1:470367/70‑1 (MQ=255) ttCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                <  1:1493211/70‑1 (MQ=255)       atggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCAccc                                                                <  1:573232/64‑1 (MQ=255)       atggcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCAccc                                                                <  1:818771/64‑1 (MQ=255)          gcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTcag                                                      <  1:4753481/71‑1 (MQ=255)          gcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcag                                                      >  1:1177052/1‑71 (MQ=255)          gcaatgCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcag                                                      >  1:759418/1‑71 (MQ=255)               gCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAgg                                                  >  1:2641234/1‑70 (MQ=255)                cATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCg                                                >  1:1579417/1‑71 (MQ=255)                  tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                              <  1:2590045/71‑1 (MQ=255)                  tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                              <  1:2464725/71‑1 (MQ=255)                  tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                              <  1:134036/71‑1 (MQ=255)                  tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGt                                              <  1:412502/71‑1 (MQ=255)                  tCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGt                                              <  1:1538831/71‑1 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                         >  1:819633/1‑58 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGt                                                         >  1:4745280/1‑58 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:3445972/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:3065629/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:4488677/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:2208098/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:2071034/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:1540128/1‑71 (MQ=255)                    gcgTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGa                                            >  1:1194027/1‑71 (MQ=255)                       tGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                         <  1:585846/71‑1 (MQ=255)                       tGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                         <  1:123033/71‑1 (MQ=255)                        gACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGAtt                                          >  1:4012189/1‑69 (MQ=255)                        gACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTg                                         >  1:1352547/1‑70 (MQ=255)                         acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                        <  1:3176995/70‑1 (MQ=255)                             ggCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcagca                                                    >  1:3303516/1‑54 (MQ=255)                              gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                  <  1:3377345/71‑1 (MQ=255)                              gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                  <  1:4739911/71‑1 (MQ=255)                              gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGa                                   >  1:3458715/1‑70 (MQ=255)                                gAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaa                                <  1:682640/71‑1 (MQ=255)                                 aaaaCGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGt                                                         >  1:4787091/1‑45 (MQ=255)                                 aaaaCGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGt                                                         >  1:1530209/1‑45 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                >  1:1796780/1‑36 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                >  1:4430077/1‑36 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                >  1:3499526/1‑36 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCGccc                                                                >  1:136106/1‑36 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAc                             >  1:2658802/1‑71 (MQ=255)                                   aaCGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAc                             >  1:3186819/1‑71 (MQ=255)                                   aaCGCAAACCGCTAAACGTCCGTCCACTACGCGccc                                                                >  1:496584/1‑36 (MQ=255)                                     cGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCg                           <  1:3433788/71‑1 (MQ=255)                                      gCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgc                          <  1:878710/71‑1 (MQ=255)                                       cAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                              >  1:2735705/1‑66 (MQ=255)                                          aCCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                       >  1:2724159/1‑70 (MQ=255)                                          aCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcagca                                                    >  1:2200860/1‑41 (MQ=255)                                          aCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcagca                                                    >  1:3895394/1‑41 (MQ=255)                                          aCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcagca                                                    >  1:3460666/1‑41 (MQ=255)                                           ccGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagc                      >  1:3527216/1‑70 (MQ=255)                                            cGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagcaa                    <  1:2379679/71‑1 (MQ=255)                                               tGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaa                        >  1:4059266/1‑64 (MQ=255)                                                 aaCGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTg               <  1:2281317/71‑1 (MQ=255)                                                 aaCGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTg               <  1:3185737/71‑1 (MQ=255)                                                  aCGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTg               <  1:2735887/70‑1 (MQ=255)                                                   cgtccgtcCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaa             <  1:594691/71‑1 (MQ=255)                                                      ccgtccACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGa              >  1:2218942/1‑67 (MQ=255)                                                       cgtccACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg         <  1:2056586/71‑1 (MQ=255)                                                       cgtccACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg         <  1:199239/71‑1 (MQ=255)                                                         tccACTACGCGCCCAGCTCGTCAGCAGGCGGTGAt                                           >  1:350401/1‑35 (MQ=255)                                                         tccACTACGCGCCCAGCTCGTCAGCAGGCGGTGAt                                           >  1:1348373/1‑35 (MQ=255)                                                         tccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACg         >  1:3536930/1‑69 (MQ=255)                                                               aCGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGcc    >  1:1247116/1‑68 (MQ=255)                                                               aCGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaa                        >  1:74563/1‑48 (MQ=255)                                                               aCGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAAc          >  1:3752038/1‑62 (MQ=255)                                                                cgcaCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGa  <  1:3969918/69‑1 (MQ=255)                                                                   |                                                                  TTCGCAATGGCAATGCATCGCGTGACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCGA  >  W3110S.gb/4123809‑4123940

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