Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,581,079	C→T	7.0%	G47S (GGT→AGT)	ypeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,581,079	0	C	T	7.0%	176.0 / 3.7	57	G47S (GGT→AGT)	ypeB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/39);  new base T (2/2);  total (16/41)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCG  >  Exported/1581011‑1581143                                                                     |                                                                 cTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCAcc                                                                  >  1:3310573/1‑69 (MQ=255) cTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                                                                <  1:814559/71‑1 (MQ=255)    cAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCt                                                             >  1:1241856/1‑71 (MQ=255)     aGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCtt                                                            >  1:3124139/1‑71 (MQ=255)     aGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCt                                                             >  1:1176370/1‑70 (MQ=255)      ggggCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTc                                                           <  1:2029208/71‑1 (MQ=255)             gACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:3112732/70‑1 (MQ=255)              aCCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCTAGCTTCCCTGCCgg                                                   <  1:3249508/71‑1 (MQ=255)              aCCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCtt                                                            <  1:1816490/62‑1 (MQ=255)                cAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGa                                                 <  1:1839227/71‑1 (MQ=255)                cAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGa                                                 <  1:1842592/71‑1 (MQ=255)                 aGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGaa                                                <  1:1645088/71‑1 (MQ=255)                   ttGCGTTAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCg                                                                 <  1:3547312/52‑1 (MQ=255)                   ttGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTg                                                       <  1:2077343/62‑1 (MQ=255)                   ttGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTg                                                       <  1:3718624/62‑1 (MQ=255)                     gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTc                                             <  1:2474887/70‑1 (MQ=255)                       gTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATc                                          >  1:3232657/1‑71 (MQ=255)                       gTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATc                                          >  1:1453463/1‑71 (MQ=255)                          aGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCttt                                       <  1:1108169/71‑1 (MQ=255)                            ccccTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGc                                     >  1:1872775/1‑71 (MQ=255)                                  ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:1118798/70‑1 (MQ=255)                                  ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:2754972/70‑1 (MQ=255)                                    ttttAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAAc                             <  1:3323468/71‑1 (MQ=255)                                    ttttAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAAc                             <  1:3330887/71‑1 (MQ=255)                                      ttAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCgg                                   >  1:92985/1‑63 (MQ=255)                                      ttAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCaa                               <  1:259607/67‑1 (MQ=255)                                      ttAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCaa                               <  1:1480697/67‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:3246332/44‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:2600029/44‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:3122903/44‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:2183611/44‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc                                                     <  1:1190374/44‑1 (MQ=255)                                       tAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACca                           <  1:2454243/70‑1 (MQ=255)                                        aaTCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACca                           >  1:1750902/1‑69 (MQ=255)                                        aaTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACcac                          <  1:3280561/70‑1 (MQ=255)                                        aaTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACcac                          <  1:346752/70‑1 (MQ=255)                                         aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACcaca                         <  1:2246571/70‑1 (MQ=255)                                          tCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCtt                                        <  1:2557688/54‑1 (MQ=255)                                           cagcagCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                        >  1:2023935/1‑69 (MQ=255)                                              cagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                    >  1:3458962/1‑70 (MQ=255)                                                gcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATAc                 <  1:2845194/71‑1 (MQ=255)                                                gcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATAc                 >  1:3261751/1‑71 (MQ=255)                                                 cGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGTAACTCATCTTTGCGGGCAAACCACAGCAGATAct                <  1:708014/71‑1 (MQ=255)                                                 cGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATAct                <  1:728153/71‑1 (MQ=255)                                                 cGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATAct                <  1:1911419/71‑1 (MQ=255)                                                  gTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATActc               >  1:240309/1‑71 (MQ=255)                                                   tGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATActcc              <  1:2317372/71‑1 (MQ=255)                                                   tGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATActcc              <  1:1109653/71‑1 (MQ=255)                                                       ttGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcagc                       >  1:3171281/1‑58 (MQ=255)                                                        tGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCgg         >  1:3718176/1‑71 (MQ=255)                                                          cATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                        >  1:2010542/1‑54 (MQ=255)                                                          cATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGc        <  1:3626102/70‑1 (MQ=255)                                                           aTTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGc        <  1:1074373/69‑1 (MQ=255)                                                           aTTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAg      <  1:2951171/71‑1 (MQ=255)                                                             tAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGt    >  1:383419/1‑71 (MQ=255)                                                               gCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:3205793/41‑1 (MQ=255)                                                               gCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:3534510/41‑1 (MQ=255)                                                               gCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCg  <  1:2994883/71‑1 (MQ=255)                                                                    ccGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:1854675/36‑1 (MQ=255)                                                                    ccGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                <  1:1660091/36‑1 (MQ=255)                                                                    ccGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcagc                       <  1:2402675/45‑1 (MQ=255)                                                                     |                                                                 CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCG  >  Exported/1581011‑1581143

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