Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,4410GC25.6% 187.6 / 23.2 133intergenic (+207/+53)ECDH1_RS12810/ECDH1_RS12815tRNA‑Tyr/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (68/22);  new base C (6/25);  total (81/52)
Fisher's exact test for biased strand distribution p-value = 5.48e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

CGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCA  >  DH1/2591293‑2591590
                                                                                                                                                    |                                                                                                                                                     
cGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCGGAACTCTCAAGCg                                                                                                                                                     >  1:1090790/1‑151 (MQ=255)
cGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCTGAAGTCTCAATCg                                                                                                                                                     >  5:117061/1‑151 (MQ=255)
    gAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAGTCTCAACCGAATc                                                                                                                                                 >  3:1090155/1‑151 (MQ=9)
     aaGGGTGAGAACCTTCGACTAAGGTTCGATTCGAGCGAAAGCGAGAGAACGTTGCCACAGGCAACGCCCCGAAGGGTGATGAGCGAAGCGAAGAATAATCCTTCCCCCACCACCATCACTTTCAATACATCCTGAAATCAAAACCGAATcc                                                                                                                                                >  4:931015/1‑151 (MQ=17)
       ggATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCGGAAATCTCAAACGAATCCGc                                                                                                                                              >  6:779981/1‑151 (MQ=9)
          agagaAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAACCTTCCCCCACCACCATCACTTTCAAAAGACCCGGAATTCTCAACCGAATCCGCaa                                                                                                                                            >  3:1086285/2‑150 (MQ=11)
            agaAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAATCTCAAGCGAATCCGCAATc                                                                                                                                          >  8:592980/1‑150 (MQ=16)
                ccTTCGACGAAGGTTCGATTCGCGCGCCAGCGAGAGACCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAata                                                                                                                                     <  3:20360/150‑1 (MQ=255)
                 cTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGGCCCTGAACTCTCAAGCGAATCCGCAATCAAatat                                                                                                                                    >  5:883452/1‑151 (MQ=255)
                 cTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAATTCTCAAACGAATCCGCAATCAAata                                                                                                                                     >  4:363092/1‑150 (MQ=16)
                    cccccaacctcccaCTCGACCGCCACCGAGAGAGCTTTTCCTAACGCAACGACCCCAACGGCTAAGCTCGCATCCCTGAGTAATCCTTCCCCCACCACAATCACTTTCAAAATTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCt                                                                                                                                 <  2:9188/138‑1 (MQ=9)
                     gACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCTCCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATAttatt                                                                                                                                >  8:942467/1‑148 (MQ=255)
                       ccAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGcc                                                                                                                              >  5:532611/1‑151 (MQ=255)
                          aGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCa                                                                                                                            >  5:146659/1‑150 (MQ=255)
                            gTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAAATCTCAACCGAATCCGCAATCAAATATTCTGCCCAAGc                                                                                                                         >  6:566444/1‑151 (MQ=255)
                                gACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCgggg                                                                                                                     >  1:44852/1‑151 (MQ=255)
                                gACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCgggg                                                                                                                     >  1:50742/1‑151 (MQ=255)
                                      cgcgaacccGAGAGAGCGTTGCCTCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGAt                                                                                                               <  3:932675/143‑1 (MQ=255)
                                        cgcCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGagctcgg                                                                                                              >  3:651988/1‑144 (MQ=255)
                                         gcCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAgg                                                                                                                 >  8:705493/1‑146 (MQ=25)
                                           cAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATgagac                                                                                                          >  3:984681/1‑150 (MQ=40)
                                               gagagaGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCtt                                                                                                      >  4:708367/1‑151 (MQ=25)
                                                    agCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGGAGCCCTGAGTAATGCGTCCCCCACCACCAGCACTTTGAAAAGGCCCTGAACTCTCAAGCGAATCCGCAATCAAGGATTCTGCCCGGGCGGGGGAGGAGGAGCAGCttggccg                                                                                                 >  3:1231024/1‑146 (MQ=11)
                                                    agCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGAcc                                                                                                 >  5:1160697/1‑150 (MQ=25)
                                                      cGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGAcaa                                                                                                >  8:276127/1‑148 (MQ=16)
                                                       gTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGAcca                                                                                                >  1:143688/1‑147 (MQ=39)
                                                          ggctcctgaaaCGACCCGAAGGGAGATGCTCGAAGCGATGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGtt                                                                                           <  8:33130/142‑1 (MQ=14)
                                                           ggacacgCAAAGAACCGAACGGTTAAGCTCTCAGCTCTGTGTAATCCTTCACCAAGCACGATCTATTTTAAAATTAGCTTTACTCTAAATCGTATGAGAAATTAAATATTTTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTc                                                                                          <  8:903403/145‑1 (MQ=9)
                                                              cAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTAGataggaagaccgg                                                                                       >  8:825452/1‑139 (MQ=11)
                                                                    aCGACCCGACGGGCGAAGCGCGAAGAGCTGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGc                                                                                 <  8:507614/151‑1 (MQ=255)
                                                                     ggaacggtgggggGAAGGGCGGAGAACGGAGGAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGc                                                                                 <  1:1029105/141‑1 (MQ=14)
                                                                     cGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTAGAGCg                                                                                >  2:769573/1‑151 (MQ=255)
                                                                         ccGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGTCac                                                                            >  3:206608/1‑150 (MQ=255)
                                                                             aGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCaagctc                                                                        >  1:927943/1‑146 (MQ=11)
                                                                             aGGGCGAAGCGCGCAGCGCATAGTAATACGTACCACAACACAAACACTTTCAAAAGAACATGAAATCTCAATCGAATCCGCAAAAAAAAAAAAAACACATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGCGCGGCAGCg                                                                          >  6:432963/1‑149 (MQ=2)
                                                                                   aaGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTAGAGCGGCAGCGAGAGAGCg                                                                  >  2:1041184/1‑151 (MQ=11)
                                                                                      cgctcAGCTCTTAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTg                                                               <  7:605503/151‑1 (MQ=255)
                                                                                        cttcgaaggttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTGACAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGcc                                                             >  6:999234/14‑151 (MQ=11)
                                                                                        cgccgcGCGGACGAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGcc                                                             <  3:335621/151‑1 (MQ=255)
                                                                                          cAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGc                                                           >  7:556836/1‑151 (MQ=255)
                                                                                          aaGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGCCAGCGAGAGAGCGTTGCCGc                                                           >  4:1249945/2‑151 (MQ=255)
                                                                                               ggtacgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                      <  5:515566/145‑1 (MQ=255)
                                                                                                gttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCaa                                                     <  3:690822/146‑1 (MQ=255)
                                                                                                     aaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAcc                                                >  7:110647/1‑151 (MQ=255)
                                                                                                        ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                             <  4:789221/151‑1 (MQ=255)
                                                                                                        ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                             <  7:530027/151‑1 (MQ=255)
                                                                                                           tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAg                                           <  2:438856/150‑1 (MQ=255)
                                                                                                           tCCCCCACCACCATCACTTTCAAAAGTAACAGAAAGCTCTAGCGAATCCGCAATGAAATATTATGCCCAAGCGGTTAAGGATGAGAAGATGCGAAGAGGGTTCGGCTGGAGCGCCAGCCAGACAGCGTTGCCGCAGGCCCCGACCCGAAgg                                          >  5:738581/1‑151 (MQ=2)
                                                                                                            cccccACCACCATCACTTTCAAAAGTCAAGGAACACTCAAGCGCCGCCGAAATCAAATAATCTGCCAATGCGGGCAAGGATGAGAAGCTTCTACGAAGATTCGTCTATAGCGGCAGCGAGCGAGAGTTGAAGCAGACAACGACCCGAAggg                                         >  6:1151154/1‑151 (MQ=2)
                                                                                                             ccccACCACCTTCAGTTTCAAAAGTCCCTGAAGTTTCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                        <  8:1032298/151‑1 (MQ=11)
                                                                                                                caccaccaTCATTTTAAAAAGTCCCTTAACTCTCAAGCTAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGTATTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGTGATAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                     <  8:800048/151‑1 (MQ=14)
                                                                                                                caccaccaTCACTTTCAACAGGCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGACAAGCTTCGACTAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                     <  7:200230/151‑1 (MQ=255)
                                                                                                                caccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                      >  7:169125/1‑150 (MQ=255)
                                                                                                                caccaccaTCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                     <  1:748417/151‑1 (MQ=255)
                                                                                                                caccaccaTCACTTTCAAAAGGCCCTGAACGCACAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                     <  3:541517/151‑1 (MQ=14)
                                                                                                                  ccaccaGCACTTACAAAAGGCCCTCAACTCCCAACCGAAGCCGCAATAAAATATTCGCCCCATGCGGGGAAGGATGAGAAGCTTCTACTACGGTGCGGCCCGAGCGCCCGCGAGAGAGCGTTGCCGCAGCCAACGACCCGAAGGGCGAAgc                                   <  5:242629/151‑1 (MQ=11)
                                                                                                                   cccctgccgatacAAAAGGCCCTGAACTCCCCACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTACGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                  <  1:815153/139‑1 (MQ=11)
                                                                                                                   cacctTCACTTTCAATAGTCCCTGAACTATCAATCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                  <  8:1071916/151‑1 (MQ=255)
                                                                                                                   caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                  >  5:1030460/1‑151 (MQ=255)
                                                                                                                   caccaTCACTTTCAAAAGTCCCTGAAATCTCAAACGAATCCACATATAAATATTCTCCGCATGCGGAGAAGTATGAGAAGATTGGAAGAAGGGTCTGCTCGCGCGAAAGGTATAGAGCGTTGCCGCGGGCGACGACCGGCAGGGCGAGgcg                                  >  2:669451/1‑151 (MQ=2)
                                                                                                                    accaTCACTTTCAAAAGCCCCTCACCTCCCACCCCAATCCGCAATAAAATATTCTGCCCATGCGGGTAAGTATGAGAAGCTTCGACTAAGGTGCGGCTCCCGCGCCCGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                 <  7:154488/151‑1 (MQ=11)
                                                                                                                    accaTCACTTTCAAAAGACCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGTCCGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                 <  5:783599/151‑1 (MQ=255)
                                                                                                                      caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGAAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                               <  2:817808/151‑1 (MQ=255)
                                                                                                                         cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                            >  8:580206/1‑151 (MQ=255)
                                                                                                                         cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                            >  6:150601/1‑151 (MQ=255)
                                                                                                                         cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                            >  5:358363/1‑151 (MQ=255)
                                                                                                                          aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGgagacgaaccgaccc                           >  5:345604/1‑140 (MQ=14)
                                                                                                                            tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                         >  1:1017532/1‑151 (MQ=255)
                                                                                                                            tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGACGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                         >  8:241037/1‑151 (MQ=255)
                                                                                                                             ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                        >  1:286852/1‑151 (MQ=255)
                                                                                                                              tCAAAAGTCCCTGAACTTTCAACCGAATCCTCATTTAAATTTTCTGACTTTGCGGGGAAGTATGTGAAGCTTGGAATATGGTTCGGCTGGAGCGTCAGCGTTATAGCGTTGCCGCATGCAACGACCCGAAGGGCGAAGCTCGCAGCGCTg                        <  6:33478/150‑1 (MQ=2)
                                                                                                                              tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                        >  2:1234523/1‑150 (MQ=255)
                                                                                                                              tCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                        <  5:1199746/150‑1 (MQ=255)
                                                                                                                               cAAAAGTCCCTTAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGa                       <  4:216200/150‑1 (MQ=255)
                                                                                                                               cAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGATCCGAAGGGCGAAGCGCGCAGCGCTg                        >  3:653766/1‑149 (MQ=255)
                                                                                                                               cAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                        >  6:227691/1‑149 (MQ=255)
                                                                                                                               cAAAAGTCCCTGAACTCGCAACCGAATCCCCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGa                       <  5:551470/150‑1 (MQ=255)
                                                                                                                               aaaacGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCGCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                      <  5:1000680/150‑1 (MQ=255)
                                                                                                                                aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                     >  6:270435/1‑151 (MQ=255)
                                                                                                                                aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                     >  2:619118/1‑151 (MQ=255)
                                                                                                                                aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                     >  8:105366/1‑151 (MQ=255)
                                                                                                                                aaaaGTCACTGAAATTTCAATCTAATCATAAATCAAATTTTCTGCACATGCGGGGAAGGATGAGAAGCTTCGACTAAGTTTGGTCTCGAGCTTCAGATAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGTCGAAGCGCGCAGCGCTGAGt                     <  4:677695/151‑1 (MQ=14)
                                                                                                                                 aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                     >  6:612975/1‑150 (MQ=255)
                                                                                                                                 aaaGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                    <  4:748124/151‑1 (MQ=255)
                                                                                                                                  aaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                    >  7:333420/1‑150 (MQ=255)
                                                                                                                                   aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                  >  5:401601/1‑151 (MQ=255)
                                                                                                                                   aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                  >  2:238705/1‑151 (MQ=255)
                                                                                                                                     gcccTCAACTCTCCCGCGAATCCCCAATCATATCTGCTGCCCATGCGGGGAAGGCTGATACGGATCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                <  1:897610/150‑1 (MQ=14)
                                                                                                                                       cgtcaaCTCGCAACCGAATCCGCAATAAAATATTGTGCCCTTGCGGGGAAGGATGAGAAGCTTCGACTACGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt              <  7:331239/147‑1 (MQ=11)
                                                                                                                                        cTGAATTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTTGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                      <  4:588720/142‑1 (MQ=255)
                                                                                                                                        cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAgagatcggac             >  3:589861/1‑142 (MQ=255)
                                                                                                                                          gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc           >  1:176318/1‑151 (MQ=255)
                                                                                                                                          gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcc            >  6:142072/1‑150 (MQ=255)
                                                                                                                                           aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc          >  7:329637/1‑151 (MQ=255)
                                                                                                                                           aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc          >  2:591650/1‑151 (MQ=255)
                                                                                                                                            aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc         >  5:1234780/1‑151 (MQ=255)
                                                                                                                                               ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc      >  1:629201/1‑151 (MQ=255)
                                                                                                                                                  aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc   >  4:572306/1‑151 (MQ=40)
                                                                                                                                                   aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc   >  1:183974/1‑150 (MQ=40)
                                                                                                                                                    gCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacca  >  2:147821/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCA  >  DH1/2591293‑2591590

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: