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,6620GA98.8% 251.9 / NA 82L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (41/39);  minor base C (1/0);  total (42/40)
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

TAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATC  >  NC_000913_3_hsa_tpiA/4111516‑4111807
                                                                                                                                                  |                                                                                                                                                 
tAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAcc                                                                                                                                                 <  1:1334083/149‑1 (MQ=255)
  aaGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGAGTCATCGGCACAAAGGCTTTGCTTTTTTTTTTTTTTTTAAACCTTCAAGCGATTCAGATAGCGCCAGCTTATTCGGTTCAACAGCGAAGGTCACcccc                                                                                                                                               >  2:1390041/1‑149 (MQ=255)
    ggAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCtt                                                                                                                                             <  1:19614/149‑1 (MQ=255)
    ggAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAACGAAGATCAACCCCtt                                                                                                                                             >  2:1301100/1‑149 (MQ=255)
     gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAGCCCttt                                                                                                                                            >  1:80407/1‑149 (MQ=255)
     gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCttt                                                                                                                                            >  1:1347277/1‑149 (MQ=255)
       aGATGGCCGCCCCGCAGGGCAGCAGGGCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                          <  1:60357/149‑1 (MQ=255)
        gATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCg                                                                                                                                         >  2:843064/1‑149 (MQ=255)
         aTGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGc                                                                                                                                        <  2:1337945/149‑1 (MQ=255)
          tGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTcccc                                                                                                                                       >  2:900057/1‑149 (MQ=255)
                  ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGACTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCATTTCGCCGTTGTcgg                                                                                                                               >  2:391102/1‑147 (MQ=255)
                   cGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                              >  2:948910/1‑149 (MQ=255)
                   cGCAGGGCACCACATCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCACATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCcgc                                                                                                                              >  2:1244708/1‑149 (MQ=255)
                     cAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGa                                                                                                                            <  1:373815/149‑1 (MQ=255)
                             cagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATaa                                                                                                                    <  2:17675/149‑1 (MQ=255)
                                gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                                                                                                                 >  2:1026214/1‑149 (MQ=255)
                                 tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                <  1:1344377/149‑1 (MQ=255)
                                   tgtgAAACAGTATAGAGATTAATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTAAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCCCCGTTGTCCGCGACAACATGACGCaag                                                                                                              >  1:523709/1‑147 (MQ=255)
                                    gtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                                                                                                             <  2:1056472/149‑1 (MQ=255)
                                       aaaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACACCATAACGCAGTGCAc                                                                                                          >  2:68210/1‑149 (MQ=255)
                                           aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                >  1:1092601/1‑139 (MQ=255)
                                           aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                                                                                                                <  2:1092601/139‑1 (MQ=255)
                                           aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                      <  2:450945/149‑1 (MQ=255)
                                            gTATAGAGATTCGTTGCGACAAGGGCTTTCCTTTTTGGCATTGATTCAACGCTTCAAGCGATTCAGAGAGCGTCAGCGTAATCGGTTCAACAGCGAAGGTCAACCCCTTTGCGCCGTTGTCCGCGACAACATAGCGCAGTGCAcgtgtt                                                                                                     <  2:966136/149‑6 (MQ=255)
                                               taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtc                                                                                                  <  2:784353/149‑1 (MQ=255)
                                                  agaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg                                                                                               <  1:525639/149‑1 (MQ=255)
                                                    aTTCATCGGCACAAAGGCTTTGCTTTTAGTCATTTATTCAAACCTTCAATCGATTCAGATAGCGCCAGCTTAACCGGTTCAACATCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             <  2:968239/149‑1 (MQ=255)
                                                      tCATCGGCACAAAGGCTTTGCTTTTTTTTTTTTTTTCAAACCTTCAAGCGATTCAGATAGTGGTAGTTTTGTCGGTTCAACAGCGAAGGTCAACCCCTTTTCTCCGTTGTCCGCGACAACATAACGCTGTTCATCTTTTTTCTCGgtgt                                                                                           >  2:596571/1‑149 (MQ=255)
                                                        aTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaa                                                                                         >  2:1346179/1‑149 (MQ=255)
                                                         tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAccc                                                                                                                                                <  2:605222/93‑1 (MQ=255)
                                                         tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAAccc                                                                                                                                                >  1:605222/1‑93 (MQ=255)
                                                         tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             <  2:1211022/144‑1 (MQ=255)
                                                         tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                             >  1:1211022/1‑144 (MQ=255)
                                                         tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaat                                                                                        >  2:646178/1‑149 (MQ=255)
                                                           ggCACAAAGGCTTTGTCTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaataa                                                                                      <  1:321033/149‑1 (MQ=255)
                                                               caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATTGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTCTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtt                                                                                  <  2:426824/149‑1 (MQ=255)
                                                                aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgttt                                                                                 <  2:394373/149‑1 (MQ=255)
                                                                       ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACCACATAACGCAGTGCACCTTCTGTCTCGGTGTACTACCGTTTGTTTTTc                                                                          >  2:463363/1‑149 (MQ=255)
                                                                               tGTCATCTGTTCAAACCTTCAAGCGATCCAGATCTTTCCTGCTTAACCGTTTCTACAGCCACGGTCAACCCCTTTTCGCCGTTGCCCGAGCCAACATACCGCGGTGCTCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGACGt                                                                  <  2:21232/149‑1 (MQ=255)
                                                                                 tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                                >  1:520271/1‑149 (MQ=255)
                                                                                  cATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtag                                                               <  1:558549/149‑1 (MQ=255)
                                                                                   aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                              >  1:23473/1‑149 (MQ=255)
                                                                                   aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                              >  1:1362439/1‑149 (MQ=255)
                                                                                   aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACACCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATACCGTTTGTTTTTCCCCGCCGTtcgt                                                              >  2:847434/1‑149 (MQ=255)
                                                                                    tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             <  1:476963/149‑1 (MQ=255)
                                                                                    tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             >  1:994166/1‑149 (MQ=255)
                                                                                    tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             >  1:50127/1‑149 (MQ=255)
                                                                                           aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTg                                                      >  1:250698/1‑149 (MQ=255)
                                                                                            aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                     >  2:1406903/1‑149 (MQ=255)
                                                                                                  cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                               >  2:1129047/1‑149 (MQ=255)
                                                                                                   aaGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCGGttttt                                              >  1:532213/1‑149 (MQ=255)
                                                                                                    aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTATTAGCGACTGCAGtttttt                                             >  2:57418/1‑149 (MQ=255)
                                                                                                     gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                            >  1:1002346/1‑149 (MQ=255)
                                                                                                      cGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                           <  1:1026214/149‑1 (MQ=255)
                                                                                                        aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                         <  1:367929/149‑1 (MQ=255)
                                                                                                           cAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                      <  1:428685/149‑1 (MQ=255)
                                                                                                            aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                     <  1:1346179/149‑1 (MQ=255)
                                                                                                            aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                     <  1:634761/149‑1 (MQ=255)
                                                                                                            aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                     <  1:51471/149‑1 (MQ=255)
                                                                                                              aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAAGAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                   >  1:576652/1‑149 (MQ=255)
                                                                                                              aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCCTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                   <  1:590003/149‑1 (MQ=255)
                                                                                                               tAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCCACCCCTTTTCGCCGTTGTCCGCGACAACATAACGTGGTGCACCTTCTGCCACCCAGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                  <  1:1406903/149‑1 (MQ=255)
                                                                                                                  gggcAGCTTAGTCGGTAAAACCGCGAAGGGCAGCTCCGTTCCCTGGTTGTCAGTCAAAACTTAACGCAGTGCCCTTTCTGTCTCGGTGTAATAACGTTTGTTTTGCCCCGCCGTTAGTAGCGCCTGCTGTTTTTTCTGGCTTTGTGttt                               <  2:69756/146‑1 (MQ=255)
                                                                                                                   gcCAGCTTAATTGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTCGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                              <  2:812513/149‑1 (MQ=255)
                                                                                                                   gcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                              >  1:897714/1‑149 (MQ=255)
                                                                                                                   gcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTACCGCCTGCAGTTTTTTCTGGCTTCGTGtttt                              >  1:900013/1‑149 (MQ=255)
                                                                                                                    ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGc                                                                                                                                        >  2:1118361/1‑42 (MQ=255)
                                                                                                                    ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGc                                                                                                                                        <  1:1118361/42‑1 (MQ=255)
                                                                                                                        cTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCCGTTTTTTCTGGCTTTGTGTATTGGTCa                         >  1:909750/1‑149 (MQ=255)
                                                                                                                           aaTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                      <  1:900057/149‑1 (MQ=255)
                                                                                                                              gggTGCAACAGCGAAGGGCAACCCTTTTTCGCCGTGGTCCGCGCCAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                   <  2:576652/148‑1 (MQ=255)
                                                                                                                                gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                 <  1:948910/149‑1 (MQ=255)
                                                                                                                                 ttCACCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAAGATACCCCAGTGCACCTTCTGTCTCGGTGTAATACCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCCTTAATgggg                <  1:1086937/149‑1 (MQ=255)
                                                                                                                                    aaCAGCGAAGGTCAACCGCTTTTCGCCGTTGTCCGCGAAAACATAACGCAGTGCACCTTCTTTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTGCAGTTTTTTCCGGCTTTGTGTTTTGTTCCTTCATcggtttg             >  2:1279482/1‑145 (MQ=255)
                                                                                                                                    aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCGTTCATGGGGtag             >  1:1255/1‑147 (MQ=255)
                                                                                                                                    aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa             >  1:1282611/1‑149 (MQ=255)
                                                                                                                                     aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             <  2:1073848/100‑1 (MQ=255)
                                                                                                                                     aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                             >  1:1073848/1‑100 (MQ=255)
                                                                                                                                        gCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGGGTTTTGGTCATTAATGGGGGGAgtgt         >  2:75219/1‑149 (MQ=255)
                                                                                                                                            aGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg     <  1:1163528/149‑1 (MQ=255)
                                                                                                                                             ggTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTCGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGa    <  1:341394/149‑1 (MQ=255)
                                                                                                                                               tCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc  >  2:7082/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATC  >  NC_000913_3_hsa_tpiA/4111516‑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: