Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,1690T.39.6% 109.0 / 79.9 54intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (2/30);  new base . (1/20);  total (4/50)
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.38e-01

AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGG  >  DH1/936020‑936315
                                                                                                                                                      |                                                                                                                                                    
aaTGGCTCTTTTTTATGGGAGACCAGAAACAAAAAAACACGTTTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCt                                                                                                                                                      <  6:763634/150‑1 (MQ=255)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                                                                                     <  8:638300/151‑1 (MQ=255)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                                                                                     <  4:550446/151‑1 (MQ=255)
 gcgatctttttttATGGGAGCCCAGAAACAAATAAAGATCTGTTAGGGTGTTAAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                                                                                    <  8:1023808/144‑1 (MQ=17)
 aTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                                                                                    >  8:821134/1‑151 (MQ=255)
   ggCTCTTTTTTTTTGGAGACCAGTTAAAATAAATGACCTTTTAGGTTGTTCAATTATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttac                                                                                                                                                  <  4:877732/151‑3 (MQ=14)
   ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                                                                                   <  5:283552/150‑1 (MQ=255)
   gcgtcTTTTTTATGGGAGCCCAGAAAAAAAAAAACACCCTTTAGGTTGTTAAATTATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                                                                                   <  8:993564/146‑2 (MQ=17)
     cgcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACt                                                                                                                                                <  4:44344/151‑1 (MQ=255)
      tctattttATGGGAGACCAGTAAAAAAAAAACACGCGTTATGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                                                                                    <  8:717640/142‑1 (MQ=18)
      gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttacatg                                                                                                                                               >  7:718434/1‑146 (MQ=255)
      gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtc                                                                                                                                              <  4:766333/151‑1 (MQ=255)
       cTTTTTTATGGGAGACCAGTAAAAAAAAATGATCATTTAGGTTGTTAATTAATTGGTTGCGGGGGACGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttactt                                                                                                                                                <  2:523524/149‑5 (MQ=14)
        ttttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                                                                                    >  7:1025090/1‑144 (MQ=255)
          ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                                                                                            <  8:720990/150‑1 (MQ=255)
          ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                                                                                            <  1:720300/150‑1 (MQ=255)
           ttttTGGGAGCCCAGAAACCAAAAATGCCCGGTTAGCGTGTTCATTAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTAGCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcat‑ca                                                                                                                                         <  5:150760/150‑1 (MQ=14)
           tttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat‑c                                                                                                                                          <  5:585548/150‑1 (MQ=255)
           tttATGCGAGACCAGAAACAATAACTGAGGCTTCAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcat‑ca                                                                                                                                         <  7:270500/150‑1 (MQ=14)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat‑ca                                                                                                                                         <  8:1078869/150‑1 (MQ=255)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat‑ca                                                                                                                                         <  3:1011771/150‑1 (MQ=255)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat‑ca                                                                                                                                         <  5:1086591/150‑1 (MQ=255)
            ttATGGGAGACCAGAAACAAAAAAACACCCGATAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑Caa                                                                                                                                        <  5:3583/150‑1 (MQ=255)
            gtATGGGAGCCCAGCAACAAAAACACACCCGTTATGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑Caa                                                                                                                                        <  7:822015/149‑1 (MQ=17)
             tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑Caa                                                                                                                                        <  7:427481/150‑1 (MQ=255)
             tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑Caa                                                                                                                                        <  6:333194/150‑1 (MQ=255)
             tATGCGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTATATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑Caaa                                                                                                                                       <  5:683190/150‑1 (MQ=255)
              tttgtagaCCTGATACTAAAAAAGATCTTTTTGGTTGTTAAATTATTGTTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCGGCGTAACCGTACTGC‑TTTACTTCAT‑CAAA‑t                                                                                                                                     <  2:382401/145‑1 (MQ=11)
              ttGGGAGACCAGTAAAAACAACACACTCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑t                                                                                                                                     <  2:926769/149‑1 (MQ=17)
              aTGGTAGACCAGAAACAAATAATGTTTTTTTAGGTTGTTAAATTATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑t                                                                                                                                     <  8:115597/150‑1 (MQ=11)
              aTGGGAGACCCGCAACAAAAACACACCCTTTAGGGTGTTCATTAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑t                                                                                                                                     <  5:595711/150‑1 (MQ=17)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGTCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑Caaa                                                                                                                                       <  4:497895/150‑1 (MQ=255)
                tggAGACCAGAAACAAAAAAACAACCGTTAGTGTTTTCAATAATTGGTTGCGGGGGCCTGATTT‑GAACCGACGACCTTCGTGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACAGCTTTTACTTCAT‑CAAA‑tt                                                                                                                                    <  5:626079/149‑1 (MQ=18)
                gggAGACCAGCAAACAAAAAACACCAGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑ttt                                                                                                                                   <  1:497441/150‑1 (MQ=17)
                gggAGACCAGATACAATAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGTCGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCGACGCCGCGTCTGCGTACTGCTTTTACTTCAT‑GAAA‑tt                                                                                                                                    <  4:866630/150‑1 (MQ=18)
                gggAGACCAGAAACAACAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑ttt                                                                                                                                   <  5:796606/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCGGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑ttt                                                                                                                                   <  5:879285/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑ttt                                                                                                                                   <  8:798945/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑ttt                                                                                                                                   <  3:104253/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑tt                                                                                                                                    <  4:912976/150‑1 (MQ=255)
                gggACACCAGAAACCAAAAAACCCCCGTTTGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑tttt                                                                                                                                  <  3:726178/151‑2 (MQ=14)
                 ggAGACCTGATACAAAAAAACACTCTTTTGGTTGTTATATTATTGGTTGCGGGGGCCGGATTT‑GAACGGACGACGTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑tttta                                                                                                                                 <  8:153591/151‑3 (MQ=11)
                 ggAGACCAGAACCAAAAAAACCCGCGTGAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑tttta                                                                                                                                 <  5:438227/151‑3 (MQ=14)
                 ggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTg                                                                                                                                  <  7:816921/151‑1 (MQ=255)
                     aCCTGAAACAAAAAAACACCCTTTAGGTTGTTAAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑ttttaattg                                                                                                                             <  4:477093/151‑7 (MQ=14)
                     aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTg                                                                                                                              <  2:940963/151‑1 (MQ=255)
                      ccAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTgg                                                                                                                             <  2:498949/151‑1 (MQ=255)
                       cAGCAACAACTAAACACGCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAA‑ttttaattggt                                                                                                                           <  2:886139/151‑9 (MQ=11)
                       cAGAAACAACAAAACACCCTTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcat‑c                                                                                                                                          <  7:859148/137‑1 (MQ=255)
                       cAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGt                                                                                                                            <  8:175150/151‑1 (MQ=255)
                       cAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtaat‑c                                                                                                                                          >  8:858232/1‑137 (MQ=255)
                        acaaccAAAAAAACACGCGTCAGCTTGTTCATTTATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAATTTTAATTGGtt                                                                                                                           <  5:593323/146‑1 (MQ=11)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGtt                                                                                                                           <  8:824962/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGtt                                                                                                                           <  8:508518/151‑1 (MQ=255)
                         gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTg                                                                                                                          <  7:284201/151‑1 (MQ=255)
                         ccaaCAAAAAAACCCCCTTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCAT‑CAAATTTTAATTGGTTg                                                                                                                          <  5:543345/149‑1 (MQ=11)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGc                                                                                                                         <  3:480166/151‑1 (MQ=255)
                            aCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑TAACATACGAACAACTCGTTAAGACAGAGACGAGCGACAAGGCTGCGCCACCCCGCGTCTCCGTACTTCATTTAC‑TCATGCAAA‑TTTTATTGGTTgggg                                                                                                                       >  7:421702/1‑151 (MQ=16)
                             cAAAAAAACACCCGTTAGGGTGTTCAATACTTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCggg                                                                                                                      <  2:953756/151‑1 (MQ=255)
                             aaaaaaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCggg                                                                                                                      <  4:628186/150‑1 (MQ=255)
                              aaaaaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCgggg                                                                                                                     <  1:831928/151‑1 (MQ=255)
                                  aaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTAGAGCGCAAGCCAGGCGGG‑TATGAGAACCCCGAGATACCAGGCTGATCCCACCCGCGTCCCCGTAAACATATAAATAAAT‑CAAA‑GTTGGTGGGTTGCGGGGGcgg                                                                                                                 >  8:510770/1‑149 (MQ=255)
                                     cacCCTTTATGTTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtca                                                                                                                                             <  2:483760/121‑1 (MQ=11)
                                     cacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttaattca                                                                                                                                              >  1:484152/1‑115 (MQ=255)
                                                  gtTCATTAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                                                                                   <  8:654945/103‑2 (MQ=18)
                                                  gtTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                                                                                   >  7:655673/1‑103 (MQ=39)
                                                             gtttGCGGGGGCCGGATGT‑GATCGGACTACCGTGGTTTTATCAGCGCGACGAGCTACGAGGCTGCTCCACACCGCGTCACCGTACTGCTTTTAGTTGCT‑GATA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTa                                                                                       <  8:609952/148‑1 (MQ=14)
                                                                                                                                                  cTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCAGTAATGCGATAATTCATCAAATTATAtttccttgacgggtactt  >  8:422928/1‑136 (MQ=14)
                                                                                                                                                     cTTTTACTTCAT‑AAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacg                                                                            <  1:497373/74‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTT‑GAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCAT‑CAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGG  >  DH1/936020‑936315

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