Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,355,216	T→C	14.9%	Y328C (TAT→TGT)	yihE ←	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,355,216	0	T	C	14.9%	118.5 / 13.7	47	Y328C (TAT→TGT)	yihE	predicted kinase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (5/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCT  >  Exported/2355157‑2355282                                                            |                                                                   taaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:281605/1‑63 (MQ=255)   aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                         <  1:3389168/69‑1 (MQ=255)         aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:761408/71‑1 (MQ=255)         aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:3595646/71‑1 (MQ=255)         aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:2958280/71‑1 (MQ=255)          aCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                <  1:3377020/71‑1 (MQ=255)               ccagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                              >  1:4310701/1‑52 (MQ=255)               ccagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                              >  1:3298566/1‑52 (MQ=255)                cagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                          >  1:1212992/1‑71 (MQ=255)                      cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     >  1:2396406/1‑70 (MQ=255)                        gccaAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTa                                                               <  1:675536/42‑1 (MQ=38)                        gacaAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTa                                                               <  1:965228/40‑1 (MQ=38)                         ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:3989157/48‑1 (MQ=255)                         ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:3849032/48‑1 (MQ=255)                         ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:2919393/48‑1 (MQ=255)                         ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                  <  1:3343559/54‑1 (MQ=255)                          caAATCTTTTTCATGATCTACTCTCTTCGATTAACACCTatg                                                             >  1:3189067/1‑40 (MQ=37)                          caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:3149841/1‑38 (MQ=255)                          caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:1133073/1‑38 (MQ=255)                          caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:1194804/1‑38 (MQ=255)                          caAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTatg                                                             >  1:2375038/1‑40 (MQ=38)                          caAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTatg                                                             >  1:774990/1‑40 (MQ=38)                             aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:2633648/71‑1 (MQ=255)                             aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:909682/71‑1 (MQ=255)                             aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:630926/71‑1 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                         <  1:2778344/42‑1 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                >  1:3124640/1‑51 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                >  1:416352/1‑51 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:2658563/1‑71 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:814380/1‑71 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:4343580/1‑71 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:1784110/1‑71 (MQ=255)                              tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:1791323/70‑1 (MQ=255)                                tttttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAAtt                                                     >  1:3828928/1‑44 (MQ=255)                                tttttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAAtt                                                     >  1:1063275/1‑44 (MQ=255)                                  tttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     <  1:119827/58‑1 (MQ=255)                                   ttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     <  1:4044288/57‑1 (MQ=255)                                   ttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTGGAACTTTTGCCt                       >  1:3605216/1‑71 (MQ=255)                                      aTGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGtt                    <  1:3476139/71‑1 (MQ=255)                                              ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCtata               <  1:2268380/68‑1 (MQ=255)                                              ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCta                 <  1:4363735/66‑1 (MQ=255)                                              ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCta                 <  1:250299/66‑1 (MQ=255)                                                   tcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGc       >  1:3949315/1‑71 (MQ=255)                                                     cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGc       <  1:96321/69‑1 (MQ=255)                                                     cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg     >  1:1336467/1‑71 (MQ=255)                                                        ttAACACCTATGTATTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCt  >  1:2014098/1‑71 (MQ=255)                                                        ttAACACCTATGTATTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTAGCCGTCt  >  1:3388958/1‑71 (MQ=255)                                                            |                                                                   TAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCT  >  Exported/2355157‑2355282

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