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_tpiA4,111,6620GA99.1% 370.3 / NA 114L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (60/53);  minor base C (0/1);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

AAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATC  >  NC_000913_3_hsa_tpiA/4111518‑4111807
                                                                                                                                                |                                                                                                                                                  
aaGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTACTCGATTAAACAGCGAAGGTCAACCCCt                                                                                                                                               >  1:270408/1‑150 (MQ=255)
 aGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCtt                                                                                                                                              >  1:787377/1‑150 (MQ=255)
 aGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCtt                                                                                                                                              >  1:997428/1‑150 (MQ=255)
   gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCtttt                                                                                                                                            >  2:284789/1‑150 (MQ=255)
         ggCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGt                                                                                                                                      >  2:29325/1‑150 (MQ=255)
          gccgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGtt                                                                                                                                     <  2:860298/150‑1 (MQ=255)
           ccgccCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTg                                                                                                                                    <  1:732886/150‑1 (MQ=255)
              ccccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTcc                                                                                                                                <  1:29325/150‑1 (MQ=255)
                ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                              >  1:347116/1‑150 (MQ=255)
                 cGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgcg                                                                                                                             >  1:346590/1‑150 (MQ=255)
                 cGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgcg                                                                                                                             >  2:675076/1‑150 (MQ=255)
                  gCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                            <  2:418838/150‑1 (MQ=255)
                  gCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                            <  2:748203/150‑1 (MQ=255)
                    aGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAATCCCTTTTCGCCGTTTTTCGCGaca                                                                                                                          >  2:380552/1‑150 (MQ=255)
                      ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                        >  2:32034/1‑150 (MQ=255)
                      ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaac                                                                                                                        >  2:696608/1‑150 (MQ=255)
                                 tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                             >  2:436229/1‑150 (MQ=255)
                                 tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                             >  1:1041572/1‑150 (MQ=255)
                                 tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                             <  1:921182/150‑1 (MQ=255)
                                   ggAAATATTACAGAAATTCACCGGCACAAAGGCTTTGCTTTTGGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCCGTTTACCAGCGAATGTTAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                                           <  2:196959/149‑1 (MQ=255)
                                      aaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCt                                                                                                        >  1:164310/1‑150 (MQ=255)
                                       aCAGTATAGCGATTCATCGGCACAAAGGTTTTGCTTTTTGTAATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCtt                                                                                                       <  2:379481/150‑1 (MQ=255)
                                       aCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCtt                                                                                                       <  1:308311/150‑1 (MQ=255)
                                        cAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                      <  2:413107/150‑1 (MQ=255)
                                         agcttatagagTCACCGCCACAACGGTTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                     <  1:144894/143‑1 (MQ=255)
                                         aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                     <  2:167787/150‑1 (MQ=255)
                                         aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                     <  1:1029894/150‑1 (MQ=255)
                                              agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                                >  1:424769/1‑150 (MQ=255)
                                               gcgttagaTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                               <  2:771389/143‑1 (MQ=255)
                                               gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                               <  2:1095555/150‑1 (MQ=255)
                                                 gaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAATCTTCATGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGCTCAACCCCTTTTCGCCGTTGCCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             <  2:1167401/150‑1 (MQ=255)
                                                       tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                     >  1:38296/1‑136 (MQ=255)
                                                       tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                     <  2:38296/136‑1 (MQ=255)
                                                       tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaata                                                                                       >  2:36972/1‑150 (MQ=255)
                                                         ctcacaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                     <  1:1135385/148‑1 (MQ=255)
                                                          gCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACg                                                                                    >  2:925350/1‑150 (MQ=255)
                                                          gCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACg                                                                                    >  1:1023704/1‑150 (MQ=255)
                                                          gCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGCTGCACTAAct                                                                                    >  2:923932/1‑149 (MQ=255)
                                                            acaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAGCgtt                                                                                  >  1:347713/1‑150 (MQ=255)
                                                            acaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTAAAACTTTTTTTGCCTGGTGCCACCACAACATAACGCAGTGCACCTTCTGTCTCGTTGTAATAACgag                                                                                  >  2:1135385/1‑148 (MQ=255)
                                                             ccaaGGCCTTGCTGTTTGTCAATTATTAAAACCGTTAAGCTTTTCAGATAGCCCCAGCTTATTCGGTCCAACAGCGAAGGGCCACCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgttt                                                                                 <  2:351162/148‑1 (MQ=255)
                                                             caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgttt                                                                                 <  1:199203/150‑1 (MQ=255)
                                                              aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACTCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttt                                                                                >  2:949717/1‑149 (MQ=255)
                                                                aGGCTTTGCTTTTTGTCATTTATTCATACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGTCAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgtt                                                                              <  1:1068447/150‑1 (MQ=255)
                                                                 gtcgttgcattTTGTCAATTATTCAAAACGTCAAGCTTTTTTTATAGCTCCTGCTTAATCGGTTCAACAACGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGGCTCGGTGTAATAACgtttgttt                                                                             <  2:524082/146‑1 (MQ=255)
                                                                   ctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                           <  1:459277/150‑1 (MQ=255)
                                                                     ttgctttTTGTCATTTATTCAGACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                         >  2:1040102/1‑150 (MQ=255)
                                                                         ttgttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                     <  2:346590/147‑1 (MQ=255)
                                                                           tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                   <  1:36972/150‑1 (MQ=255)
                                                                            ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGATCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCATCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                  >  2:127231/1‑150 (MQ=255)
                                                                            ttGGCATTAATTAAGACCTTCAACCGATTTAAATAGCGCCAGTTTAACCGGTTAAAAAACGAAGGTGAACCTCTTTTCGCCCTTGTCTGCGACATTATAACGCGGTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                  <  2:45665/150‑1 (MQ=255)
                                                                               tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATTACGTTTGTTTTTCCCCGCCGTtag                                                               >  2:1018731/1‑150 (MQ=255)
                                                                               tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtag                                                               >  2:217679/1‑150 (MQ=255)
                                                                                 aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             >  1:1112817/1‑150 (MQ=255)
                                                                                 aTTTATTCAAACCTTCAAGCGATTCAGATAGAGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTTCTTTTCCCCTCGTTtagta                                                             >  2:554076/1‑150 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                      <  1:602465/108‑1 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                      >  2:602465/1‑108 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACGTTTGTTTTTCCCCGCCGTtagtag                                                            >  1:1100310/1‑150 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                     <  1:402287/125‑1 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                     >  2:402287/1‑125 (MQ=255)
                                                                                  tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagtag                                                            >  1:305646/1‑150 (MQ=255)
                                                                                   ttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCCTtagtttc                                                           >  2:944265/1‑147 (MQ=255)
                                                                                      ttCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAAGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                        >  1:774504/1‑150 (MQ=255)
                                                                                        cAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                      <  1:436229/150‑1 (MQ=255)
                                                                                         aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGCTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                     >  1:900702/1‑150 (MQ=255)
                                                                                          aaCCTTCAAGCTATTCATGTCGCGCCAGCTTAAGCGGTTCAACTGCGAAGGTCAACCCCTTTTCGCCGTAGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                    <  2:928598/150‑1 (MQ=255)
                                                                                          aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGt                                                                                                                                      >  2:97415/1‑69 (MQ=255)
                                                                                          aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCCCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCa                                                    >  2:176445/1‑150 (MQ=255)
                                                                                             cTTCATGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAATGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGCGCAGCTTCTGTCTCGGAGTAATAACGTTTGTCTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                 <  1:998099/150‑1 (MQ=255)
                                                                                               tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                               <  1:32034/150‑1 (MQ=255)
                                                                                                  aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttttttt                                            >  1:79089/1‑149 (MQ=255)
                                                                                                   gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATGACGTTTGTTTTTCGAGGCCGTTAGTAGCGACTGAAGgtttttgt                                           >  1:301129/1‑148 (MQ=255)
                                                                                                      aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                                <  1:651045/126‑1 (MQ=255)
                                                                                                           gATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                   >  2:817894/1‑150 (MQ=255)
                                                                                                            aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTAGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCCTtgtc                                  >  2:291981/1‑149 (MQ=255)
                                                                                                             tAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCNCCGTTAGTAGNGNCNGCAGCTTTTTCTTGCGTTgtgt                                 >  1:905569/1‑150 (MQ=255)
                                                                                                               gcggCAGCTTCATAGGGTCAAATGCGAACGAAAACTTCTTTTTGCCGTTTTCAGCGACAACATAACGCGGTGCACCTTCTGTCTCGGTGTAATATCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                               <  2:1002672/150‑1 (MQ=255)
                                                                                                                cgcCAGCTTAATCGGTTCAACGGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                              <  1:677447/150‑1 (MQ=255)
                                                                                                                 gcCCGCTTAATCCGGTCCACAGCGAAGGTCACCCCCCTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATTACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                             <  1:554076/150‑1 (MQ=255)
                                                                                                                   cgggttAATCGGTTCAACAGCGAAGGTACACCCCTTTTCCTTGTTGTCCGTGACAACATAACGCAGTGCACCTTCTGTCCCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                           <  2:203720/146‑1 (MQ=255)
                                                                                                                   cAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCGTTGGGTTTTGGt                           >  1:508939/1‑150 (MQ=255)
                                                                                                                           tgtgTTAAACAGTGAAGGTTCACCCCTTTTCGCGCTTGCCCGCTAGAACTTAACTCAGTGCTCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                   <  2:904936/147‑1 (MQ=255)
                                                                                                                            cGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                  <  1:209906/150‑1 (MQ=255)
                                                                                                                             ggTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCTGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                 <  2:347116/150‑1 (MQ=255)
                                                                                                                             ggTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                 <  1:359094/150‑1 (MQ=255)
                                                                                                                               ttCAACAGCGAACGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt               <  1:521200/150‑1 (MQ=255)
                                                                                                                                 cAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa             <  1:629248/150‑1 (MQ=255)
                                                                                                                                   aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             <  2:627714/68‑1 (MQ=255)
                                                                                                                                   aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             >  1:627714/1‑68 (MQ=255)
                                                                                                                                   aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt           >  1:410139/1‑150 (MQ=255)
                                                                                                                                       caaaGGTCAACCCCTCTTCGCCGTTGTCCGCGAAAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc       <  1:247584/148‑1 (MQ=255)
                                                                                                                                       cGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                           <  2:15288/114‑1 (MQ=255)
                                                                                                                                       cGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                           >  1:15288/1‑114 (MQ=255)
                                                                                                                                        gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg      >  2:727062/1‑150 (MQ=255)
                                                                                                                                        gAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg      <  1:185882/150‑1 (MQ=255)
                                                                                                                                          aGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTTATTATGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGa    >  1:224241/1‑150 (MQ=255)
                                                                                                                                           ggTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                  >  1:537935/1‑119 (MQ=255)
                                                                                                                                           ggTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                  <  2:537935/119‑1 (MQ=255)
                                                                                                                                            gTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc  >  1:275014/1‑150 (MQ=255)
                                                                                                                                            gTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGGGAGTGTTCGGATc  >  1:1042339/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATC  >  NC_000913_3_hsa_tpiA/4111518‑4111807

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: