Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,544,022 1 bp→AA coding (39/1389 nt) narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0211.A85.0% 178.5 / 26.6 60coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (7/2);  new base A (30/21);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_hsa_tpiA1,544,0220CA100.0% 205.0 / NA 60L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (37/23);  total (37/23)

AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543882‑1544171
                                                                                                                                             |||                                                                                                                                                     
agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aagat                                                                                                                                               >  1:1024430/1‑145 (MQ=255)
  aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATaa                                                                                                                                              <  2:177118/149‑1 (MQ=255)
             gcgcgaTAAATAGCATCCAGACAAAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTAttttt                                                                                                                                   <  2:357988/147‑1 (MQ=255)
              tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCGATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtc                                                                                                                                  >  1:496008/1‑149 (MQ=255)
               gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTCTTTATTCGCGAA‑AAA‑AGATAACGACTATTTTtct                                                                                                                                 <  1:519835/149‑1 (MQ=255)
               gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGGAGACAACTGACTGGTATCCAGAGATTCCTTCGATCAATATGTTTTCCTTTATTTTCCCAGAAGACCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtct                                                                                                                                 >  2:436570/1‑149 (MQ=255)
                cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtctc                                                                                                                                <  1:529495/149‑1 (MQ=255)
                 gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAATGACTATTTTTCTCa                                                                                                                               >  2:1031380/1‑149 (MQ=255)
                    aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCtttt                                                                                                                            >  2:451399/1‑149 (MQ=255)
                               aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                 >  2:1040162/1‑149 (MQ=255)
                                    caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                            <  2:1024430/149‑1 (MQ=255)
                                        aGGCAAGAAGTAGACAACTGACTGATATCCAGATATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                        <  1:459568/149‑1 (MQ=255)
                                         ggCAAGCAATAGAACACTGCTTGATATCCAGAGCTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                       <  2:524306/149‑1 (MQ=255)
                                            aagaagTAGACTACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    <  1:75572/149‑1 (MQ=255)
                                            aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                    >  2:32875/1‑149 (MQ=255)
                                             agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                   >  2:683540/1‑149 (MQ=255)
                                                    gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATCTTGTTCCTCACATGcacaca                                                                                            >  2:283074/1‑149 (MQ=255)
                                                    gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                            >  1:24213/1‑149 (MQ=255)
                                                     aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:102084/1‑149 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                         >  1:704673/1‑149 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATATTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                         >  1:921639/1‑149 (MQ=255)
                                                         ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGACATATTGTTCCTCCCATGCACACATTGGt                                                                                       >  1:530159/1‑149 (MQ=255)
                                                           gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     <  2:33938/149‑1 (MQ=255)
                                                           gacaaatatCCAGAGATTACTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTAAAAAAAAAAAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                     <  2:33943/144‑1 (MQ=255)
                                                              tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                                  <  1:536027/149‑1 (MQ=255)
                                                                   tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                             >  2:717876/1‑149 (MQ=255)
                                                                   tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                             >  1:870938/1‑149 (MQ=255)
                                                                    ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  1:1040162/149‑1 (MQ=255)
                                                                     cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                           >  2:973030/1‑149 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          >  1:754282/1‑149 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGGAATGAAAAAAAGa                                                                          >  1:377621/1‑149 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTGTCCCGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTCCAGTGCCATATTGTTCCTCACATGCACACATTTGTAATGAAAAAAAGa                                                                          >  2:814411/1‑149 (MQ=255)
                                                                        agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                        <  1:320299/149‑1 (MQ=255)
                                                                          aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaca                                                                      >  2:782285/1‑147 (MQ=255)
                                                                              cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                  <  1:717876/149‑1 (MQ=255)
                                                                               ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATCAAAAAAGTACAAAACACg                                                                 >  1:402104/1‑149 (MQ=255)
                                                                                tcaggCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                <  1:511215/145‑1 (MQ=255)
                                                                                  gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                              >  2:332528/1‑149 (MQ=255)
                                                                                     cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                           <  1:973030/149‑1 (MQ=255)
                                                                                            ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    >  1:957126/1‑149 (MQ=255)
                                                                                            ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    <  1:272969/149‑1 (MQ=255)
                                                                                            ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                    <  1:272973/149‑1 (MQ=255)
                                                                                                 cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                               >  1:635281/1‑149 (MQ=255)
                                                                                                  tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  1:592074/1‑149 (MQ=255)
                                                                                                  tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  2:553639/1‑149 (MQ=255)
                                                                                                  tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTAATTTTGCAGTGCCATATTGTTCGTCACATGCACACATTGGTAGGGAAAAAAAGACAAAGCACGAGGTAAGGCGCAATAGGCa                                              >  1:165825/1‑149 (MQ=255)
                                                                                                   ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCCTATTGTTCCTCACATGCACACATGGGTAATGGAAAAAAGACAAAAAACGAGGTAGGGCTGACTAGccgg                                             >  1:25623/1‑147 (MQ=255)
                                                                                                   ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                      <  2:369937/92‑1 (MQ=255)
                                                                                                   ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                      >  1:369937/1‑92 (MQ=255)
                                                                                                   ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             <  1:32875/149‑1 (MQ=255)
                                                                                                     aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACCTTGGTACACAAAAAAAGACAAAACAAGAGGGAAGGCGCAATAGCCAGtt                                           >  1:554760/1‑149 (MQ=255)
                                                                                                              gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGACAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGacat                                  >  2:217947/1‑146 (MQ=255)
                                                                                                              gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAtt                                  >  1:273814/1‑149 (MQ=255)
                                                                                                                    cGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAggct                            >  2:211902/1‑147 (MQ=255)
                                                                                                                             tGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTCTTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                   <  1:1044804/149‑1 (MQ=255)
                                                                                                                             tGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                   >  2:723500/1‑149 (MQ=255)
                                                                                                                                tttCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                <  1:723500/149‑1 (MQ=255)
                                                                                                                                    cAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            >  2:299624/1‑149 (MQ=255)
                                                                                                                                    cAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            >  1:904180/1‑149 (MQ=255)
                                                                                                                                       tCGCG‑A‑AAACAGATAACGACTATTTTTCTCATTTTGGAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa         <  2:25623/149‑1 (MQ=255)
                                                                                                                                         cgcgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACATGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt        <  1:837326/143‑1 (MQ=255)
                                                                                                                                             aaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          <  2:634378/45‑1 (MQ=255)
                                                                                                                                             aaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          >  1:634378/3‑47 (MQ=255)
                                                                                                                                              aaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa   <  2:273814/148‑1 (MQ=255)
                                                                                                                                               aaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  <  1:553639/149‑1 (MQ=255)
                                                                                                                                             |||                                                                                                                                                     
AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_hsa_tpiA/1543882‑1544171

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑