Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,547,470	C→T	75.3%	S91S (TCG→TCA)	yddL ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,547,470	0	C	T	75.3%	135.8 / 49.9	77	S91S (TCG→TCA)	yddL	predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base T (40/18);  total (49/28)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

ATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAT  >  W3110S.gb/1547402‑1547537                                                                     |                                                                    atTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCcgcgcg                                                                    <  1:3322502/70‑1 (MQ=255)    gtgtAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGc                                                                 <  1:2240990/70‑1 (MQ=255)     tgtAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg                                                               <  1:2320713/71‑1 (MQ=255)     tgtAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg                                                               <  1:1475084/71‑1 (MQ=255)     tgtAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg                                                               <  1:2324975/71‑1 (MQ=255)     tgtAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCg                                                               <  1:3317220/71‑1 (MQ=255)            taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCtt                                                        >  1:2324375/1‑71 (MQ=255)            taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCtt                                                        >  1:878054/1‑71 (MQ=255)            taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCtt                                                        <  1:3318680/71‑1 (MQ=255)            taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCtt                                                        >  1:60052/1‑71 (MQ=255)                gTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATc                                                    <  1:2058238/71‑1 (MQ=255)                 tCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCtt                                                        <  1:780020/66‑1 (MQ=255)                  cAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGt                                                  >  1:408362/1‑71 (MQ=255)                   aGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGGTGGCTTCATCGt                                                  <  1:321264/70‑1 (MQ=255)                       aTCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAg                                             >  1:1908267/1‑71 (MQ=255)                       aTCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAg                                             >  1:662191/1‑71 (MQ=255)                         cttatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGcc                                           >  1:1457711/3‑71 (MQ=255)                         cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:3664162/1‑66 (MQ=255)                         cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:2762425/1‑66 (MQ=255)                         cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:3180321/1‑66 (MQ=255)                         cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:1816325/1‑66 (MQ=255)                         cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGcc                                           >  1:1612533/1‑71 (MQ=255)                             ttGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCtt                                         >  1:95048/1‑69 (MQ=255)                               gtgATTGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCa                                                      >  1:1503994/1‑54 (MQ=39)                               gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGc                                                                 >  1:2585471/1‑43 (MQ=255)                               gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCa                                                      >  1:2185929/1‑54 (MQ=255)                               gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:1559403/1‑60 (MQ=255)                               gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:35595/1‑60 (MQ=255)                               gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:504615/1‑60 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGcc                                                                <  1:2468037/43‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGcc                                                                <  1:952059/43‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGcc                                                                <  1:2063073/43‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGcc                                                                <  1:226602/43‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGcc                                                                <  1:3521530/43‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa                                    <  1:3537199/71‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa                                    <  1:1300220/71‑1 (MQ=255)                                tgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa                                    <  1:1170408/71‑1 (MQ=255)                                tcaTGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTaaaa                                    <  1:1868970/69‑1 (MQ=255)                                   tGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCa                                 <  1:2589957/71‑1 (MQ=255)                                     ggCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               >  1:2688605/1‑71 (MQ=255)                                     ggCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGATGGATTAATAGTTATAGCCTTTAAAAt                                   >  1:2339235/1‑67 (MQ=255)                                      gCACTTAACTCAATCAGAGATTCTTGCCGCGTGAg                                                                  <  1:2503472/35‑1 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               >  1:3708615/1‑65 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               >  1:2755423/1‑65 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               >  1:125511/1‑65 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               >  1:3470822/1‑65 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTcc                           >  1:902081/1‑69 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTcc                           >  1:97743/1‑69 (MQ=255)                                           tAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTAATCGTTATAGCCTTTAAAATCATATTcc                           >  1:3230847/1‑69 (MQ=255)                                            aaCTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:507783/1‑47 (MQ=255)                                            aaCTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:997096/1‑47 (MQ=255)                                            aaCTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTta                                                >  1:3124259/1‑47 (MQ=255)                                                cAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGAc                    >  1:2942416/1‑71 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   >  1:622674/1‑51 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   >  1:517276/1‑51 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   >  1:3143138/1‑51 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   >  1:44473/1‑51 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   >  1:3645850/1‑51 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaac               >  1:3113731/1‑70 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa               >  1:3135152/1‑71 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa               >  1:3099335/1‑71 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa               >  1:3532900/1‑71 (MQ=255)                                                     agagATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa               >  1:3368576/1‑71 (MQ=255)                                                     agagATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCa                         >  1:432196/1‑61 (MQ=255)                                                     agagATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCaaaa               >  1:1590765/1‑71 (MQ=255)                                                        gaTTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAAcc             >  1:431717/1‑70 (MQ=255)                                                          ttCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGt                       <  1:3195612/58‑1 (MQ=255)                                                           tCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAAt          >  1:3386570/1‑70 (MQ=255)                                                            cTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCa        >  1:1107408/1‑71 (MQ=255)                                                            cTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCa        >  1:713045/1‑71 (MQ=255)                                                            cTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCa        >  1:42170/1‑71 (MQ=255)                                                             ttGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCa                                 <  1:2636734/45‑1 (MQ=255)                                                                ccGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:2703450/71‑1 (MQ=255)                                                                ccGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:448441/71‑1 (MQ=255)                                                                ccGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:2412867/71‑1 (MQ=255)                                                                ccGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:2329178/71‑1 (MQ=255)                                                                ccGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:3563631/71‑1 (MQ=255)                                                                ccGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTg    <  1:1544109/71‑1 (MQ=255)                                                                  gcgtgAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAt  <  1:2315384/71‑1 (MQ=255)                                                                  gcgtgAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAt  <  1:1987886/71‑1 (MQ=255)                                                                     |                                                                    ATTGTGTAATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAT  >  W3110S.gb/1547402‑1547537

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