Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,946,2630AC41.7% 142.2 / 126.8 96noncoding (2560/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base C (18/22);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGG  >  NC_000913_3_hsa_tpiA/3946118‑3946377
                                                                                                                                                 |                                                                                                                  
gCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAATAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGCTGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTa                                                                                                                >  2:590296/1‑150 (MQ=12)
  cATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                                                                              <  3:18149/150‑1 (MQ=32)
   gtCGCTCAACGGATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAg                                                                                                             <  3:353091/149‑1 (MQ=2)
    tCGCTCAACGGATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAgc                                                                                                            <  2:59153/150‑2 (MQ=2)
     cGCTCAACGGATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAgcg                                                                                                           <  1:558599/150‑3 (MQ=2)
      gCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTgg                                                                                                          >  1:289513/1‑150 (MQ=32)
       cTCAACGGATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGt                                                                                                         <  1:172524/150‑1 (MQ=2)
         cAACGGATAAAAGGTACTCCGGGGATAACAGGATGATACCGCACAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCACAGCACAGCCTGGAGGGGAGGTAGGTCCCAAGGGTATGGGTGTTCGCCATTTAAAGTGGTAc                                                                                                       >  1:323216/1‑150 (MQ=11)
          aaCGGATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTAcg                                                                                                      >  2:365040/1‑150 (MQ=2)
          aaCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                                                                      >  4:369230/1‑150 (MQ=32)
           aCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgc                                                                                                     <  3:82293/150‑1 (MQ=32)
            cGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                                                                    >  3:277398/1‑150 (MQ=32)
              gATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAg                                                                                                  >  1:223752/1‑150 (MQ=2)
              gATAAAAGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAg                                                                                                  >  1:603356/1‑150 (MQ=2)
                  aaaGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTgg                                                                                              <  2:300028/150‑1 (MQ=2)
                  aaaGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                                                                              <  2:244793/150‑1 (MQ=32)
                   aaGTTACCCCGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTggg                                                                                             <  2:457772/150‑1 (MQ=2)
                    aGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                                                                            <  3:97863/150‑1 (MQ=17)
                    aGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                                                                            >  2:79765/1‑150 (MQ=32)
                    aGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                                                                            >  4:205678/1‑150 (MQ=32)
                      ttaccccGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGttc                                                                                          <  2:588377/145‑2 (MQ=2)
                      gTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                                                                          >  4:340177/1‑150 (MQ=32)
                       tACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTa                                                                                         <  4:230210/150‑1 (MQ=32)
                        gcTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAg                                                                                        <  4:227873/149‑1 (MQ=17)
                           ccGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTCGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAAc                                                                                     >  3:75576/1‑150 (MQ=2)
                           ccGGGGATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAAc                                                                                     >  3:149081/1‑150 (MQ=2)
                           ccGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAc                                                                                     <  3:15999/150‑1 (MQ=32)
                            cGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcg                                                                                    >  3:122125/1‑150 (MQ=32)
                            cGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcg                                                                                    >  4:362616/1‑150 (MQ=32)
                              gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                                                                                  <  2:502000/150‑1 (MQ=32)
                               ggATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAAcgtcg                                                                                 <  2:123481/150‑1 (MQ=2)
                                gATAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAAcgtcgt                                                                                <  1:257224/150‑1 (MQ=2)
                                 aTAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTg                                                                               <  1:598912/150‑1 (MQ=2)
                                 aTAACAGGCTTATCTCCCCCAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTg                                                                               <  2:280325/150‑1 (MQ=2)
                                      aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACa                                                                          <  4:252181/150‑1 (MQ=32)
                                         cTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGtt                                                                       >  2:394203/1‑150 (MQ=32)
                                              aCCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc                                                                  >  3:217391/1‑150 (MQ=32)
                                               tcccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTcc                                                                 <  2:43747/147‑1 (MQ=2)
                                               tcccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTcc                                                                 >  2:86425/4‑150 (MQ=2)
                                               ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCTAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc                                                                 <  1:219079/150‑1 (MQ=17)
                                                cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc                                                                >  3:197942/1‑150 (MQ=32)
                                                cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCGGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGCACGTCGTGAGACAGTGCGGTccc                                                                >  4:387485/1‑150 (MQ=255)
                                                 gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt                                                               <  1:431023/150‑1 (MQ=32)
                                                 gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt                                                               <  3:297374/150‑1 (MQ=32)
                                                 gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt                                                               <  4:92803/150‑1 (MQ=32)
                                                 ccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGttc                                                                                          >  2:82068/2‑122 (MQ=2)
                                                 ccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGttc                                                                                          >  3:118202/2‑122 (MQ=2)
                                                 ccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGttc                                                                                          <  4:118202/122‑2 (MQ=2)
                                                 ccccAAGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGttc                                                                                          <  1:82068/122‑2 (MQ=2)
                                                     aaGAGTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCt                                                           <  3:316142/150‑1 (MQ=11)
                                                        ggTTCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCGCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTgtt                                                        <  4:270598/149‑3 (MQ=11)
                                                          ttCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGt                                                      >  1:441675/1‑150 (MQ=32)
                                                          ttCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTAGTGAGACAGTTCGGTCCCTATCTGCCGt                                                      <  1:186951/150‑1 (MQ=17)
                                                           tCACATCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTgttgtg                                                     >  4:329353/1‑145 (MQ=2)
                                                             atatCGACGGCGGTGTTTGGCACCTCCATGTCGGCTCATCACCTACTGGAGATGCAGTAGGTCAACAGGCTATGGATCTTCGAACTTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTggg                                                   <  3:195237/150‑1 (MQ=11)
                                                             acatCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCAAAGGGGTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTggg                                                   >  3:304764/3‑150 (MQ=2)
                                                               atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGgcg                                                 >  2:545042/1‑150 (MQ=32)
                                                                tCGACGGGGAGGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTCCCAAGAGTTGGGCAGTACGCCCATTAAAGCGGTATGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGGTCCTATCTGTTGTGGgcgt                                                >  2:382759/1‑149 (MQ=2)
                                                                 cGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCt                                               <  3:4665/150‑1 (MQ=32)
                                                                       ggaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaa                                         <  4:173905/146‑1 (MQ=2)
                                                                        ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGATGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAAc                                        >  4:73346/1‑150 (MQ=17)
                                                                        gaggtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGaaaat                                        >  3:353461/4‑149 (MQ=2)
                                                                          tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTg                                      >  1:513408/1‑150 (MQ=32)
                                                                           gtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGa                                     >  2:225172/1‑150 (MQ=2)
                                                                           gtTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGa                                     >  2:145338/1‑150 (MQ=11)
                                                                             tttgCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGaaa                                   <  1:454594/147‑3 (MQ=2)
                                                                              tGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAggg                                  >  2:515926/1‑150 (MQ=32)
                                                                                gCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCCTTAAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAggggg                                >  1:242011/1‑150 (MQ=17)
                                                                                gCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGCGACAGCTCGGTCTTTATCTGCCGTGGGCGCTGGAGAACTGAggggg                                <  2:481104/150‑1 (MQ=255)
                                                                                gCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAggggg                                >  4:120546/1‑150 (MQ=32)
                                                                                  aCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgc                              >  4:181314/1‑150 (MQ=32)
                                                                                     tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgctgc                           <  3:100582/150‑1 (MQ=32)
                                                                                     tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgctgc                           <  3:100620/150‑1 (MQ=32)
                                                                                     tCGATGTCGGCTCATCACATCCTGGAGATGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGGGAACTGAGGGGGgctgc                           >  4:364935/1‑150 (MQ=17)
                                                                                       gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTc                         <  4:377645/150‑1 (MQ=32)
                                                                                         tGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCt                       >  3:303070/1‑150 (MQ=17)
                                                                                           tCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAg                     >  2:220458/1‑150 (MQ=32)
                                                                                            cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGt                    <  1:181616/150‑1 (MQ=32)
                                                                                             ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCCGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCCGGCGAACTGAGGGGGGCCGCTCCTAGTa                   >  4:360637/1‑150 (MQ=11)
                                                                                             ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTa                   <  3:112498/150‑1 (MQ=32)
                                                                                                tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACga                <  4:380822/150‑1 (MQ=32)
                                                                                                 catcgCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACgag               <  4:237773/150‑1 (MQ=2)
                                                                                                  atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGTACTGAGGGGGGCTGCTCCTAGTACgaga              <  1:324959/150‑1 (MQ=17)
                                                                                                  atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGAAGAACTGAGGGTGGCTGCTCCTAGTACgaga              <  2:243856/150‑1 (MQ=17)
                                                                                                   tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCACTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGagag             >  2:576860/1‑150 (MQ=17)
                                                                                                      caTCCTGGAGCTGAAGTAGGTCACAAGTGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGAc          <  4:42032/150‑1 (MQ=17)
                                                                                                       aTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGagc         >  2:270691/1‑148 (MQ=17)
                                                                                                        tCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCg        >  2:493167/1‑150 (MQ=2)
                                                                                                          cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      <  1:216081/150‑1 (MQ=32)
                                                                                                          cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa      <  2:241091/150‑1 (MQ=32)
                                                                                                             gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTg   >  2:301335/1‑150 (MQ=31)
                                                                                                             gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGGGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGGCCGGAGTg   >  2:356066/1‑150 (MQ=14)
                                                                                                              aGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGAATgg  <  1:342975/150‑1 (MQ=2)
                                                                                                              aGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGAATgg  <  1:342986/150‑1 (MQ=2)
                                                                                                              aGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGa                                                                                                   <  4:307151/53‑1 (MQ=31)
                                                                                                              aGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGa                                                                                                   >  3:307151/1‑53 (MQ=31)
                                                                                                                                                 |                                                                                                                  
GCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGG  >  NC_000913_3_hsa_tpiA/3946118‑3946377

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