Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG28.0% 102.7 / 28.8 140D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/38);  new base G (8/25);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 3.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

AAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGA  >  DH1/2591369‑2591661
                                                                                                                                                |                                                                                                                                                    
aaGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTGGAGCGCCAGCg                                                                                                                                                 >  2:99419/1‑150 (MQ=255)
 aGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGGCCCGGAAGTCTCAAACAGATCGGAAGTCCAAAATTCTGCACAACCGTCCATGCGGGCGAAGCTTCGTGTAAGGTTCGGGTGGATCGtaatcg                                                                                                                                                 >  3:67487/1‑144 (MQ=2)
     cGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGACAGCgagag                                                                                                                                             >  1:663819/1‑149 (MQ=255)
         gcgcgcccacgggactAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGtt                                                                                                                                       <  5:883489/136‑1 (MQ=14)
         cgcctgcgccgggtctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGtt                                                                                                                                       <  5:166387/136‑1 (MQ=14)
             gcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTAGAGCGCCAGCGAGAGAGCGTTGCCg                                                                                                                                   >  8:464098/1‑151 (MQ=255)
              tcgatgtttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGc                                                                                                                                     <  6:1086176/137‑1 (MQ=255)
                   gtttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                             <  4:983032/145‑1 (MQ=255)
                      accgaTCGTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACg                                                                                                                          <  7:966702/147‑1 (MQ=255)
                       cgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCCACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                         <  5:470788/149‑1 (MQ=255)
                          aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                      >  2:183382/1‑151 (MQ=255)
                              ttCCCCCACCACCATCACTTTCAAAAGGCCCTGAAATCTCAAACGAATCCACAAGGGAATATACACCCGATGCGCTGAAGGAAGAGCCGATACGCGCAAGGTTCGGCTCGAGCGGGAGCTAGAGAGAGTTGCCGCCggggacgccgggggg                                                                                                                  >  1:973164/1‑138 (MQ=2)
                               ccatcaaccaccaTCACTTTCAAAAGTCCCTGAACTTCCAACCTAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                 <  2:139342/145‑1 (MQ=11)
                                cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                >  5:524613/1‑151 (MQ=255)
                                cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                <  4:791008/151‑1 (MQ=255)
                                cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGCTACTGCTCGAGCGTCAGCGAGAGAGCGTTGCGGGAGGCATCGCACCGAAcgg                                                                                                                >  4:175207/1‑151 (MQ=17)
                                 ccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               <  6:816103/151‑1 (MQ=255)
                                  cccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                              <  3:973118/151‑1 (MQ=255)
                                   ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                              >  4:923542/1‑150 (MQ=9)
                                        accaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCAGCCACTTCGGGGAAGGATGAGACGGTTCCACTATGGTTCTGCTAGAGCGGCAGCTAGATATCATTGAAGCAGGCAACGACCCGAAGGGCGAAgcgg                                                                                                        >  2:984766/1‑150 (MQ=14)
                                         ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                       >  2:87168/1‑151 (MQ=255)
                                         ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcc                                                                                                       >  8:84529/1‑150 (MQ=255)
                                          caTCACTTTCAAAAGTCCCTGAACTGTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      <  4:581965/151‑1 (MQ=255)
                                          caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  8:776584/1‑151 (MQ=255)
                                            tCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  8:688278/1‑149 (MQ=255)
                                             cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCCGCTCGAGCCTCAGAGAGAGAGCGTTGCCGCAGGCAACGACACGAAGGGCGAAGCGagcagc                                                                                                   >  2:914011/1‑151 (MQ=14)
                                              aCTTTCAAAAGTCCCTGAACTCTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                  <  4:418277/151‑1 (MQ=255)
                                                tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTATGCCCATGCGGGGAAGGATGAGAAGATTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgcg                                                                                                >  2:512992/1‑150 (MQ=255)
                                                 ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               >  8:433078/1‑151 (MQ=255)
                                                   cAAAAGTCCCTAAAATCTCATTTCAATCCTCTATGATATTTGGTTCTCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGa                                                                                              <  8:805029/150‑1 (MQ=11)
                                                     actggcccTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCGCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                                                                                           <  1:857140/146‑1 (MQ=11)
                                                     aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                            >  8:105681/1‑150 (MQ=255)
                                                      acgTCCCTGAACTCCCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGCGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTaa                                                                                          <  7:683626/149‑1 (MQ=255)
                                                       aggcccTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                         <  5:547208/148‑1 (MQ=255)
                                                       aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                         >  7:185865/1‑151 (MQ=255)
                                                        gTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGGATCGCCTGAGGTGACGGGGGAGCTGCAGGCAGAGAGCGTAGGCGCCGGCAACGACCCGAAGGGCGAAGAGCGCAGCGCTGAGTaggg                                                                                        >  7:75665/1‑148 (MQ=14)
                                                        gTCCCTCAACTCGCAAGCCAATCCGCAAGCAAATTTTCTGCCCATCCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATc                                                                                        <  1:567761/151‑1 (MQ=14)
                                                          cccTCAACTCCCAAGCCACCCCCCAATAAAATTTTCGCCCCTGGGGGGGACGTATGACACGCTTCCACAACGGTGCGCCGCGCGGGCCCCCCAGAGCGCGGGGCCGCCGCAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCt                                                                                      <  7:51665/151‑1 (MQ=2)
                                                            cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                    >  1:705391/1‑151 (MQ=255)
                                                             tGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcc                                                                                   >  5:828293/1‑151 (MQ=255)
                                                              caaCTCTCAAGCGAATCCGCAATCAAATATGCTGCCCATGCGGGGAAGGAGGAGAAGCTTCGACTAAGGTTCGGCGCGAGCGGCCGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                    <  5:197897/148‑1 (MQ=14)
                                                                 ctctCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc                                                                                >  1:955080/1‑150 (MQ=255)
                                                                  tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                         >  6:529078/1‑140 (MQ=255)
                                                                  aagctaatcatttctcAATCAAATATTCTGCCCATGGGGGTAAGGATGATAAGCTGCGACTAAGGTTCGCCTGGCGCGCCAGCTAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                               <  8:670894/135‑1 (MQ=11)
                                                                   ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                             >  1:351584/1‑151 (MQ=255)
                                                                   ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                             >  6:410652/1‑151 (MQ=255)
                                                                     cAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacca                                                                            <  6:1071774/150‑1 (MQ=12)
                                                                      aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                          >  6:686806/1‑151 (MQ=40)
                                                                         cGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTCGACGGTACCCCACACCGCGTAGCACCAGGCACCCCCGCGAAGGGGGCAGCGCGGCGCGCGGAGTAATCCTACCCCCACCAGCAtg                                                                       >  7:1015563/1‑150 (MQ=2)
                                                                          gAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCa                                                                      <  6:444652/151‑1 (MQ=14)
                                                                          gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcc                                                                      >  6:950294/1‑150 (MQ=38)
                                                                            aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCt                                                                    >  5:248980/1‑151 (MQ=38)
                                                                            aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCt                                                                    >  3:139999/1‑151 (MQ=38)
                                                                            aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCt                                                                    >  5:327063/1‑151 (MQ=38)
                                                                                gCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttgc                                                                 >  3:678394/1‑148 (MQ=255)
                                                                                  cttCATATCTTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCa                                                                <  4:988994/147‑1 (MQ=255)
                                                                                  aaTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCat                                                               >  5:868018/1‑149 (MQ=255)
                                                                                   aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCaa                                                               <  4:891540/149‑1 (MQ=25)
                                                                                   aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCaaa                                                              >  7:607430/1‑150 (MQ=38)
                                                                                     cAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccaccaagaatttaagacgc                                                           >  3:1040151/1‑137 (MQ=11)
                                                                                     aaaaTTTTGTGCCGTTGGGGGGAAGTTTGTGACGCTTCGACTTTGGTTGGCGTCGCGCGGCCGGGAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAg                                                            <  1:1062675/149‑1 (MQ=2)
                                                                                     aaaaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAg                                                            <  1:803834/149‑1 (MQ=18)
                                                                                      aaaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGCCCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGcc                                                          >  1:966481/1‑151 (MQ=37)
                                                                                             cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAt                                                   >  1:266537/1‑151 (MQ=40)
                                                                                             atgtccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAt                                                   <  2:1520/147‑1 (MQ=255)
                                                                                             atctccATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                    <  5:158637/146‑1 (MQ=255)
                                                                                               gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCataaagacacaaatcc                                                 >  5:142556/1‑139 (MQ=11)
                                                                                               ccgatgtcGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtc                                                 <  1:375233/144‑1 (MQ=255)
                                                                                                 ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCataatcaatcaaaacctg                                               >  1:1006051/1‑137 (MQ=11)
                                                                                                  cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGc                                                                                                                                     >  3:405079/1‑64 (MQ=255)
                                                                                                  cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCATAAgtaaacaaaagctg                                               >  5:79762/1‑137 (MQ=11)
                                                                                                  cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                       >  4:725925/1‑126 (MQ=37)
                                                                                                     gCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaaa                                            <  7:718763/150‑1 (MQ=12)
                                                                                                     gCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaaa                                            <  1:267992/150‑1 (MQ=255)
                                                                                                         ggAAGGATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGGTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGaa                                        <  4:548261/150‑1 (MQ=255)
                                                                                                          gAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATc                                      >  2:514314/1‑151 (MQ=255)
                                                                                                           aaGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATcc                                     <  4:1007678/151‑1 (MQ=255)
                                                                                                            aGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCg                                    <  8:363809/151‑1 (MQ=255)
                                                                                                            aGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCg                                    <  5:1055726/151‑1 (MQ=255)
                                                                                                              gATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCa                                  <  8:892622/151‑1 (MQ=255)
                                                                                                               aTGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCaa                                 <  4:50861/151‑1 (MQ=255)
                                                                                                               aTGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCa                                  <  4:1020615/150‑1 (MQ=255)
                                                                                                                tGAGAAGCTTCGACTAAGGTTCGGGTCGAGCGGCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCaa                                 <  8:456704/150‑1 (MQ=255)
                                                                                                                tGAGAAGCTGCGACTAGGGTTCGGCTCGCGCGGCGGGTTGTGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCaa                                 <  7:676062/150‑1 (MQ=14)
                                                                                                                ggagaAGCTGCGACTACGGGTCGGCTCGAGCGTCAGCTATAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAAt                                <  1:513409/150‑1 (MQ=255)
                                                                                                                 gagaAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCAAAAGCCAACGAATGCTCTAACgccgagcgaatt                               >  7:805589/1‑140 (MQ=14)
                                                                                                                  agaAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                              <  5:635472/151‑1 (MQ=255)
                                                                                                                    aaGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCaaa                            <  6:880147/151‑1 (MQ=255)
                                                                                                                     aGCTTCGACCAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAat                           <  1:755542/151‑1 (MQ=255)
                                                                                                                      tcttatattAAGGTTCGTCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCaa                             <  6:718386/140‑1 (MQ=255)
                                                                                                                        ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACAAAACACCCGCACATCGGAAGCGAAACCGCAAGCacaaaat                        >  3:119235/1‑145 (MQ=11)
                                                                                                                            aCTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACTTAACGACACAAAAGCTCAAACGAATCCGCAATCAATTAGTCTGc                    >  5:1043861/1‑151 (MQ=14)
                                                                                                                               aaGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGccc                  <  2:54151/150‑1 (MQ=255)
                                                                                                                               aaGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                 <  4:818615/151‑1 (MQ=255)
                                                                                                                                 ggTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAg               <  4:204117/151‑1 (MQ=255)
                                                                                                                                   ttCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCATAATACACCGAATACTGAAACGCATCCCGAATAAACAATTCTGccaccccg             >  7:718028/1‑145 (MQ=14)
                                                                                                                                     cGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcctta                                                              >  3:568722/1‑95 (MQ=2)
                                                                                                                                          cGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGaa        <  2:266306/149‑1 (MQ=255)
                                                                                                                                            aGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCTAAACCAAACCAATACTCATACTAATCGCCAATGATAGACTCTGGCCCCGCGGGTAAGGAt    >  5:378996/1‑151 (MQ=2)
                                                                                                                                              gctCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCTGAATTCTCAACCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATga  <  1:102770/149‑1 (MQ=21)
                                                                                                                                                |                                                                                                                                                    
AAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGA  >  DH1/2591369‑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: ‑