Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,544,0211.A72.6% 152.0 / 56.8 62coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (10/7);  new base A (24/21);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_hsa_tpiA/1543875‑1544158
                                                                                                                                                    |                                                                                                                                              
gttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcg‑cgaa                                                                                                                                               <  1:1226130/148‑3 (MQ=255)
     tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑a‑aa‑aag                                                                                                                                      >  2:1089874/1‑147 (MQ=255)
         gttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATaa                                                                                                                                  <  2:701624/148‑1 (MQ=255)
          ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAAc                                                                                                                                 >  2:851923/1‑149 (MQ=255)
            aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGa                                                                                                                               <  1:641447/149‑1 (MQ=255)
             aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑A‑AA‑AGGATAACGAc                                                                                                                              >  1:246619/1‑149 (MQ=255)
                gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAt                                                                                                                           >  1:889920/1‑149 (MQ=255)
                   aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTGCCTTTATTTTCCCAAAAGTCCGGATTTTCTGGTTTCCAGTCGCCA‑C‑A‑AA‑A‑GCTAACGACTAtttt                                                                                                                        >  1:551220/1‑149 (MQ=255)
                    cTGCGCTAAATAGAATCCAGACACAGAAGGCAAGAAGTAGACAAATGACTGATATCCAGAGATTTCTTCGAGCAAGATGTTTTCCTTTATTTTCTCAGAAGGCCGGATTTTCTTGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAttttt                                                                                                                       <  1:547324/149‑1 (MQ=255)
                     tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTAtttttt                                                                                                                      >  2:90432/1‑148 (MQ=255)
                     tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGTTTTCCAGTCG‑CGAA‑A‑AA‑A‑TATAACGACTATTTTtc                                                                                                                      >  2:852486/1‑149 (MQ=255)
                             aTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTg                                                                                                              <  1:720697/149‑1 (MQ=255)
                                 cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGt                                                                                                          <  1:470667/149‑1 (MQ=255)
                                    ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                       >  1:585965/1‑149 (MQ=255)
                                        acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCatat                                                                                                   >  1:699609/1‑149 (MQ=255)
                                           aaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                <  2:706229/148‑1 (MQ=255)
                                              caGGCAATAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                             <  2:1239860/148‑1 (MQ=255)
                                                 gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAAGGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                          >  1:351405/1‑149 (MQ=255)
                                                   aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        <  1:952673/149‑1 (MQ=255)
                                                   aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        >  1:130387/1‑149 (MQ=255)
                                                   aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGCTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                        >  1:446951/1‑149 (MQ=255)
                                                      aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                     >  1:214481/1‑149 (MQ=255)
                                                           gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                >  2:1117308/1‑149 (MQ=255)
                                                           gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                >  1:1117744/1‑149 (MQ=255)
                                                            aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACATGCACACAt                                                                               >  1:212989/1‑149 (MQ=255)
                                                              aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CC‑C‑AGAA‑A‑GATCAAGATGTTTGTTGTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             <  1:294347/149‑1 (MQ=255)
                                                              aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTAACAGAAGACAGGCTTTTCTGCTTTCCAGGCC‑AGAA‑A‑AA‑A‑GATAACGACAATTATTTCCATTTTGCAGCGCCAAATTTTTCCTCACATGCACACATTg                                                                             >  2:294347/1‑149 (MQ=255)
                                                                  gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                         <  2:1117744/149‑1 (MQ=255)
                                                                     gcatatCCCGGGATTTCCTCCAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCAGCCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                      <  1:118455/147‑1 (MQ=255)
                                                                       atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa                                                                    >  2:694606/1‑149 (MQ=255)
                                                                          tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATTaaaaaa                                                                 >  2:122497/1‑149 (MQ=255)
                                                                          tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                 >  1:706219/1‑149 (MQ=255)
                                                                           ccAGATATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                <  1:261301/149‑1 (MQ=255)
                                                                           ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                <  1:66971/149‑1 (MQ=255)
                                                                             agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTAttt                                                                                                                         <  1:90241/90‑1 (MQ=255)
                                                                             agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                              <  1:92899/149‑1 (MQ=255)
                                                                             agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTAttt                                                                                                                         >  2:90241/1‑90 (MQ=255)
                                                                                 aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                                          <  1:737513/149‑1 (MQ=255)
                                                                                      ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                     >  1:1015291/1‑149 (MQ=255)
                                                                                           gCAATATGTTATCCTTTATTTACCCAAAAGACCGGCTTTTCTGTTTTCCATTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑a                                               >  2:1045574/1‑149 (MQ=255)
                                                                                             aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAg                                             >  2:602707/1‑149 (MQ=255)
                                                                                                  gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCa                                        <  1:324909/149‑1 (MQ=255)
                                                                                                   ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATTCACACATGTGTAATGCAAACAAGACAAAACCCGAGGA‑AAGGCGCaa                                       >  2:1130845/1‑149 (MQ=255)
                                                                                                   ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCaa                                       >  2:488745/1‑149 (MQ=255)
                                                                                                   ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTGTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCaa                                       <  1:544442/149‑1 (MQ=255)
                                                                                                       ccTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGc                                   <  1:177699/149‑1 (MQ=255)
                                                                                                         tttATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCAGTCG‑CG‑A‑A‑AA‑ATGATAACGACTATTTTTCTCCTTTTGCAGTGCCATATTGTTCCTCACATGCACTCATTGGTCATGAAAAAATGACAAGACACGAAGT‑ACAGCGCAACAgtcg                                 >  2:312262/1‑146 (MQ=255)
                                                                                                         tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGccg                                 >  1:566497/1‑148 (MQ=255)
                                                                                                         tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCa                                 >  2:1065993/1‑149 (MQ=255)
                                                                                                             ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGtta                             >  1:1011879/1‑149 (MQ=255)
                                                                                                               ttCCCAGGAGGCCGGATTTTCTGGTTTCCAGTCG‑AG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGttatt                           <  1:532564/149‑1 (MQ=255)
                                                                                                                  ccAGAAGGCCGGAGTTTCTGGTTACCAGAAG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGa                        <  2:319947/149‑1 (MQ=255)
                                                                                                                     gAAGGCGGGTTTTTCTGGTTTTGAATCG‑TGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGGGT‑AAGGCGCAATAGCCAGTTATTAGAAtt                     <  2:24504/149‑1 (MQ=255)
                                                                                                                       aGGCCTGATTGTCTTGTTTCAAGGCT‑CG‑A‑A‑AAGA‑GATAACGACTATTTTTCTCAATTTGCAGTGCCATATTGTTCCTCATATGCACACAATGGTAAAGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTaa                   <  2:601697/149‑1 (MQ=255)
                                                                                                                          ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             >  2:597364/1‑89 (MQ=255)
                                                                                                                          ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                             <  1:597364/89‑1 (MQ=255)
                                                                                                                             gATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa             >  1:817975/1‑149 (MQ=255)
                                                                                                                                 ttCTGGTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg         <  2:343702/149‑1 (MQ=255)
                                                                                                                                 ttCTGGTTTCAAGTCT‑CG‑A‑ATAA‑A‑GATAACGGCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTTACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATCGCCAGTTATTAGAATTGAGGATGACTTg         <  2:615720/149‑1 (MQ=255)
                                                                                                                                      gTTTCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAA‑AAGCGGTGATGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa    >  1:1029902/1‑149 (MQ=255)
                                                                                                                                       tttCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAACATGATTAAACACGAGAT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa   >  1:603090/1‑149 (MQ=255)
                                                                                                                                       tttCCAGTCG‑CGAA‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa   <  1:1141241/149‑1 (MQ=255)
                                                                                                                                        ttCCAGTCG‑CG‑AAA‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGGAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg  <  2:889920/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AA‑A‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_hsa_tpiA/1543875‑1544158

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