Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	285,318	C→A	intergenic (+150/‑77)	yagF → / → yagG	CP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	285,318	0	C	A	98.1%	174.1 / NA	52	intergenic (+150/‑77)	yagF/yagG	CP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (25/26);  minor base G (0/1);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGCGGCTTCGT  >  NC_000913_3_hsa_tpiA/285176‑285459                                                                                                                                               |                                                                                                                                              aaGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTaa                                                                                                                                        >  1:776465/1‑150 (MQ=255)   gagaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAAtt                                                                                                                                      <  4:82853/150‑1 (MQ=255)    agaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTg                                                                                                                                     >  4:262472/1‑150 (MQ=255)      aTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGaa                                                                                                                                   >  1:37907/1‑150 (MQ=255)        aaaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGGTATATCGTGGGCGTCACAAAAGCAAAAGAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATa                                                                                                                                 >  1:645044/1‑150 (MQ=255)             cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATATCGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCa                                                                                                                            <  2:776465/150‑1 (MQ=255)             cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCa                                                                                                                            <  4:447251/150‑1 (MQ=255)                        tatTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCAt                                                                                                                 <  4:205015/150‑1 (MQ=255)                           ttCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCAATTaa                                                                                                              >  3:458893/1‑150 (MQ=255)                               cGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACg                                                                                                          >  3:354597/1‑150 (MQ=255)                                     ggtcTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAGCACGAtgttg                                                                                                    <  2:133389/147‑1 (MQ=255)                                      ggCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGAtgttgt                                                                                                   >  3:84088/1‑150 (MQ=255)                                          ttgttgTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAAtt                                                                                               >  1:783620/1‑150 (MQ=255)                                             ttgttTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCgg                                                                                            <  2:37907/150‑1 (MQ=255)                                                  tttGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCacac                                                                                       <  4:81337/150‑1 (MQ=255)                                                  tttGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGAGATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCacac                                                                                       >  3:105966/1‑150 (MQ=255)                                                  tttGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGAGTCACAAAAGCAAAATAACATAACTAATAAACAATAAGTTACCACTGATGGTTAATTTAATAGCTCATTACACTGCATTAACACGATGTTGTAATTCGGcaaat                                                                                       >  1:728796/1‑147 (MQ=255)                                                      aaaTTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCt                                                                                                                              <  2:684947/107‑1 (MQ=255)                                                      aaaTTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCt                                                                                                                              >  1:684947/1‑107 (MQ=255)                                                       aaTTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACAt                                                                                  <  4:534161/150‑1 (MQ=255)                                                            attaATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAgg                                                                              <  2:133645/149‑1 (MQ=255)                                                            attaATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAgg                                                                              >  1:133645/1‑149 (MQ=255)                                                              taATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgtg                                                                           <  2:97677/150‑1 (MQ=255)                                                              taATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgtg                                                                           >  3:269922/1‑150 (MQ=255)                                                                 tttAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTaa                                                                        <  4:172331/150‑1 (MQ=255)                                                                       cAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtga                                                                  <  2:158113/150‑1 (MQ=255)                                                                             gTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGc                                                            >  3:329712/1‑150 (MQ=255)                                                                               tGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCaa                                                          <  2:238372/150‑1 (MQ=255)                                                                                  gaTATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTa                                                       <  2:829420/150‑1 (MQ=255)                                                                                      tatCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCa                                                   >  1:743169/1‑150 (MQ=255)                                                                                      tatCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATAGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCa                                                   >  3:438582/1‑150 (MQ=255)                                                                                              gCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTc                                                                                                                             >  4:492038/1‑68 (MQ=255)                                                                                              gCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTc                                                                                                                             <  3:492038/68‑1 (MQ=255)                                                                                                    caAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAAtt                                     <  2:273157/150‑1 (MQ=255)                                                                                                     aaaaGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACg                                                             >  2:795128/1‑125 (MQ=255)                                                                                                     aaaaGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACg                                                             <  1:795128/125‑1 (MQ=255)                                                                                                             aaTAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAAtt                                                                                               <  4:194847/83‑1 (MQ=255)                                                                                                             aaTAACGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAAtt                                                                                               >  3:194847/1‑83 (MQ=255)                                                                                                                aaCGTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTTATGACGCAAATAACCATGAAAGACAAAATTGGCTACGGGCTg                         <  4:269922/150‑1 (MQ=255)                                                                                                                   gTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAac                                                                                                             <  2:16510/63‑1 (MQ=255)                                                                                                                   gTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAac                                                                                                             >  1:16510/1‑63 (MQ=255)                                                                                                                   gTAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGga                      <  4:354597/150‑1 (MQ=255)                                                                                                                    tAATTCGGAAATAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGgag                     >  2:43510/1‑150 (MQ=255)                                                                                                                               tAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGca                                                                                          <  3:415163/70‑1 (MQ=255)                                                                                                                               tAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGca                                                                                          >  4:415163/1‑70 (MQ=255)                                                                                                                               tAAGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGc          >  4:389643/1‑150 (MQ=255)                                                                                                                                 aGATATGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGCgg        <  4:329712/150‑1 (MQ=255)                                                                                                                                     ttgctcATTGGTTGTTAATTGAATAGCTGATTTCAGTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGTTGTCACTCTTATGACGCAGTTAACCATGAAAGACAAAATCGGCTACGGGCTGGGATACACCTCCTGCGGCTTc    <  2:728796/145‑1 (MQ=255)                                                                                                                                     atGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAgg                                                                              <  3:354940/76‑1 (MQ=255)                                                                                                                                     atGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAgg                                                                              >  4:354940/1‑76 (MQ=255)                                                                                                                                     atGACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGCGGCTTc    <  4:120054/150‑1 (MQ=255)                                                                                                                                       gACCATTGATGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGCGGCTTCgc  >  2:443036/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                              AAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTAACCATGAAAGACAAAATTGGCTACGGGCTGGGAGACACCGCCTGCGGCTTCGT  >  NC_000913_3_hsa_tpiA/285176‑285459

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