Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	3,316,538	T→G	Q330P (CAG→CCG)	nusA ←	transcription termination/antitermination L factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	3,316,538	0	T	G	100.0%	142.6 / NA	42	Q330P (CAG→CCG)	nusA	transcription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/16);  total (26/16)

AGCCTTCTTCTACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAA  >  NC_000913_3_hsa_tpiA/3316395‑3316680                                                                                                                                                |                                                                                                                                               aGCCATCATCTACCAGAACAGGCGCGAAGTCTTCGTCGATGTCGAGATACTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGaa                                                                                                                                           <  1:286013/149‑1 (MQ=255)   ccTTCTTCTACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGc                                                                                                                                         >  1:510352/1‑149 (MQ=255)    cttcttCTACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATGGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGcc                                                                                                                                        >  1:148841/1‑149 (MQ=255)      tcttctACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAg                                                                                                                                      <  1:416712/149‑1 (MQ=255)       cttctACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGa                                                                                                                                     >  2:273109/1‑149 (MQ=255)        ttctACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGAc                                                                                                                                    >  2:32986/1‑149 (MQ=255)            aCCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCAc                                                                                                                                >  1:233218/1‑149 (MQ=255)               aGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGtt                                                                                                                             >  1:487536/1‑149 (MQ=255)               aGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCATGATGCTTAGCCTGCAGGTCGTCAACGGTCAGCACTGTGAGGTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGtt                                                                                                                             >  2:699435/1‑149 (MQ=255)                   cAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGa                                                                                                                         <  1:273109/149‑1 (MQ=255)                     gTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACCCACGTTCTGAcc                                                                                                                       >  2:482842/1‑149 (MQ=255)                         cgAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTa                                                                                                                   >  2:497565/1‑149 (MQ=255)                                            gagcTATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTGCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCACATCGCCTGCGCCa                                                                                                >  2:648430/1‑149 (MQ=255)                                            gagaTATTTGGTGAAGGTGTCGATCGCTACGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCCGCCGCACGATCTGACAGGTACCGCCAATCGCGCGCGccg                                                                                                >  2:248137/1‑148 (MQ=255)                                             agaTATTTGGTGAAGGTGTCGATCGCTGCGTACGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCCACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAACAAGACGCACGTTCTGACCGTTACGGCCAATAGCCTgcgcacg                                                                                               >  2:415283/1‑146 (MQ=255)                                                  tttGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTAc                                                                                          >  1:337997/1‑149 (MQ=255)                                                   ttGGTGATGGTTTTGGTCGCTGCGTGCGCTTCCGCCTGATGCTCAGCCTGCAGGTCGTCAACCGTCATGACGTTGAGTTCCCCACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTAcc                                                                                         <  2:194992/149‑1 (MQ=255)                                                      gTGAAGGGGTCGATCGCTGCGTGCGCTTCCGTCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGc                                                                                      <  2:165592/149‑1 (MQ=255)                                                            gtgtCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTCCCGGCTTCAAc                                                                                >  2:574226/1‑149 (MQ=255)                                                               tCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGc                                                                             <  2:423210/149‑1 (MQ=255)                                                                    cGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGatat                                                                        <  2:487536/149‑1 (MQ=255)                                                                          gtgcgCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATgg                                                                  >  1:421418/1‑149 (MQ=255)                                                                                gtCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGtt                                                            <  2:99113/148‑1 (MQ=255)                                                                                     ccTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCt                                                       <  1:441905/149‑1 (MQ=255)                                                                                       tGATTCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTc                                                     <  2:413089/149‑1 (MQ=255)                                                                                         aTGCTTAGCCCACATGGCGTCAAAGGTCATCACGTTGAGTTCCGAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGAGATCCATGGTGTGTTTATCTTCAt                                                   <  2:86212/149‑1 (MQ=255)                                                                                           gCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATcc                                                 >  2:573192/1‑149 (MQ=255)                                                                                             ttAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATccac                                               <  2:720842/149‑1 (MQ=255)                                                                                                  cTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGa                                          <  1:580351/149‑1 (MQ=255)                                                                                                       ggTCGCCAACGGTCATCACGTTGAGTTCGCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCTAATCGCTTGCGCCAGGTTTCCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGaa                                     <  2:324047/149‑1 (MQ=255)                                                                                                        gtcgtcAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAg                                    <  1:291208/149‑1 (MQ=255)                                                                                                         tcgtcAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAgcaca                                   >  2:183340/1‑145 (MQ=255)                                                                                                                  gTCATGCCGTTGAGTTCCCAACCGCGCAGCGGCGAAGCCAGACGCACGTTCGGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGc                          <  2:570229/149‑1 (MQ=255)                                                                                                                        aCGTTGAGTTCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACGACGATAGAAGCAACGTCTGCCGGTGc                    >  1:92530/1‑149 (MQ=255)                                                                                                                                ttCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCAACCCCGCTAGAATCAGCGTCTGCCCGTGCCATTGCGt            >  2:13425/1‑149 (MQ=255)                                                                                                                                ttCCCAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGAACATTGTTGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGt            >  1:537178/1‑149 (MQ=255)                                                                                                                                   ccAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGAACGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATCGAAGCAACGACTGCCGGTGCCATTGCGCTaa         >  1:54705/1‑149 (MQ=255)                                                                                                                                    cAACCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAAt        >  1:567802/1‑149 (MQ=255)                                                                                                                                     aaCCGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATc       >  2:550540/1‑149 (MQ=255)                                                                                                                                       ccGCTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTTATCAc     >  2:422812/1‑149 (MQ=255)                                                                                                                                          cTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGaa  >  1:619571/1‑149 (MQ=255)                                                                                                                                          cTCAGCGGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGaa  >  2:515652/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                               AGCCTTCTTCTACCAGAACAGTCGCGAAGTCTTCGTCGATGTCGAGATATTTGGTGAAGGTGTCGATCGCTGCGTGCGCTTCCGCCTGATGCTTAGCCTGCAGGTCGTCAACGGTCATCACGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAA  >  NC_000913_3_hsa_tpiA/3316395‑3316680

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