Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,111,662	G→A	L179L (CTG→TTG)	yiiQ ←	DUF1454 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,111,662	0	G	A	100.0%	113.3 / NA	36	L179L (CTG→TTG)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

AATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGA  >  NC_000913_3_hsa_tpiA/4111514‑4111805                                                                                                                                                     |                                                                                                                                                aaTAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCTGCACAAAGGGTTTGGTTTGTGTCATTTATTCAACCCTTCAAGCGATGTGGAAAGCGGCAGCTTAATCGGTTCAACAGCGAGGGTAAg                                                                                                                                                 >  1:96522/1‑149 (MQ=255)    gaaggaaggGGGGCGCCCGGGGGGGCGGGGGGTCGGTGAAACAGTATAGAGATTCATCGGCACCAAGGCTTCGCTTTTTCTCGTGTATTGAGACCTGCTAGCGTTTCAGGTAGCTCCAGCGTAATCGGTTCAACAGCGAAGGTCAAccc                                                                                                                                              <  2:16859/148‑1 (MQ=255)        gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTCTAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCttt                                                                                                                                          >  1:45841/1‑149 (MQ=255)                          gggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTCGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTGCAACAGCGAAGGTCAACCCCTTTTCTCCTTTGTCCGCGcca                                                                                                                        >  2:358121/1‑149 (MQ=255)                             cagcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaaca                                                                                                                     <  1:342013/149‑1 (MQ=255)                              agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACAt                                                                                                                    >  1:409648/1‑149 (MQ=255)                                  ggTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACg                                                                                                                <  1:341692/149‑1 (MQ=255)                                   gTCTGTGCAACAGTATAGAGATTCATCAGCACCAAGGCATTGCTTTATGTCATTTATTCAAACCTTCCAGCCATTCAGATAGCTTCAGCTTACTCGGTTCAACAGCCACGGTCAACCCCTCTTCGCCGTTGTCCGCCCCAACACAAcac                                                                                                               >  2:16653/1‑147 (MQ=255)                                   gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCCGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGc                                                                                                               >  2:16354/1‑149 (MQ=255)                                    tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                              <  1:357223/149‑1 (MQ=255)                                    tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                              <  1:280465/149‑1 (MQ=255)                                        tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGc                                                                                                          <  1:116177/149‑1 (MQ=255)                                        tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGc                                                                                                          >  1:150303/1‑149 (MQ=255)                                             cAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTACGCGACAACATAACGCTGTGCACCtt                                                                                                     >  1:395604/1‑149 (MQ=255)                                                 ataGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCGGCGACAACATAACGCAGTGCACCTTCTGt                                                                                                 >  1:171806/1‑149 (MQ=255)                                                    gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCATCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                              <  2:320234/149‑1 (MQ=255)                                                    gagaGTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                              <  1:398413/149‑1 (MQ=255)                                                       aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTCGgt                                                                                           <  1:187758/149‑1 (MQ=255)                                                            tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtatt                                                                                      >  1:278095/1‑147 (MQ=255)                                                                 acaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGGTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTTAACCCATTTTCGACGTTGTCCGCGACGCCATCACGCAGTGCACCTTCGGTCTCGGTGTAATGACgt                                                                                 >  1:290926/1‑149 (MQ=255)                                                                   aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCTTTTTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgttt                                                                               >  2:155440/1‑149 (MQ=255)                                                                            gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGAAAACATAACGCAGTGCACCTTCTGTCTCGGTGTAGTAACGTTTGtttttgac                                                                      >  1:217225/1‑146 (MQ=255)                                                                                 ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCAAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                 <  1:226866/149‑1 (MQ=255)                                                                                       tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACGTTTGTTTTTCCCTTCCGttatta                                                           >  2:66362/1‑149 (MQ=255)                                                                                            tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                      <  1:23615/149‑1 (MQ=255)                                                                                              aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTGCTAGCGAcga                                                    >  1:620/1‑147 (MQ=255)                                                                                                       aGCGATTCAGATAGAGCCAGCTTAGTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                           >  1:371938/1‑149 (MQ=255)                                                                                                            ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTCATAACGTTTGTTTTTCCACGCCTTTACTAGCTGCTGCAGTTTTTTCTGGc                                      >  2:50424/1‑149 (MQ=255)                                                                                                               aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                         >  1:281300/1‑111 (MQ=255)                                                                                                               aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                         <  2:281300/111‑1 (MQ=255)                                                                                                                gATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTCATCACGTTTGTTTTTCCCCGCAGTTAGTTGCGACTGCAGTTTTTTCCGTCTTtg                                  >  2:237324/1‑149 (MQ=255)                                                                                                                       ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                           >  1:238410/1‑149 (MQ=255)                                                                                                                         ggggaaaTGGGTTCAACAGCGAAGGTCAACCCGTTTTGGCTGTTGTCCGTGACTAAATACTGCAGTGCACCTTCTGTCTCGGTGTAATATCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                         <  2:413724/144‑1 (MQ=255)                                                                                                                                    ttCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg              <  1:81525/149‑1 (MQ=255)                                                                                                                                         cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                   <  2:141352/75‑1 (MQ=255)                                                                                                                                         cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                   >  1:141352/1‑75 (MQ=255)                                                                                                                                                ggTCAACCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTCTACTAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGa  <  1:100569/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                AATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGA  >  NC_000913_3_hsa_tpiA/4111514‑4111805

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