Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	709,001	A→G	52.7%	Y102H (TAT→CAT)	yceI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	709,001	0	A	G	52.7%	15.8 / 128.2	91	Y102H (TAT→CAT)	yceI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/31);  new base G (20/28);  total (32/59)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGA  >  Exported/708934‑709066                                                                    |                                                                  tCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                <  1:2594353/71‑1 (MQ=255)  cGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGatat                                                                  <  1:3594096/68‑1 (MQ=255)  cGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGttt                                                                >  1:1927218/1‑70 (MQ=255)   gCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:1060986/70‑1 (MQ=255)   gCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:411735/70‑1 (MQ=255)   gCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:1617600/70‑1 (MQ=255)   gCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:139822/70‑1 (MQ=255)    ccGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGc                                                             <  1:1634278/71‑1 (MQ=255)      ggTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCgg                                                           <  1:919089/71‑1 (MQ=255)        tAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTa                                                         <  1:2283292/71‑1 (MQ=255)        tAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTa                                                         <  1:1918443/71‑1 (MQ=255)        tAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTa                                                         >  1:3413423/1‑71 (MQ=255)         aaTATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAt                                                        >  1:721510/1‑71 (MQ=255)          atatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                       <  1:2376849/71‑1 (MQ=255)           tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                       >  1:1815239/1‑70 (MQ=255)           tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTg                                                      <  1:1212837/71‑1 (MQ=255)           tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTg                                                      <  1:1452016/71‑1 (MQ=255)           tatCAAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTg                                                      <  1:2190381/71‑1 (MQ=255)              ccAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTa                                                         <  1:2375935/65‑1 (MQ=255)                aGTTCGTCAGCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTa                                                         <  1:1176712/63‑1 (MQ=255)                aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                 <  1:323455/71‑1 (MQ=255)                aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                 <  1:1276296/71‑1 (MQ=255)                 gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaaa                                                <  1:944833/71‑1 (MQ=255)                 gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaaa                                                <  1:3727285/71‑1 (MQ=255)                  ttCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAAt                                               <  1:1349773/71‑1 (MQ=255)                   tCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGa                                                     >  1:854511/1‑64 (MQ=255)                   tCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                <  1:662131/53‑1 (MQ=255)                    cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCgg                                                           >  1:74650/1‑57 (MQ=255)                    cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCgg                                                           >  1:1496799/1‑57 (MQ=255)                    cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATc                                              <  1:2360931/70‑1 (MQ=255)                    cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCt                                             >  1:1721518/1‑71 (MQ=255)                     gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                            <  1:3491357/71‑1 (MQ=255)                      tCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGc                                           <  1:2326931/71‑1 (MQ=255)                        aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCAc                                         >  1:3522481/1‑71 (MQ=255)                          cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTg                                       >  1:2215594/1‑71 (MQ=255)                          cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                              >  1:3036201/1‑48 (MQ=255)                          cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                              >  1:3203150/1‑48 (MQ=255)                          cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                              >  1:2985337/1‑48 (MQ=255)                          cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTg                                       <  1:2475020/71‑1 (MQ=255)                           gTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:2975816/46‑1 (MQ=255)                           gTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:749371/46‑1 (MQ=255)                           gTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGtttt                                                               <  1:358679/46‑1 (MQ=255)                           gTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGc                                      <  1:896104/71‑1 (MQ=255)                               ttcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGc                                                             >  1:1760379/1‑44 (MQ=255)                               ttcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGa                                                     >  1:2385108/1‑52 (MQ=255)                               ttcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAg                                                    >  1:601410/1‑53 (MQ=255)                               ttcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAg                                                    >  1:2433830/1‑53 (MQ=255)                               ttcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAg                                                    >  1:2298915/1‑53 (MQ=255)                                tcttcACGCTGGTGGAGGTGAATGTTGCCTTTGGATATTTTGCGGt                                                          <  1:1878995/46‑1 (MQ=255)                                tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGa                                    <  1:3519897/68‑1 (MQ=255)                                tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGa                                    <  1:3016628/68‑1 (MQ=255)                                tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAg                                  <  1:3547182/70‑1 (MQ=255)                                tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAg                                  <  1:1466457/70‑1 (MQ=255)                                 cttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaaa                                                <  1:3564375/55‑1 (MQ=255)                                 cttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAg                                  <  1:2061993/69‑1 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGAt                                <  1:440346/68‑1 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGt                              <  1:1071758/70‑1 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGAt                                <  1:2372950/68‑1 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGAt                                <  1:2644201/68‑1 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGtt                             >  1:1694985/1‑71 (MQ=255)                                    cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGtt                             >  1:3554428/1‑71 (MQ=255)                                      cGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTg                                                              >  1:1091123/1‑36 (MQ=255)                                      cGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGaa                                                 <  1:1674806/49‑1 (MQ=255)                                      cGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGaa                                                 <  1:3339259/49‑1 (MQ=255)                                      cGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTa                           <  1:1495878/71‑1 (MQ=255)                                        cTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATc                         <  1:1699923/71‑1 (MQ=255)                                           gtggAGGTGAATGTTGCCTGTGGATGTTTTGCGGTGTTGAGGAAATCTGCACTGCGAAGATGTTTATcg                        <  1:2651637/69‑1 (MQ=255)                                           gtggAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATcg                        <  1:397961/69‑1 (MQ=255)                                            tggAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGt                     <  1:541507/71‑1 (MQ=255)                                              gAGGTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGtt                             <  1:563293/61‑1 (MQ=255)                                                ggTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                 <  1:326580/71‑1 (MQ=255)                                                ggTGAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                 <  1:1899849/71‑1 (MQ=255)                                                   gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                     <  1:1885734/48‑1 (MQ=255)                                                   gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt              >  1:3135783/1‑71 (MQ=255)                                                   gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt              >  1:2789070/1‑71 (MQ=255)                                                    aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt              <  1:363194/70‑1 (MQ=255)                                                    aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt              <  1:2021972/70‑1 (MQ=255)                                                    aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt              <  1:2975749/70‑1 (MQ=255)                                                      tGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt           >  1:3642065/1‑71 (MQ=255)                                                       gTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt           >  1:2040806/1‑70 (MQ=255)                                                       gTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt           >  1:553372/1‑70 (MQ=255)                                                        ttGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGc                                      <  1:2328984/42‑1 (MQ=255)                                                          gCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGa                                    >  1:1165880/1‑42 (MQ=255)                                                             tgtgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCTCGTTc                   >  1:239903/1‑56 (MQ=255)                                                             tgtgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTc                   >  1:2185429/1‑56 (MQ=255)                                                             tgtgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATc    >  1:2456201/1‑71 (MQ=255)                                                              gtgGATGTTTTGCGGTATTGAGGAAATCTGCACTGc                                      >  1:175187/1‑36 (MQ=255)                                                               tgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATc    <  1:949817/69‑1 (MQ=255)                                                               tgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATc    <  1:339419/69‑1 (MQ=255)                                                               tgGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATc    <  1:3202611/69‑1 (MQ=255)                                                               tgGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTGGGCGTGATTAGTATc    <  1:3017987/69‑1 (MQ=255)                                                               tgGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGa  >  1:1268736/1‑71 (MQ=255)                                                                    |                                                                  TCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGA  >  Exported/708934‑709066

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