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%	122.5 / NA	38	L179L (CTG→TTG)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/20);  total (18/20)

AAGATGGCC‑GCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTG  >  NC_000913_3_hsa_tpiA/4111522‑4111799                                                                                                                                              |                                                                                                                                          aaGATGGCCTGCCCCGCTGGGCTGCAGCTCTGTGAAACAGTATAGAGATTCATCAGTACAAAGGCTTTGCTTCTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAAGAGCGAAGGTCTACCCCttt                                                                                                                                    <  2:607978/149‑1 (MQ=255)            ccccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTc                                                                                                                         <  2:311859/149‑1 (MQ=255)                    ggCAGCAGGTCTGTGAAACATTATAGAGATTCATCGCCACAAAGGCTTTGCTTTTTAACATCGAGACTAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                 <  1:168413/149‑1 (MQ=255)                                        gTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCATACCTTCAAGCGATTTAGATAGCGCCAGCTTAATCGGTGCAACAGCGAAGGTCAATCCCTTTTCGCCGTTGTCCGCGACAACATGATGCTGTGCACCTTCt                                                                                             <  2:651926/149‑1 (MQ=255)                                          atgggggttCAGCGGCACAAAGGCTGTGCTTGTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACCGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGt                                                                                           <  1:678432/142‑1 (MQ=255)                                            agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtct                                                                                         >  2:132830/1‑149 (MQ=255)                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg                                                                                      >  2:624420/1‑148 (MQ=255)                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg                                                                                      <  1:624420/148‑1 (MQ=255)                                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                     <  1:122181/149‑1 (MQ=255)                                                 ttCAACGGGACAAAGGCTTTGCTTTTTGTAATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAACGTCAATCCCTTCTCTGCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgtg                                                                                    <  2:3670/149‑1 (MQ=255)                                                    aTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaa                                                                                 >  2:646812/1‑149 (MQ=255)                                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTTTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaat                                                                                >  2:100034/1‑149 (MQ=255)                                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                    >  2:649085/1‑113 (MQ=255)                                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                    <  1:649085/113‑1 (MQ=255)                                                       ggCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGATGtaataa                                                                              >  1:77936/1‑149 (MQ=255)                                                         cacaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACg                                                                            <  1:256357/149‑1 (MQ=255)                                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                        >  2:660753/1‑149 (MQ=255)                                                                               aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                      >  1:483609/1‑149 (MQ=255)                                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGAAAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCGTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagtc                                                     >  1:120842/1‑148 (MQ=255)                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTg                                              >  1:177125/1‑149 (MQ=255)                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTg                                              >  2:570979/1‑149 (MQ=255)                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                              >  2:463845/1‑149 (MQ=255)                                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCCTTAGTAGCGACTg                                              >  2:328598/1‑149 (MQ=255)                                                                                         aCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                            <  1:132830/149‑1 (MQ=255)                                                                                            ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                         <  1:463845/149‑1 (MQ=255)                                                                                                aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                     >  1:463116/1‑149 (MQ=255)                                                                                                 gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                    <  1:646812/149‑1 (MQ=255)                                                                                                  cGATTCCGATAGCGTCAGTCTAAGCGGTTCAGCCGCTAATGTCAACCACTTTTTGCCGTTGTCCGCGACAACATAACGCGGTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                   <  2:709099/149‑1 (MQ=255)                                                                                                  cGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                   <  2:481646/149‑1 (MQ=255)                                                                                                  cGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                   >  1:481646/1‑149 (MQ=255)                                                                                                        agctcgcgcCAGCTTAATCGGTTCAACCGCGAAGGTCAACCCCTTTTCGCCGTTGTCTGCGACAACATAACGCAGTGCACCTTCTGTCTCTGTGTAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                             <  1:113197/144‑1 (MQ=255)                                                                                                          agaGCGCCAGCTTAACCGGTTCAACAGCCAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                           <  1:83562/147‑1 (MQ=255)                                                                                                          aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                           <  1:660753/149‑1 (MQ=255)                                                                                                          aTAGCGCCAGCTGAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                           <  1:71124/149‑1 (MQ=255)                                                                                                                                aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa     <  1:656542/149‑1 (MQ=255)                                                                                                                                 aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg    >  2:440509/1‑149 (MQ=255)                                                                                                                                  cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACCTAATACAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGCgt   >  1:642992/1‑149 (MQ=255)                                                                                                                                   aGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCCCGGTGTAATACCGTTTGTTTTTCCCCGCCGTTATTAGAGCATCCAGTCTTTTCTGTCTTTGTGTTGTGGTCATTAATTGGGCGAgtc  >  2:22536/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                          AAGATGGCC‑GCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTG  >  NC_000913_3_hsa_tpiA/4111522‑4111799

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