Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,188,555 G→A R1075H (CGT→CAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,188,5550GA98.5% 227.2 / ‑4.3 68R1075H (CGT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (39/28);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

ACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGC  >  NC_000913_3_hsa_tpiA/4188416‑4188698
                                                                                                                                           |                                                                                                                                               
aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTg                                                                                                                                       >  2:442442/1‑150 (MQ=255)
     aGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAAt                                                                                                                                  <  1:483222/150‑1 (MQ=255)
       tAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCg                                                                                                                                >  1:744816/1‑150 (MQ=255)
          gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTg                                                                                                                             >  1:809075/1‑150 (MQ=255)
             gTTTAGAACGCTTTACCCGCATGATCCACGGCCAAACCATTACGCGTCAGACCGGCGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTTGATTCCGCAGAACGTACCTCAGTTGGTAACGATCTGCATTCGGCACTGAAAATCGTTGATg                                                                                                                          <  2:902834/150‑1 (MQ=255)
              ttttGTACGCTTTACTGACATGCTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGACAACCGTTGATGc                                                                                                                         >  2:119850/1‑150 (MQ=255)
              ttttGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAACATCTGCATCCGGCCCTGAAACTCGTTGATgg                                                                                                                         >  2:167879/1‑149 (MQ=255)
                      gCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGGAAATCGTTGATGCTCAGGGTa                                                                                                                 >  1:822294/1‑150 (MQ=255)
                       cTTTACTGACAGGAGCGACGGCCAGACCATGACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTaa                                                                                                                <  2:107620/150‑1 (MQ=255)
                           aCTGACATGATCGACGGCCAACCCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAATGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgac                                                                                                            <  2:149924/150‑1 (MQ=255)
                             tGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAATCCGGCACTGAAAATCGTTGCTGCTCAGGGTAACGACGt                                                                                                          >  2:746244/1‑150 (MQ=255)
                             tGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACGGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGt                                                                                                          >  1:661966/1‑150 (MQ=255)
                               aCATGATAGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTc                                                                                                        >  1:465632/1‑150 (MQ=255)
                                cATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                       >  1:819438/1‑150 (MQ=255)
                                  tGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCCCTGAAAATCGTTGATGCTCCGGGGAACGACGTTCTGa                                                                                                     >  1:509317/1‑150 (MQ=255)
                                  tGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                     <  1:187331/150‑1 (MQ=255)
                                   gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                    >  2:954251/1‑150 (MQ=255)
                                   gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                    <  1:666253/150‑1 (MQ=255)
                                   gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGCTCTGGATTCCGCAGAACGTACCGCAGGTGCTAAAAATCTGCATCCGGCACTGAAAATCGTTGAAGATCAAGGTAACGAAGTTCTgcc                                                                                                    >  1:472308/1‑148 (MQ=255)
                                   gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTAGCTGGTGGTTCTGGATTCCGCAGAGCGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGtcatgct                                                                                                    >  2:954202/1‑144 (MQ=255)
                                            ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           <  1:1019538/150‑1 (MQ=255)
                                            ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           <  1:517137/150‑1 (MQ=255)
                                              aGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCg                                                                                         >  1:966760/1‑150 (MQ=255)
                                                 ccATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGata                                                                                      <  1:629342/150‑1 (MQ=255)
                                                  cATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGatat                                                                                     <  2:744816/150‑1 (MQ=255)
                                                     tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGcc                                                                                  >  2:50053/1‑150 (MQ=255)
                                                      aCGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCa                                                                                 <  1:50053/150‑1 (MQ=255)
                                                          gTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACACCGTTCTGTTCCCAGGTCCCGATATGACAccga                                                                             >  1:866663/1‑149 (MQ=255)
                                                          gTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTAGATGCTCAGGGTAACGAGGTTNTGANCCCAGGTAGCGATATGCCCgcgc                                                                             >  1:907744/1‑150 (MQ=255)
                                                                          gACCGGTCTGTCTTCGCTGGTGGTTCTGGATGCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCggg                                                             >  2:485088/1‑150 (MQ=255)
                                                                             cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                    >  1:739688/1‑92 (MQ=255)
                                                                             cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                    <  2:739688/92‑1 (MQ=255)
                                                                              ggTGTGGCTTGGGTGGTGGTTCTGGATTCCGCAGAACGGACCGCAGGTGGGACAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTaaa                                                         <  2:97201/150‑1 (MQ=255)
                                                                              ggTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTaaa                                                         <  2:809075/150‑1 (MQ=255)
                                                                               gtctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATc                                                                                                                                 <  1:569742/77‑1 (MQ=255)
                                                                               gtctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATc                                                                                                                                 >  2:569742/1‑77 (MQ=255)
                                                                               gtctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTGACGACGTTCTGATCCCAGGTACAGATATGCTAGCGCAGTACTTACTGCCGGGTAAAg                                                        >  2:799398/1‑150 (MQ=255)
                                                                               gtctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAg                                                        >  2:828239/1‑150 (MQ=255)
                                                                                      ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCCGCGCAGTACTTCCTGCCGGGTAAAGCGATTGt                                                 >  1:376030/1‑150 (MQ=255)
                                                                                                 ttCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAg                                      >  2:36024/1‑150 (MQ=255)
                                                                                                 ttCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAg                                      >  2:940791/1‑150 (MQ=255)
                                                                                                      gATTCTGCAGAACCTACGGCAGGTGGTAAATATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGc                                 <  2:497636/150‑1 (MQ=255)
                                                                                                        ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGCAGGTGACGt                               >  2:82400/1‑150 (MQ=255)
                                                                                                        ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGt                               >  2:358225/1‑150 (MQ=255)
                                                                                                         tCCGCGGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                              <  1:358225/150‑1 (MQ=255)
                                                                                                          ccGCAGAACGTACCGCAGGTAGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTAc                             <  2:365650/150‑1 (MQ=255)
                                                                                                           cGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                            <  2:641421/150‑1 (MQ=255)
                                                                                                           cGCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                            <  1:828239/150‑1 (MQ=255)
                                                                                                           cGCAGAACGTACCGCAGGGGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGTAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                            <  2:969014/150‑1 (MQ=255)
                                                                                                            gCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgcg                                                                              <  1:80690/99‑1 (MQ=255)
                                                                                                            gCAGAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgcg                                                                              >  2:80690/1‑99 (MQ=255)
                                                                                                               gAACGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATc                        <  1:442442/150‑1 (MQ=255)
                                                                                                                aaCGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTg                                          >  2:314494/1‑131 (MQ=255)
                                                                                                                aaCGTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTg                                          <  1:314494/131‑1 (MQ=255)
                                                                                                                   gTACCGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGCTCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGTCGGGTAAAGAGATTGTTCAGCTGGAAGATGGCATACAGATCAGct                    >  2:682382/1‑150 (MQ=255)
                                                                                                                       cGCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAATGTTCAGCTGGCAGACGGCGTACAGATCAGCTCTgg                >  2:70595/1‑150 (MQ=255)
                                                                                                                        gCAGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt               <  1:640937/150‑1 (MQ=255)
                                                                                                                          aGGTGGTAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTGCCGATATGACAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGAACAGATCATCTCTGGTGa             >  1:712258/1‑150 (MQ=255)
                                                                                                                            gtggtAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGaca           >  1:987003/1‑150 (MQ=255)
                                                                                                                                tAAAGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAcccg       >  1:6481/1‑149 (MQ=255)
                                                                                                                                 aaaGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTg      <  1:51589/150‑1 (MQ=255)
                                                                                                                                 aaaGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCAAATGTTCAGATGGAAGATAGAGTACATAACAGACCTAGTGGAACCCTg      >  2:975784/1‑150 (MQ=255)
                                                                                                                                  aaGATCTGCATCCGGCACTGGAAATCGTTGACGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTgg     <  2:787865/150‑1 (MQ=255)
                                                                                                                                  aaGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTgg     <  2:788811/150‑1 (MQ=255)
                                                                                                                                  aaGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTgg     >  2:391726/1‑150 (MQ=255)
                                                                                                                                  aaGATCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTgg     >  2:391717/1‑150 (MQ=255)
                                                                                                                                    gATCTGCATCCGGCACCGATATTTTTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGAGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgcg   >  2:891401/1‑150 (MQ=255)
                                                                                                                                     aTCTGCATCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgcgc  >  1:104089/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                               
ACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGC  >  NC_000913_3_hsa_tpiA/4188416‑4188698

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