Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG27.1% 269.1 / 59.1 297D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (114/64);  new base G (21/45);  total (163/134)
Fisher's exact test for biased strand distribution p-value = 1.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

GACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591363‑2591660
                                                                                                                                                      |                                                                                                                                                    
    gtatgcgcgaacgccgcAGCGCGGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGc                                                                                                                                                 <  7:1825924/137‑1 (MQ=17)
                    cAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGc                                                                                                                                 >  3:1553909/1‑151 (MQ=255)
                     aGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTTGAGCGCCAGCgag                                                                                                                                              >  3:769159/1‑137 (MQ=255)
                            tggaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGc                                                                                                                                    <  4:1253676/137‑1 (MQ=255)
                            aGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACg                                                                                                                         <  6:1598479/151‑1 (MQ=255)
                              tAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGAc                                                                                                                       >  7:969527/1‑151 (MQ=255)
                                  ccTTCCCCCACCACCATCACTTTCAAAACTAAGGGAAAGCTCTAGCTAATCCGGAATCAAGTATTCTCACAATTCTGGGAACGATGAGAAGAGTCGACTAACGTTGGCCTGCATCGGCATCGAAACACCGTTGCCGCCGGcgccgccgcgg                                                                                                                   >  5:1951803/1‑141 (MQ=2)
                                   cTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                   >  3:1137702/1‑150 (MQ=255)
                                     tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAg                                                                                                                 <  6:623710/150‑1 (MQ=255)
                                      cccccACCACCATCACTTTCAAAAGTCCCTGAACTTTCAACCGAATCCGAAATAAAATATTCTGCCCATGCGGGGAAGGATGATAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                               <  8:1410747/151‑1 (MQ=255)
                                      cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                               <  5:1596525/151‑1 (MQ=255)
                                       ccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              >  4:1395820/1‑151 (MQ=255)
                                       ccccACCACCATCACTTTCAAAAGTCCCTGAACTCCCCACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              <  3:17874/151‑1 (MQ=255)
                                        cccACCACCATCACTTTCAAAAGTAACGGAACTCTCAAGCGAATCCTCAATCAAATATTCTTCCCAAGCGCGGAAGTATGAGAAGCTTCGACTAAAGTTCGACTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCCACGCCCCGAAGGGCg                                                                                                             >  1:1020075/1‑151 (MQ=11)
                                          caccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                            >  8:169588/1‑150 (MQ=255)
                                           cccaccaccaCTTTCATCCGTCCTTTAACTCCCAACCTAATCCGAAATAAAATTTTCTGCCCATGGGGGTAAGGATGATAAGCTTCGACTATGGTTCGGCTCGAGCGCCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAg                                                                                                          <  6:448143/150‑1 (MQ=11)
                                           accaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGTGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                            >  8:1604133/1‑149 (MQ=255)
                                           accaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                            >  3:968306/1‑149 (MQ=255)
                                            ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTAGAGCGACAGCGAGAGAGCGTTGCCTCAGGCAACGACCCGAAGGGCGAAgc                                                                                                         >  4:552215/1‑151 (MQ=255)
                                             caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                        >  1:1873562/1‑151 (MQ=255)
                                             aaccaTCACTTTCAAAAGTCCCTGAACTTACAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                        <  6:1514464/150‑1 (MQ=255)
                                              accaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGagctccgaagagc                                                                                                       >  1:922530/1‑139 (MQ=255)
                                                  ccACTTTCAAACGCCCCGCAACTCCCACGCCAACCCGCAAGAAAATATGCTGCCCATGCGGGGATGGATGAGAAGCTTCGACTACGGTGCGGCGCGCGCGGCAGCGAGAGAGCGTGCCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAg                                                                                                   <  5:1799416/150‑1 (MQ=9)
                                                   cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                  <  4:1848511/151‑1 (MQ=255)
                                                    aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 >  4:926906/1‑151 (MQ=255)
                                                    aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 >  7:356024/1‑151 (MQ=255)
                                                    aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 >  4:385081/1‑151 (MQ=255)
                                                     cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTGGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                                                                                       >  2:1795678/1‑144 (MQ=255)
                                                       ttcaaaccgcccTGACCTCACAAGCGAATCCGCATGCAAATATGCGGCCCCGGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                               <  3:903340/141‑1 (MQ=14)
                                                       ttCAAAAGTCCCTGAAGTTTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                               <  2:1040085/150‑1 (MQ=255)
                                                       ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                               >  6:301577/1‑150 (MQ=255)
                                                       atcacacgtaccTGAACGCTCATGCAAATCCGCAGTGACATCTTCGGCCCTTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                <  5:217433/139‑1 (MQ=14)
                                                        tCAAACGGCCCTCAACTCTCAACCGAACCCGCAATCAAATATTCTGCCCCTGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                              <  3:384429/150‑1 (MQ=11)
                                                         cAACAGTCCCTGAACTCCCCACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                            <  7:1697236/151‑1 (MQ=255)
                                                          atacgctcttcctatctaAAGCTAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           <  2:220889/137‑1 (MQ=11)
                                                           aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           >  7:1005179/1‑150 (MQ=255)
                                                            aaGTCCCTGAACTCTCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGTGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTaa                                                                                         <  2:1470075/151‑1 (MQ=255)
                                                             aGTCCCTGAACTGTCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                        <  2:593458/151‑1 (MQ=255)
                                                              gTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATc                                                                                       >  4:234011/1‑151 (MQ=255)
                                                                 ccTGAACTCTCAAGCTAATCCTCACTCAAATATTCTGACCATGCGTTGAAGGATGAGAATCTTCGACTACGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                    <  6:1680178/151‑1 (MQ=17)
                                                                 ccTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTTGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGAGCtcggaagagca                                                                                     >  5:865429/1‑140 (MQ=11)
                                                                   tGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcc                                                                                  >  7:554335/1‑151 (MQ=255)
                                                                    gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  6:374472/1‑151 (MQ=255)
                                                                    gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  3:1284796/1‑151 (MQ=255)
                                                                    gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  4:474962/1‑151 (MQ=255)
                                                                     aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc                                                                                >  7:1245452/1‑151 (MQ=255)
                                                                       ctctCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                   >  4:1205771/1‑146 (MQ=255)
                                                                        tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                             >  8:211220/1‑151 (MQ=255)
                                                                        tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                             >  5:1691341/1‑151 (MQ=255)
                                                                         ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                            >  3:1076915/1‑151 (MQ=255)
                                                                           gatattaatttctaaTCAAATATTGTGCCGATGGGGGGAAGGATGAGATGCTTCGACTAAGGTTCGTCTGGAGCGGCTGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                          <  8:1989255/138‑1 (MQ=11)
                                                                           ccatcGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                          <  2:1186850/147‑1 (MQ=9)
                                                                           cAAGCGATTCCGCATTTAAATATTCTGCCCTTGCGGGGATGTATGAGAAGGTTCGACTAAGGTTCGCCTGGTGCGTCAGCGAGAGAGCTTTGCGGGAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacca                                                                           <  8:648012/150‑1 (MQ=11)
                                                                            aatcGAATCCTCAATCAATTTTTCTGCCCATGGGGGGATGGATGAGAAGCTTCGACTAAGGTTGGTCTTGAGGTCTAGGGAGTGAGGGTTGCGGCAGGCAAGGACCCGAAGGGGGAAGCGCGCAGCGCTGAGTATTCCTTGCCggaccac                                                                          <  6:14324/147‑1 (MQ=2)
                                                                              agcaacactaattaaaaTATTCTTCCCATCCTCGTAGGTATTCGTAGCTTCGACTAAGGTTCGTCTGGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                       <  2:1969077/137‑1 (MQ=11)
                                                                                  aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                        >  1:1428431/1‑114 (MQ=40)
                                                                                   tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                  >  7:1664639/1‑151 (MQ=38)
                                                                                   tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCACccaccacc                                                                         >  6:823/1‑144 (MQ=38)
                                                                                   tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCtt                                                                  >  4:891967/1‑151 (MQ=255)
                                                                                    ccGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttc                                                                 >  4:734951/1‑150 (MQ=37)
                                                                                      gttctCATATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCa                                                               <  3:459870/147‑1 (MQ=255)
                                                                                      gCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCCCCATCGCttcct                                                               >  1:943040/1‑148 (MQ=255)
                                                                                        aaTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCat                                                              >  2:1534448/1‑149 (MQ=255)
                                                                                          tCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttccttaa                                                            >  2:676497/1‑144 (MQ=16)
                                                                                          tCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACaaaa                                                            >  6:247654/1‑150 (MQ=255)
                                                                                           cAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttgcttaaat                                                          >  3:134983/1‑143 (MQ=255)
                                                                                              atatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGcccc                                                       <  1:538762/151‑1 (MQ=25)
                                                                                                atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCagatcgcaagc                                                     >  7:645964/1‑141 (MQ=255)
                                                                                                atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCataaaaaaac                                                      >  2:1634684/1‑144 (MQ=255)
                                                                                                  tCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttccttaaaaaatc                                                      >  5:907215/1‑136 (MQ=255)
                                                                                                    tGCCCAAGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACaaaaaaaatcaaaa                                                  >  2:760026/1‑140 (MQ=11)
                                                                                                     gCCCATGCGGGTAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtc                                                <  6:1283309/151‑1 (MQ=255)
                                                                                                     gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                            >  5:1105209/1‑123 (MQ=37)
                                                                                                     gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACtaaaaaaaacaaaaac                                                >  6:1760696/1‑139 (MQ=11)
                                                                                                     gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                       >  3:830624/1‑128 (MQ=21)
                                                                                                     gCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTGCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtc                                                <  5:1459832/151‑1 (MQ=255)
                                                                                                      cccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTCATCCTTCCCccacc                                                                            >  1:2059922/1‑122 (MQ=37)
                                                                                                      cccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtct                                               <  7:499656/151‑1 (MQ=255)
                                                                                                      cccATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtct                                               <  5:671234/151‑1 (MQ=18)
                                                                                                      cccATGCGGGGAAGGATGAGAAGCTTCCACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtct                                               <  5:508273/151‑1 (MQ=9)
                                                                                                       ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCa                                                                                                                                <  6:1065163/69‑1 (MQ=255)
                                                                                                         aTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGAACCGGAATTCTCa                                             >  4:1712980/1‑150 (MQ=255)
                                                                                                           gCGGGTAAGGATTATAAGCTTCTACAAAGGTTCGGCTCGAGCGCCAGCGATAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                          <  8:1078848/151‑1 (MQ=11)
                                                                                                           gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                          <  3:250992/151‑1 (MQ=255)
                                                                                                           gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                          <  2:907370/151‑1 (MQ=255)
                                                                                                             ggggAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGa                                        <  2:308363/151‑1 (MQ=255)
                                                                                                                  aGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCCCCATCCCTTTCATAAGCACAGCAATCCTCTAACGCAGCCg                                   >  1:39986/1‑151 (MQ=14)
                                                                                                                        agaAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGGGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATc                              <  8:887749/150‑1 (MQ=255)
                                                                                                                        agaAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCCCCACCAGCCCTTTCAAAAGCCCCGGAATTCTCAAACGCATCCGCAATc                              >  3:559580/1‑150 (MQ=255)
                                                                                                                             cgacgtgtccggTTCGGCTCTAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAatat                        <  7:1610169/141‑1 (MQ=255)
                                                                                                                              ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCTAAAGCAAACCAAAACTCAAACGAAGCCGCAATGaactcgt                       >  1:1688444/1‑146 (MQ=14)
                                                                                                                               tCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCCCCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATAtt                       >  1:1483741/1‑150 (MQ=255)
                                                                                                                               gcGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTc                      <  7:138335/150‑1 (MQ=255)
                                                                                                                                 gACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                    <  2:479584/151‑1 (MQ=255)
                                                                                                                                     aaGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAAAGCGCAGCGCTGAGTAATCCTTACCCCACCAACAAGAATTTAAAAAGCCCCGGAAATCTCAAACGAATAAGAAATCAAATATTCTGccc                 >  2:1335578/1‑150 (MQ=14)
                                                                                                                                       ggTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                >  7:664791/1‑149 (MQ=255)
                                                                                                                                             gCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCgggg         >  3:1972673/1‑150 (MQ=255)
                                                                                                                                                cGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAg      <  2:414989/150‑1 (MQ=255)
                                                                                                                                                cGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAgg     <  1:1246331/151‑1 (MQ=38)
                                                                                                                                                 gAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCTAAAGCACCGGAAAGCTCTAACGCCGCCGGAATGATAGATTCTGCCCCAGCGGTGaacgg    >  1:7618/1‑148 (MQ=14)
                                                                                                                                                      |                                                                                                                                                    
GACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591363‑2591660

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: