Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,111,482	G→T	intergenic (+797/+115)	cdh → / ← yiiQ	CDP‑diacylglycerol phosphotidylhydrolase/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,482	0	G	T	100.0%	86.2 / NA	28	intergenic (+797/+115)	cdh/yiiQ	CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/13);  total (15/13)

ATAGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGGGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCA  >  NC_000913_3_hsa_tpiA/4111346‑4111624                                                                                                                                         |                                                                                                                                               ataGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTaa                                                                                                                                    >  1:221026/1‑149 (MQ=255)   aGGCAGTAGGGGGAGCACAAGCCACCTGGGTGGCGGTCGGCACCTTGGCCGCGTTCAGGGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCCCTTTTAtt                                                                                                                                  <  2:6934/149‑1 (MQ=255)      cAGTAGGGGGAGCACAAACCCCCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATtctc                                                                                                                               <  2:88949/149‑1 (MQ=255)               gAGCACAAACCACCCCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGATCACCTTCCAGCTTCCCCCACCGAATAACGTCCTGGAGTGCGCCATTTTAATTCTCCACGCTtat                                                                                                                      <  1:388351/149‑1 (MQ=255)                         caccTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATaa                                                                                                            >  2:442014/1‑149 (MQ=255)                                      gccggcACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGgccg                                                                                               >  1:131645/1‑149 (MQ=255)                                                tGGCCGCGGTCAGAGTGCCGATGAGCTCCCCCAGACTCGGCTGCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAggg                                                                                     <  2:108146/149‑1 (MQ=255)                                                          cAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATGCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCt                                                                           <  2:221026/149‑1 (MQ=255)                                                                  cGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACa                                                                   <  1:89655/149‑1 (MQ=255)                                                                    aTGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGt                                                                 >  2:483761/1‑149 (MQ=255)                                                                     tGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCCTCTTCCAGTGTCCTCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGta                                                                <  1:112425/149‑1 (MQ=255)                                                                     tGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGta                                                                <  2:539338/149‑1 (MQ=255)                                                                       aGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCAAGATTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGtata                                                              <  1:442014/149‑1 (MQ=255)                                                                            ccccAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGAt                                                         <  1:483761/149‑1 (MQ=255)                                                                            ccccAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAAATGGCCGCCCCGCGGGGCAGCGGGTCTGTGAAACAGTCTagccat                                                         >  1:523707/1‑145 (MQ=255)                                                                              ccAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTc                                                       >  1:319718/1‑149 (MQ=255)                                                                                 gACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATc                                                    >  2:401095/1‑149 (MQ=255)                                                                                    tctGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCTACGAAGAACTTCCTGTAGTGAGCCATTTTAATTCTCCAGGCTGATAAGCGAATAAAGGAAGATGGCGGCCCCTCAGGGCAGCAGGCCTGTGAAACAGAATAGAGATTCATCgga                                                 >  2:304578/1‑148 (MQ=255)                                                                                     ctGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGca                                                >  1:251699/1‑149 (MQ=255)                                                                                                        ttCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGt                             >  1:128913/1‑149 (MQ=255)                                                                                                                     cAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAAc                >  2:580481/1‑149 (MQ=255)                                                                                                                     cAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAAc                >  2:76006/1‑149 (MQ=255)                                                                                                                         gaagaaCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGCCTTTGCTTTTTGTCATTTATTCAAACCTTc            >  2:161752/1‑149 (MQ=255)                                                                                                                         gaagaaCTTCCTGGAGTGCGCCATTTTAACTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTc            <  1:514728/149‑1 (MQ=255)                                                                                                                                  ccTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTCTAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTc   <  2:522818/149‑1 (MQ=255)                                                                                                                                  ccTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTc   <  2:131645/149‑1 (MQ=255)                                                                                                                                   cTGGAGTGCGCCATTTTACTTCTCCACGCTTATAAGCGAATAAAGGAAGATGTCCGCCCCGCAGAGCAACAGGTCTGTGAAACAGTATAGAGATTCATCGGATCAAAGGCTTTGCTGCTTGTCATTTATTCAAACCTTCATACGAttta  >  2:101278/1‑147 (MQ=255)                                                                                                                                   cTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTTCAACCTTCAAGCGATTcc  >  2:391244/1‑148 (MQ=255)                                                                                                                                         |                                                                                                                                               ATAGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGGGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCA  >  NC_000913_3_hsa_tpiA/4111346‑4111624

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