Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,419,670	G→T	intergenic (‑214/+95)	ydaG ← / ← racR	Rac prophage; uncharacterized protein/Rac prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,419,670	0	G	T	100.0%	116.7 / NA	34	intergenic (‑214/+95)	ydaG/racR	Rac prophage; uncharacterized protein/Rac prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/18);  total (16/18)

TGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAG  >  NC_000913_3_hsa_tpiA/1419526‑1419814                                                                                                                                                 |                                                                                                                                                 tGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAAtttttt                                                                                                                                              <  1:80323/149‑1 (MQ=255)              agggCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTaa                                                                                                                                 <  2:276760/148‑1 (MQ=255)                 gctggctACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATa                                                                                                                               <  1:181772/148‑1 (MQ=255)                 gctggctACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATa                                                                                                                               >  2:181772/1‑148 (MQ=255)                              ccccTTTCCACTGCCTGGTAACTTGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGAcag                                                                                                                 <  2:174458/149‑1 (MQ=255)                               cccTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGAcagc                                                                                                                >  1:3129/1‑149 (MQ=255)                                  ttACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGa                                                                                                             <  2:354911/149‑1 (MQ=255)                                    aCCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGaaa                                                                                                           <  2:346063/149‑1 (MQ=255)                                            cTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggg                                                                                                   >  2:297038/1‑149 (MQ=255)                                             tGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATgggg                                                                                                  >  1:510711/1‑149 (MQ=255)                                              ggTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATggggg                                                                                                 >  1:535979/1‑149 (MQ=255)                                                         gTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAAATATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAAc                                                                                      <  2:293202/149‑1 (MQ=255)                                                                 cGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGAGATGGGGGGGTTAAATAACGGAATCCa                                                                              >  1:404744/1‑149 (MQ=255)                                                                 cGGCGTTCTGGGGGGCGGGGTGGGTGGTATGCTGGAACTATGGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTGATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCa                                                                              <  1:84593/149‑1 (MQ=255)                                                                  ggCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAAGGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAg                                                                             <  2:410766/149‑1 (MQ=255)                                                                    cGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGa                                                                           <  1:250054/149‑1 (MQ=255)                                                                        cTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGttt                                                                       >  1:46027/1‑149 (MQ=255)                                                                          gaggggCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGTAGAAAGTGATGGGGGGGGTAAATAACGGAATCCAGGAGTTTTc                                                                     >  2:167412/3‑149 (MQ=255)                                                                             ggggggggTGGGTGGTATGGCTGAACGATAGGAAATGTCTAATTGATTGTCAAGAGGCTATGCCTAAGTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGt                                                                  <  2:198433/149‑1 (MQ=255)                                                                               gCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCa                                                                <  1:332750/149‑1 (MQ=255)                                                                                cGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGCCAGCAGACAGTGATGGGGCGGCTAAATAACGGAATCCAGCAGTTCTCCGTCAg                                                               >  2:347947/1‑149 (MQ=255)                                                                                  gggtgggtGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGAc                                                             <  1:159231/149‑1 (MQ=255)                                                                                     tgggtggTATGCTGGAAATATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCat                                                          <  2:42461/149‑1 (MQ=255)                                                                                         tggtATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAAc                                                                                      >  1:174502/1‑117 (MQ=255)                                                                                         tggtATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAAc                                                                                      <  2:174502/117‑1 (MQ=255)                                                                                         tggtATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATAtaa                                                      >  1:568746/1‑149 (MQ=255)                                                                                           gtATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATAtaagt                                                    >  1:434495/1‑149 (MQ=255)                                                                                                     attataGGTAATGCCGAATTGATTGTCAATAGGCTATGGCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCa                                          <  2:419936/147‑1 (MQ=255)                                                                                                            gTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTtg                                   >  2:476288/1‑149 (MQ=255)                                                                                                              aaTGCCTAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACTTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTtgtg                                 <  1:539799/149‑1 (MQ=255)                                                                                                                    tAATTGATTGTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTc                           >  2:338432/1‑149 (MQ=255)                                                                                                                        tgattgTCAATAGGCTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGc                       <  2:568746/149‑1 (MQ=255)                                                                                                                                      cTATGCCTAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCGGAt         >  1:268700/1‑149 (MQ=255)                                                                                                                                            cTAATTTTTTGAGAGTAACCTAATCGGTGATGTAGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAaga   <  1:415797/149‑1 (MQ=255)                                                                                                                                             tAATTTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAagag  >  1:502801/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 TGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTCAATAGGCTATGCCTAATGTTTTGAGCGTAACCTAATAGGTGATGGCGACAGCAGAAAGTGATGGGGGGGTTAAATAACGGAATCCAGGAGTTTTCCGTCAGACCATATAAGTTTAAGTTCCAGTTTTTGTGATGTTCTGGCTTTTCCGTTCAGATTCAAGAG  >  NC_000913_3_hsa_tpiA/1419526‑1419814

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