Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 1,158,559 C→T 100% Q231* (CAG→TAG)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,158,5590CT96.9% 209.7 / ‑1.9 64Q231* (CAG→TAG) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (36/26);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGC  >  NC_000913_3_hsa_tpiA/1158417‑1158701
                                                                                                                                              |                                                                                                                                              
tcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGtttt                                                                                                                                         >  3:349519/1‑150 (MQ=255)
 cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTc                                                                                                                                        <  3:315482/150‑1 (MQ=255)
 cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGCTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGGGTGGTCTGCACCACATCGGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGattac                                                                                                                                        >  1:251584/1‑148 (MQ=255)
  aGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTCAGAAGAAAGCGTTGGATCTGACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTAACTTTCCAGATACAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTcc                                                                                                                                       >  2:172911/1‑150 (MQ=255)
         tgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGa                                                                                                                                >  1:43100/1‑150 (MQ=255)
                cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCcg                                                                                                                         <  2:100205/150‑1 (MQ=255)
                cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCcg                                                                                                                         >  1:246516/1‑150 (MQ=255)
                  gAACTCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCcgcg                                                                                                                       <  4:269713/150‑1 (MQ=255)
                    ccccGGTAGTTGCTTTTAGCATTTACGGTGTCATCGAACTTTGCCTGGTCCCGTTTGGTCTGCACCACATCGGGAACGTACCGTTCCAGATGCAGATTGGTGAATACACCATCGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGtc                                                                                                                     <  4:207909/149‑2 (MQ=255)
                      ccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTtat                                                                                                                   >  1:450719/1‑150 (MQ=255)
                       cGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCAAGATGCAGCTTGGTGCATCCCCCCACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTtata                                                                                                                  >  1:529326/1‑150 (MQ=255)
                       cGGTAGTGGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCGGGTACCGTTTGGTCTGCGCCCCATCTGGAACGTACCTTTCCAGATGCGGATTGGTGACTACAGCAGCGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTtaga                                                                                                                  <  2:138509/150‑3 (MQ=255)
                            gTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCg                                                                                                             >  1:531363/1‑150 (MQ=255)
                             ttGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCgg                                                                                                            >  2:404436/1‑150 (MQ=255)
                              tGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCggg                                                                                                           <  2:162460/150‑1 (MQ=255)
                               gCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGt                                                                                                          <  2:139903/150‑1 (MQ=255)
                                         tttACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTg                                                                                                >  4:174477/1‑150 (MQ=255)
                                           tACGGTTTCATCGAAAGTTGCTGGGTAGCGTTTTGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCg                                                                                              >  2:564979/1‑150 (MQ=255)
                                            aCGGTTTCATCGAACGTTGCCTGGTACCGCTTGGTGTCGACCACATCTTGTACGTACCTTCCATGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGCTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACGCGACttagg                                                                                             <  2:469013/150‑5 (MQ=255)
                                              ggTTTCATCGAAAGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGt                                                                                           >  1:306732/1‑150 (MQ=255)
                                                    aTCGAACGTTGCCTGGTACCGTTTGGAGTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAAAGCAGCAGGTTAGGATTTCCACGGCGACACTCCGCGTTCTATGGCGGGTGACCCGACTGCGGAAAAACTg                                                                                     >  4:161869/1‑150 (MQ=255)
                                                     tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                                                                    >  3:376482/1‑150 (MQ=255)
                                                     tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                                                                    >  4:372110/1‑150 (MQ=255)
                                                        aaCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtg                                                                                 >  1:126031/1‑150 (MQ=255)
                                                          cGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggt                                                                               <  1:568572/150‑1 (MQ=255)
                                                            ttGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTAGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggtgg                                                                             >  3:280038/1‑150 (MQ=255)
                                                               ccTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAGCGCAGCAGGTTAGGTTTTCCGCGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCtt                                                                          <  2:81904/150‑1 (MQ=255)
                                                                         tttGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaa                                                                >  3:42491/1‑150 (MQ=255)
                                                                         tttGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTAAATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaa                                                                >  1:55168/1‑150 (MQ=255)
                                                                               cTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTAc                                                          <  3:287714/150‑1 (MQ=255)
                                                                               cTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTAc                                                          <  1:316773/150‑1 (MQ=255)
                                                                                tGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACg                                                         >  3:160059/1‑150 (MQ=255)
                                                                                   accacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTc                                                      >  3:21856/1‑150 (MQ=255)
                                                                                    ccacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGTCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCt                                                     >  4:290039/1‑150 (MQ=255)
                                                                                    ccacATCCGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGTAGGTTCGGTTTTCCACGGTGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCt                                                     <  1:438799/150‑1 (MQ=255)
                                                                                      ccaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGc                                                   <  1:278993/149‑1 (MQ=255)
                                                                                           tGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCt                                              >  2:409314/1‑150 (MQ=255)
                                                                                           tGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCt                                              <  2:450719/150‑1 (MQ=255)
                                                                                             gAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc                                            <  1:181505/150‑1 (MQ=255)
                                                                                             gAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc                                            <  2:182127/150‑1 (MQ=255)
                                                                                                cGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGc                                         >  1:236871/1‑150 (MQ=255)
                                                                                                cGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGc                                         >  4:236730/1‑150 (MQ=255)
                                                                                                 gTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTTCCGCa                                        >  2:375588/1‑150 (MQ=255)
                                                                                                   aCCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAAt                                      <  2:43100/150‑1 (MQ=255)
                                                                                                         ccAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTAt                                <  1:512460/150‑1 (MQ=255)
                                                                                                         ccAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTAt                                <  1:326456/150‑1 (MQ=255)
                                                                                                         ccAGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTAt                                <  1:359038/150‑1 (MQ=255)
                                                                                                           aGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCt                              >  2:271713/1‑150 (MQ=255)
                                                                                                           aGATGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGCACGGTCTGCCAGCTGCCGCAATTGCTATCt                              >  2:515395/1‑150 (MQ=255)
                                                                                                             aTGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTgg                            <  4:145423/150‑1 (MQ=255)
                                                                                                             aTGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTgg                            <  3:131197/150‑1 (MQ=255)
                                                                                                              tGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGc                           >  1:48183/1‑150 (MQ=255)
                                                                                                              tGCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGc                           >  1:337257/1‑150 (MQ=255)
                                                                                                               gCAGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCa                          <  2:246516/150‑1 (MQ=255)
                                                                                                                 aGATTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCAct                        >  3:80722/1‑150 (MQ=255)
                                                                                                                   aTTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCActct                      >  3:120709/1‑150 (MQ=255)
                                                                                                                   aTTGGTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCGGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCActct                      >  1:379213/1‑150 (MQ=255)
                                                                                                                      ggTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCt                   >  1:437384/1‑150 (MQ=255)
                                                                                                                       gTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTa                  >  2:141485/1‑150 (MQ=255)
                                                                                                                       gTGAATACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTa                  >  2:20111/1‑150 (MQ=255)
                                                                                                                          aaTACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAAc               >  2:155757/1‑150 (MQ=255)
                                                                                                                            tACACCAACGCAGCAGGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCa             <  4:349519/150‑1 (MQ=255)
                                                                                                                                      cagcagGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACcgcg   >  1:329152/1‑150 (MQ=255)
                                                                                                                                       ggcagGTTAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACcgcgc  <  2:529326/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGC  >  NC_000913_3_hsa_tpiA/1158417‑1158701

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: