Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	2,168,532	G→T	17.1%	intergenic (‑226/‑180)	yegR ← / → yegS	uncharacterized protein/phosphatidylglycerol kinase, metal‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	2,168,532	0	G	T	17.1%	100.3 / 5.4	45	intergenic (‑226/‑180)	yegR/yegS	uncharacterized protein/phosphatidylglycerol kinase, metal‑dependent
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base T (6/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTA  >  NC_000913_3_hsa_tpiA/2168387‑2168679                                                                                                                                                  |                                                                                                                                                    aTTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAATATGCCACGATTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                                                                                                                                  >  2:100684/1‑149 (MQ=255) aTTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                                                                                                                                  <  2:312919/149‑1 (MQ=255)   tataATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGTTATATGCTCTTTTTTGAAAtt                                                                                                                                                >  2:407922/1‑149 (MQ=255)        taaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCg                                                                                                                                           <  1:252734/149‑1 (MQ=255)           aTGTTCTGTACCACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTTTATTGAAGATGCCAATGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTg                                                                                                                                        <  1:640619/149‑1 (MQ=255)           aTGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTg                                                                                                                                        <  1:520311/149‑1 (MQ=255)            tattCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCATTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTgg                                                                                                                                       <  2:418216/147‑1 (MQ=255)              ttCTGTACAACATTTCCTACATAAGTAGGCAACAGTGTCATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCa                                                                                                                                     >  2:521264/1‑149 (MQ=255)                   tACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAAATCCATGGGGGCTTTAATATATTATATTGGAGATGACACTGTTTACTTGCATCTTAGGTATATGTGCTTTTTTACAATTTATCGCTGGCCGCCGt                                                                                                                                >  2:648899/1‑149 (MQ=255)                     caacaTTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTa                                                                                                                              >  2:217771/1‑149 (MQ=255)                                     aGTTGGAATTACGGACATTCAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATACTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTCTCGGTGGCAGCCGTTAGTATTCGCTGTCCCCa                                                                                                              <  1:691102/149‑1 (MQ=255)                                       tAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTTTTTTTTTTCTGTGGCAGCCGTTAGTATTCGCTGTCCCCAtt                                                                                                            >  2:225623/1‑149 (MQ=255)                                       tAGGAACTATGGACAATGTGGCCCTTCAGGTTAATTCCATGGGGGCTTTAATATATCATATTGAAGATGTCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCAtt                                                                                                            <  2:129724/149‑1 (MQ=255)                                           acttgcGGACATTGAGGTCTTTCAGTGTAACTCTATGGGTGTTTTAATATATTATATTGAAGGTGCCACTGTTTAGTTGAATATTAGGTCTATGTATTTTTTTTAAATTTGTCGGTGGCAGGCGTTAGTATTCGCTGTCCCCATTGCaa                                                                                                        <  2:184773/144‑1 (MQ=255)                                           aaTTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGGTGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTGTCGGTGGCAGCCGTTAGTAATCACTGTCCCCAtagtaa                                                                                                        >  2:474198/1‑144 (MQ=255)                                                 ggACATTGAGGCCCTTCAGGGTAACTCCATGGGGCCTTTAATAAATTCTCTTGAAGCTCCCAATCATTAAATAAATAATAGGTATATGCCCTTTTTTTATATTTATCTGTGCCAGCCGCTAGAATTCTCCTTCCCCCTTACACGCCCCt                                                                                                  >  2:246008/1‑149 (MQ=255)                                                     aTTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCtggtg                                                                                              <  1:521264/149‑1 (MQ=255)                                                     aTTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCtggtg                                                                                              <  1:488935/149‑1 (MQ=255)                                                         aGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTaa                                                                                          <  1:591067/149‑1 (MQ=255)                                                           gCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTTTTTTTTTTCTGTGTCAGCAGTTAGTTTTCGCTGTCCCCATTGCAAGGTCCTGGTGGTaatc                                                                                        >  1:60246/1‑147 (MQ=255)                                                             ccTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCAc                                                                                      <  1:394786/149‑1 (MQ=255)                                                              cTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTTAATTTTTTTTTTGGCATCCGTTATTATATACTGTCCCCCTTGCAAGCTCCTGGTGTGAACCACt                                                                                     >  2:352989/1‑149 (MQ=255)                                                               ttCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTg                                                                                    >  2:600675/1‑149 (MQ=255)                                                                        aCTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCa                                                                           >  2:69751/1‑149 (MQ=255)                                                                          tCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATa                                                                         >  2:396885/1‑149 (MQ=255)                                                                                   gCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTTTTTTTTTTCTGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCCCTGAATCCTCCATACTTGAACTg                                                                >  1:237345/1‑149 (MQ=255)                                                                                           atattatattGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTa                                                        <  1:598870/149‑1 (MQ=255)                                                                                            tattatattGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTTTTTTTTTTTGGTGGCAGCCGTTAGTATTCGCTTTCCCCATTGCAAGCTCCTGGTGGTAACCACTGCTTCCTCCTTTCTTTAACTGTCttttttt                                                       >  2:57949/1‑147 (MQ=255)                                                                                                 tattGATGATGTCATTGTTTAGTTGTATATTAGGTATATGCTCTGTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCGCATTGCTAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTAtcctcc                                                  <  2:235112/149‑1 (MQ=255)                                                                                                    tGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGAc                                               <  1:114853/149‑1 (MQ=255)                                                                                                         aTGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAACTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTGTCCTCCGACTTTCa                                          >  1:639938/1‑149 (MQ=255)                                                                                                           gCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATc                                        <  2:607555/149‑1 (MQ=255)                                                                                                                  tttAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTCAGTCTACACTGCCCCCCTTGCAAGCTCCTGGTGGTAAGCGCTGAAACCTCCATAATTGAAATGAATTTCTATCCTCCGACTTTCATccagctgc                                 >  2:240950/1‑143 (MQ=255)                                                                                                                     tggtGAATGTTAGGTATATGTCCTTTTTTGAATTGTATGGGTGGGCGGGGTTAGTAATCGCTGTGCCCATTGAAAGTACCTGGTGGTAACCACTGAATTCTTCATACTTGAACTGACTTTGTATCCTCCGACTTTCATCCTGTTCTGAc                              <  2:171367/146‑1 (MQ=255)                                                                                                                      gTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATAGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACt                             >  1:552481/1‑149 (MQ=255)                                                                                                                       ttGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTTTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATACTCCGCCTTTCATCCTGTTCTGAcgc                            >  1:333291/1‑147 (MQ=255)                                                                                                                       ttGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTc                            >  1:616916/1‑149 (MQ=255)                                                                                                                       ttGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTc                            >  1:71594/1‑149 (MQ=255)                                                                                                                          aaTATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCAc                         <  1:273696/149‑1 (MQ=255)                                                                                                                              ttAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTTCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCttt                     >  1:554460/1‑149 (MQ=255)                                                                                                                               taagTATATAATCTTATATGAAATTTATCGGTGGCAGCCGATAGTAATCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCtttt                    <  2:652850/146‑1 (MQ=255)                                                                                                                               tAGGTATTGGCTCTTTTTTTACATTCATCGGTCGCAGCCGTTCGTATCCGCTGTCCCCATTGCAAGCTCCTCGTGGTAACCACTGAATCCTGCTTACTTTAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCtttt                    <  1:205491/149‑1 (MQ=255)                                                                                                                                cggTATTTGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTg                   <  1:623403/148‑1 (MQ=255)                                                                                                                                aGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTg                   <  1:76275/149‑1 (MQ=255)                                                                                                                                                 tGAAATTTCTCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAGTCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACGCtc  >  2:389460/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTCCCCATTGCAAGCTCCTGGTGGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTA  >  NC_000913_3_hsa_tpiA/2168387‑2168679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑