Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG26.5% 202.4 / 46.8 250D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (88/59);  new base G (26/27);  total (130/120)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

CGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGA  >  DH1/2591367‑2591661
                                                                                                                                                  |                                                                                                                                                    
cGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTGGAGCGCCAGc                                                                                                                                                  >  2:1475348/1‑151 (MQ=255)
              cgcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGcc                                                                                                                                    >  3:545100/1‑151 (MQ=255)
                 cgaccggtctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCa                                                                                                                                 <  1:1105929/142‑1 (MQ=255)
                 aGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGTCTAGAGCGGCAGCGAGAGAGCGTTGCCGc                                                                                                                                  >  4:26882/1‑150 (MQ=11)
                     ggttctAATCCTTCCCCCACCACCATCTCTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                             <  2:1112974/146‑1 (MQ=11)
                       gtcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAAc                                                                                                                           <  3:992426/147‑1 (MQ=255)
                         gTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACg                                                                                                                          >  7:29584/1‑150 (MQ=255)
                           aaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTATCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAc                                                                                                                        <  4:1515506/150‑1 (MQ=255)
                                 gcccccACCACCATCACTTTCAAAAGGCCCTCAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCGCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                 <  5:475403/150‑1 (MQ=255)
                                  cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAATCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                 <  6:540552/150‑1 (MQ=255)
                                  cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCAGCgagag                                                                                                                                             >  5:127129/1‑122 (MQ=255)
                                   ccctacccagctCACTTTCAAAAGTCCCTGAACTACCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               <  3:16609/140‑1 (MQ=11)
                                   ccccACCACCATCACTTTCAAAAGTCCCTGAACTATCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               <  4:882917/151‑1 (MQ=255)
                                   ccccACCACCATCACTTTCAAAAGGCCCTGAACTCTCAAGCGAAACCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTACGGGTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               <  1:280358/151‑1 (MQ=255)
                                     ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                              >  5:739092/1‑150 (MQ=255)
                                      tcctatcttCACTTTAAAAAGTCCCTTAACTTTCAAGCGAATCCGAAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCTACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                             <  4:1637741/142‑1 (MQ=14)
                                      caccaccaTCACTTACAAACGGCCCTGAACTCACAACCCAAGCCCCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCGCGAGCGCCCGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGaa                                                                                                            <  3:285335/151‑1 (MQ=14)
                                       accaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                             >  5:1285988/1‑149 (MQ=255)
                                       accaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                             >  2:61228/1‑149 (MQ=255)
                                        ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                          >  7:595368/1‑151 (MQ=255)
                                        ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                          >  3:1038263/1‑151 (MQ=255)
                                         caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                         >  5:1680871/1‑151 (MQ=255)
                                           ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                       >  7:1093035/1‑151 (MQ=255)
                                            caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  4:1354845/1‑151 (MQ=255)
                                            caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  3:213511/1‑151 (MQ=255)
                                             atcacatgcaccactCCCTCAACTCTCAAGCGCACGCGCCTTCAAACCTGCGCCCCTGGCGGGTAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                     <  3:1544190/137‑1 (MQ=2)
                                               aaCTTTCAAAAGTCCCTCAACTCCCAAGCGAATCCGCAATAAAATATTCTGCCAATGCGGGGAAGCATGAGAAGCTTCTACTAAGGTTCGGCTCGCCCGCCACCCAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                   <  1:1049215/150‑1 (MQ=11)
                                                  tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                >  6:1228859/1‑151 (MQ=255)
                                                   ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               >  3:1288443/1‑151 (MQ=255)
                                                   ttCAAAAGGCCCTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGCCGCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                <  3:261556/150‑1 (MQ=11)
                                                     cAAATGTCCCTGAACTTTCATCCGAATCCGCAATTAAATATTGTGCCAATGGGGGGAAGGATGAGAAGGTTCGACTAAGGTTGGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGa                                                                                              <  4:301639/150‑1 (MQ=11)
                                                         aggcccTGAACTCGCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                         <  3:1470043/148‑1 (MQ=255)
                                                             ccTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                     >  7:1028595/1‑151 (MQ=255)
                                                             ccTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                     <  3:968970/151‑1 (MQ=255)
                                                              cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                    >  4:1462688/1‑151 (MQ=255)
                                                               tGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcc                                                                                   >  4:731609/1‑151 (MQ=255)
                                                                  aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc                                                                                >  6:20373/1‑151 (MQ=255)
                                                                    tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                              >  6:1462473/1‑151 (MQ=255)
                                                                    tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                              >  3:429510/1‑151 (MQ=255)
                                                                    tctcAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                              >  2:1250254/1‑151 (MQ=16)
                                                                     cccAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGTATGATAAGTTTCTACTATGTTTCGGCTCGAGCGTCTGCGATATAGCGTTGCCCCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                               <  2:823518/147‑1 (MQ=11)
                                                                       cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGACTaacggaaccgccccccg                                                                           >  7:401509/1‑136 (MQ=14)
                                                                        aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                          >  8:443168/1‑151 (MQ=40)
                                                                        aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                          >  3:1452996/1‑151 (MQ=40)
                                                                         aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                            >  1:929713/1‑132 (MQ=25)
                                                                          gCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                        >  6:1583615/1‑151 (MQ=39)
                                                                            gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcg                                                                      >  3:1351746/1‑150 (MQ=38)
                                                                            gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcc                                                                      >  5:107029/1‑150 (MQ=38)
                                                                              aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacagatccg                                                                    >  1:935512/1‑144 (MQ=255)
                                                                              aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCt                                                                    >  8:208152/1‑151 (MQ=38)
                                                                               tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAGGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                   >  5:1343906/1‑151 (MQ=255)
                                                                               tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCtt                                                                   >  8:1624506/1‑151 (MQ=38)
                                                                                ccGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAGGCTTCGACTAAGGTTCGCCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                   <  6:1343240/150‑1 (MQ=255)
                                                                                ccGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttg                                                                  >  7:74142/1‑150 (MQ=38)
                                                                                 cGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcc                                                                 >  7:1552391/1‑149 (MQ=12)
                                                                                   cAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCaa                                                               <  4:670622/151‑1 (MQ=18)
                                                                                   cAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCa                                                                >  6:435697/1‑150 (MQ=38)
                                                                                       cAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAGCAcggaacagccc                                                           >  5:395550/1‑141 (MQ=255)
                                                                                       cAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                          <  6:401962/136‑1 (MQ=11)
                                                                                       aaaaTATTCTGCCCATGCGGGGAAGGATGAGAAGGTTCGACTAAGGTGCGGCTCGCGCGCCAGCGAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGc                                                           <  7:181589/150‑1 (MQ=255)
                                                                                         aaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttccttaaata                                                          >  2:863807/1‑141 (MQ=9)
                                                                                           tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCtaaacagaac                                                       >  3:1088230/1‑145 (MQ=255)
                                                                                           tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttacttaaaaaaac                                                       >  2:183928/1‑139 (MQ=255)
                                                                                            atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCg                                                       <  3:909928/150‑1 (MQ=39)
                                                                                            atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGACAGAGCGTTGCCGCAGGCAACGACCCGAAGCGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtatataccaagagcc                                                      >  4:577271/1‑137 (MQ=11)
                                                                                             ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTGCTAAAgaacaca                                                      >  5:238561/1‑144 (MQ=255)
                                                                                              tCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttccttaaaaaatcaa                                                     >  8:29425/1‑136 (MQ=11)
                                                                                              tCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGa                                                     <  2:651152/150‑1 (MQ=38)
                                                                                                tGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCataaataaacaaaa                                                   >  8:1273419/1‑140 (MQ=11)
                                                                                                 gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAtt                                                  <  8:1054765/150‑1 (MQ=255)
                                                                                                 gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCAAAAGcaaacgaatac                                                 >  7:942373/1‑141 (MQ=255)
                                                                                                 ccgatgtcGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAtt                                                  <  6:568032/143‑1 (MQ=255)
                                                                                                   ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                               >  1:643429/1‑119 (MQ=38)
                                                                                                    cctgcgctgatggATGTGAAGTTTCTACTAAGTCTCGGTTGGCGGGTTAGCGAGAGTGCTTTGGCAGAGGAATCGACCCGACGGGGGAAGGTCGTAGCTGTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtct                                                <  4:1242407/138‑1 (MQ=9)
                                                                                                      tGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                   <  1:150623/96‑1 (MQ=25)
                                                                                                       ggggggAAGGATGAGACGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaaa                                            <  2:253152/148‑1 (MQ=255)
                                                                                                       gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                           <  8:120285/151‑1 (MQ=255)
                                                                                                       gCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                           <  5:605179/151‑1 (MQ=255)
                                                                                                       ccggtcttggtttataaGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                           <  8:101522/136‑1 (MQ=14)
                                                                                                        gggggAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTT‑CCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGa                                         <  6:1087098/150‑1 (MQ=255)
                                                                                                         ggggAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               >  2:1164151/1‑81 (MQ=40)
                                                                                                           ggAAGGATGAGAAGCTTCGACTAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAAt                                       <  5:1381943/151‑1 (MQ=25)
                                                                                                            gAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGATTGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATc                                      <  2:1333331/151‑1 (MQ=255)
                                                                                                             aaGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATcc                                     <  1:25451/151‑1 (MQ=255)
                                                                                                              aGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCg                                    <  1:837019/151‑1 (MQ=255)
                                                                                                                        gCTTCGACAAAGGTGCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAata                          <  5:1130375/151‑1 (MQ=255)
                                                                                                                          ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACAAAAGACCACAAAAACTCAAACGCATCCGAAAtaataaaa                         >  6:792043/1‑143 (MQ=14)
                                                                                                                            tgACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCAAAAGAAACGCAATGCTCAAACGAATACGCAAtatataaatct                      >  1:116625/2‑141 (MQ=14)
                                                                                                                            cgacgaAGGTTCGACTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCt                      <  1:68754/151‑1 (MQ=255)
                                                                                                                             gACTAAGGTTCGTCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                     <  6:680983/151‑1 (MQ=255)
                                                                                                                              aCTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACAAAAGGACAGGAAATCTCAAACGAATCCGGAATaaaagaatatgc                    >  4:831294/1‑143 (MQ=14)
                                                                                                                                tcagttgcgacTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                     <  4:561610/138‑1 (MQ=255)
                                                                                                                                 aaGGGTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                 <  8:1521481/151‑1 (MQ=255)
                                                                                                                                  aGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCaa                <  3:1562051/151‑1 (MQ=255)
                                                                                                                                   ggTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                 >  1:790305/1‑149 (MQ=255)
                                                                                                                                   ggTTCGACTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAg               <  2:12555/151‑1 (MQ=255)
                                                                                                                                    gtgcgccTCGCGCGCCAGCTAGAGAGGGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGc              <  7:1207554/145‑1 (MQ=14)
                                                                                                                                    gTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGc              <  1:1310353/151‑1 (MQ=255)
                                                                                                                                          cgcGAGCGCCAGGGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGaa        <  1:1265127/149‑1 (MQ=255)
                                                                                                                                                cGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCGCCAATTCTCATACGAAGCCGGAATGTTATAGTCTGCCACCGCGGGTAAGGaggc  >  8:1428184/1‑148 (MQ=14)
                                                                                                                                                  |                                                                                                                                                    
CGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGA  >  DH1/2591367‑2591661

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