Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,387,662	T→C	47.4%	Y94H (TAT→CAT)	elaD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,387,662	0	T	C	47.4%	‑1.2 / 53.3	38	Y94H (TAT→CAT)	elaD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/11);  new base C (6/12);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATTGTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAGATGATC  >  W3110S.gb/2387593‑2387719                                                                      |                                                          atTGTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtata                                                          >  1:2287220/1‑71 (MQ=255) atTGTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtata                                                          >  1:1541348/1‑71 (MQ=255)  ttGTTAGAACCCCAACAAATTAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtata                                                          >  1:1466224/1‑70 (MQ=255)  ttGTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtata                                                          >  1:127059/1‑70 (MQ=255)  ttGTTAGAACCCCAACAAATGAGAACCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtata                                                          >  1:553513/1‑70 (MQ=255)    gTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtatata                                                        <  1:2705776/70‑1 (MQ=255)    gTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATaa                                                       >  1:1277353/1‑71 (MQ=255)    gTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATaa                                                       >  1:1652085/1‑71 (MQ=255)     ttAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATtatata                                                        <  1:1111157/69‑1 (MQ=255)       aGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCt                                                    >  1:1927712/1‑71 (MQ=255)        gAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCt                                                    <  1:2629398/70‑1 (MQ=255)        gAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCt                                                    <  1:1799669/70‑1 (MQ=255)             ccAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGcaca                                              <  1:835265/71‑1 (MQ=255)             ccAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGcaca                                              <  1:1638237/71‑1 (MQ=255)              caacaaATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACat                                             <  1:1347046/71‑1 (MQ=255)              caacaaATGAGAACCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGcaca                                              >  1:1606940/1‑70 (MQ=255)                  aaaTGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACata                                            >  1:536375/1‑68 (MQ=255)                  aaaTGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACata                                            >  1:2729555/1‑68 (MQ=255)                   aaTGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGc                                        <  1:2479142/71‑1 (MQ=255)                       agaCCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGcaacc                                    <  1:231881/71‑1 (MQ=255)                           ccTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACt                                <  1:1924256/71‑1 (MQ=255)                            cTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGcaacc                                    <  1:1598163/66‑1 (MQ=255)                            cTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAAAtatata                                                        <  1:360843/46‑1 (MQ=255)                              aTCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATc                                          <  1:2287109/58‑1 (MQ=255)                              aTCGATTATCCGACTTATCAAAAATCCTTGTTTAAAATTATATAATCTCGCACATATCGcaacc                                    <  1:919524/64‑1 (MQ=255)                               tCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGc                                        <  1:1398336/59‑1 (MQ=255)                               tCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGcaacc                                    <  1:1407145/63‑1 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACata                                            >  1:1605463/1‑53 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAAt                                                      >  1:2646461/1‑43 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATc                                                     >  1:658564/1‑44 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATc                                                     >  1:2322660/1‑44 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATc                                                     >  1:1669051/1‑44 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTcccccc                          >  1:1478920/1‑71 (MQ=255)                                 gATTATCCGACTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTccccc                           >  1:1173350/1‑70 (MQ=255)                                      tCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCaa                      <  1:512472/70‑1 (MQ=255)                                          aCTTATCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATca                  <  1:192104/70‑1 (MQ=255)                                            ttATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGcaacc                                    <  1:599320/50‑1 (MQ=255)                                              aTCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACAt               <  1:2018180/69‑1 (MQ=255)                                               tCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATca                  <  1:2804234/65‑1 (MQ=255)                                               tCAAAAATGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATca                  <  1:142302/65‑1 (MQ=255)                                                  aaaaTGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAAt                    <  1:102359/60‑1 (MQ=255)                                                   aaaTGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAg        <  1:1493144/71‑1 (MQ=255)                                                     aTGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAgat      >  1:2653737/1‑71 (MQ=255)                                                      tGCTTGTTTAAAATTACATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAgatg     <  1:325416/71‑1 (MQ=255)                                                        cTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAgatgat   >  1:1150796/1‑71 (MQ=255)                                                           gTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAGATGAtc  <  1:2285812/69‑1 (MQ=255)                                                                      |                                                          ATTGTTAGAACCCCAACAAATGAGACCCTATCGATTATCCGACTTATCAAAAATGCTTGTTTAAAATTATATAATCTGGCACATATCGCAACCAACTCCCCCCTCAAATCACATGATTCAGATGATC  >  W3110S.gb/2387593‑2387719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑