Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,8090TC75.0% 124.5 / 52.6 96T29A (ACA→GCA) borDprophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base T (9/14);  new base C (32/37);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATG  >  NC_000913/578562‑579045
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
cAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGGTTCCTTTGGTGCTACTGctgtctg                                                                                                                                                                                                                                            <  1:330480/251‑1 (MQ=255)
 aGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGTTGCTACTGctgtctg‑t                                                                                                                                                                                                                                          >  1:346719/1‑251 (MQ=255)
     agagCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCtg‑tttgt                                                                                                                                                                                                                                      <  1:246569/251‑1 (MQ=255)
        gttgcccATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGtttt                                                                                                                                                                                                                                   >  1:71946/6‑251 (MQ=255)
        gCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACCCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTTCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGtttt                                                                                                                                                                                                                                   <  2:355053/251‑1 (MQ=255)
           gtacatagcgaggggccatcaaTGCAATATTTTTGAGCAATACACACGCGCTTACAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGTAATCCATTTACGAATGTTTGCGGGGTTTCTGTTTT‑ACCAACATTTTTTGCGCCGCCCCAAATTTTAGCTGCATCGGCTGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATTGGTGATGTTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGaa                                                                                                                                                                                                                                <  2:382255/230‑1 (MQ=255)
               atatCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGCATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGt                                                                                                                                                                                                                            >  1:411008/1‑251 (MQ=255)
                 atCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTCCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTa                                                                                                                                                                                                                           <  1:131902/250‑1 (MQ=255)
                   cGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTaaac                                                                                                                                                                                                                        >  1:439443/1‑250 (MQ=255)
                         gCAACTCATGCAATTATTGTGAACAATACACACGCGCTTCCAGCGGCGTATAAATACCTAAAGTAATAAAACCGCGCAATCCATTTACGACTGTTTTCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTTCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGTTGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAGCAGTAAACGTCttt                                                                                                                                                                                                                  >  2:428243/1‑249 (MQ=255)
                           aaCTCATGCCATTATTTTGCGAAATACACACGCGCCTCCAGCGGAGTCTAACCGCCTAAAGTAATAACCCCGAGCAATCCATTTACGAGTGTTTGCGGGGTTTCTGTTTT‑ACCAACATTTTCTGAGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGCAGAAAT‑GATGGGTGTTGGTTTCCTTTGGTGCTTCTGCTGTCTG‑TTTGTTTTGAACAGTAATTGTCTGttt                                                                                                                                                                                                                >  2:451266/1‑250 (MQ=255)
                             cTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTg                                                                                                                                                                                                                >  1:418547/1‑249 (MQ=255)
                                    aaTTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCACTTCCAGAAACGAAGAAAT‑GATGGTTGATGGTTTCCTTTGGTGCTACTGCTGCCGT‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATcc                                                                                                                                                                                                       >  2:4792/1‑251 (MQ=255)
                                          tgtGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGACTGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGAGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACCTACTGtaat                                                                                                                                                                                                 >  2:312373/3‑251 (MQ=255)
                                             gAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAg                                                                                                                                                                                              <  1:117971/251‑1 (MQ=255)
                                              aGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGc                                                                                                                                                                                             <  1:184126/251‑1 (MQ=255)
                                                    acacacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGTGTGATGGTTTCCTTTGGTGCTACTGCTGCCgg‑tttgt                                                                                                                                                                                                                                      <  1:356267/204‑1 (MQ=255)
                                                    acacacGAGCTTCCAGCGGAGTATAAATGCCTTAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGACACAAATTTTGGCTGCATCGAAAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCATTTGGTGCTACTGATGCCgg‑tttgt                                                                                                                                                                                                                                      >  2:356267/1‑204 (MQ=255)
                                                     cccactcgtttCCAGCGGGTTATAAATGCCTAATGTAATAAAACCGACGACTCCTTTTACGAATGTTTGCTGGGTTTCTTTTTT‑AACATCATTTTCTGCGCCGCTACAAATTTTAGCTGCATCGTCAGTTTTCTTCTGCCCAATTCCAGATACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCCGGGCCa                                                                                                                                                                                      <  2:390001/242‑1 (MQ=255)
                                                     cacccgctcttccatggtaGTAAAAATTCCTAAAGTAATAAAACCGTGCAATCCATTTACGAATGTTTGATGGGTTTCTGTTTT‑AAAACCATTTTCAGCGCCGCCACAAATGTTATCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAGACGAAGAAAT‑GATGGGTGATGGTTGACTTTGGTGCTGCTGCTGTCTG‑TTTGTTATGAACAGTAAAGGTCTGTTGAGCACTTCCTGTAATAAGCAGGGCCa                                                                                                                                                                                      <  2:323021/233‑1 (MQ=255)
                                                      acacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGAGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGcc                                                                                                                                                                                       >  1:229617/1‑249 (MQ=255)
                                                       cacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCCGGGCCAgc                                                                                                                                                                                    >  2:110172/1‑251 (MQ=255)
                                                               tCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                            >  1:250253/1‑251 (MQ=255)
                                                                ccAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                            <  1:318392/250‑1 (MQ=255)
                                                                 taGCGGAGAATTAATGACGAAAGTAATAAAACCGAGCAATCCATTTACGAACGTTTGCTGGGTTTCTGTTTT‑AACAAAATTTTCTGCGCCGACACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGCAAAGAAGAAAT‑GATGGGAGGTGGTTTCTTTTGTTGCGACTGCTGCCGG‑TTGGTTTTGAGTAGTAAACGGCTGTTGTGCACATCCTGTAATAAGCTGGGCCAGAGCAGTAgat                                                                                                                                                                           >  2:14814/2‑248 (MQ=255)
                                                                 aagctgaatataAATGCCTAAAGTAATAAAACCGAGCAATCAATTTACGAACGATTGATGGGATTCTATATC‑AACAACAATTCCTGCGCCGCCACAAATTTTGGCAGCCTCGACAGTCATCTTCTGCCCCATTCCAGCAACGATGAAATTGA‑CAATGATGATATCATTTGTTGCTACTGCTGACGA‑TTTGTTTGGAACACTAAACTTCTTTTGAGCACAACCTGTCATAATCAGGCCCAACGCacacgatg                                                                                                                                                                          >  2:129546/9‑244 (MQ=255)
                                                                  aGCGGAGTAAAAATGCCTAAAGTAATAAAACCGTGCAATCCATTTACGAATGTTTGCTGGGTGTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTGTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAg                                                                                                                                                                         <  2:110212/251‑1 (MQ=255)
                                                                      gAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                     >  1:183170/1‑251 (MQ=255)
                                                                      gAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTGTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                     <  1:300043/251‑1 (MQ=255)
                                                                       tgTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTACTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGAATCGACAGTTTTCTTCTGCCAAATTCAAGAAACGAAGAAAT‑GATGGGTGATGGTTTACTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCATGGAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                     >  2:120353/2‑250 (MQ=255)
                                                                        gTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAt                                                                                                                                                                   <  1:446150/251‑1 (MQ=255)
                                                                         tataAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAACTTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GCTGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACATCTGTTGAGCACATCCTGTAATAAGCAGGGCCCGCTCAGTAGCGAGTAGCAtt                                                                                                                                                                  >  2:398027/1‑251 (MQ=255)
                                                                         tataAATGACTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGCTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGCTTTCCTTTGGTGCTACTGCTGACTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCACGCGCAGTAGCGAGTAGCAtt                                                                                                                                                                  >  2:73562/1‑251 (MQ=255)
                                                                                   tAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATgg                                                                                                                                                        >  1:190103/1‑251 (MQ=255)
                                                                                    aaaGTAATAAACCCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTGCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGgt                                                                                                                                                       <  2:439443/251‑1 (MQ=255)
                                                                                      aGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTACTGCGCCGCCACAAATTTTGGATGCATCGACAATTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGAAGGGCCAGCGCAGTAGCGTGTAGCATTTTTTTCATGgtgt                                                                                                                                                     >  2:106712/1‑251 (MQ=255)
                                                                                       gTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGtt                                                                                                                                                    <  1:398027/251‑1 (MQ=255)
                                                                                          ataaAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTAtt                                                                                                                                                 >  2:229662/1‑251 (MQ=255)
                                                                                           taaAACCGAGCATTCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AAAAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATACTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTc                                                                                                                                                >  1:181737/1‑251 (MQ=255)
                                                                                            aaaaCCGAGCACTCCATTTACGCCTGTTTGCTGGGGTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCACATCCAGAAACGACGAAAT‑GCTGGGAGATGGTTTCCTTTGATGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTTTTGAGAACATCCTATAATAAGCAGGGCCAGCGAAGTAGCGAGTAGAATTGTTTTCATGGTGTTATTcc                                                                                                                                               >  2:447769/1‑251 (MQ=255)
                                                                                                ccGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGAt                                                                                                                                           <  1:184910/251‑1 (MQ=255)
                                                                                                 cGAGCAATCCATTTCAGAATAGTTGCTTGGTGTCTGTTTC‑AACAACATTTTCTGCTCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGTAGTAGCGAGTAGCATTTTTTTCATTGTGTTATTCCCGATg                                                                                                                                          <  1:220455/251‑1 (MQ=255)
                                                                                                      aacccttttACAAATGTTTGCGGGGGTTCTGTTTT‑AACAAAATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCtttt                                                                                                                                     <  2:23962/245‑1 (MQ=255)
                                                                                                         ccATTTACGAATGTGTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCAACGCCACCAATTTTGGCTGCTTCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAA‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCAGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATACTGTAATAAGCTGGGCCAGCGCCGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCttttttt                                                                                                                                  >  2:275235/1‑249 (MQ=255)
                                                                                                          aatttacgaattctggttGGGCTTCGGTTTT‑AATCACATTTTCTGCGCCGCTACAATGTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCTTTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                                                                                                 <  1:11834/234‑1 (MQ=255)
                                                                                                             ttACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTCGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGtt                                                                                                                              >  2:224594/1‑251 (MQ=255)
                                                                                                              tACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACACATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCCGAAACGACGACAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGCACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCCTGGTGTTATTCCCGATGCTTTTTGCAGTTc                                                                                                                             >  1:163670/1‑251 (MQ=255)
                                                                                                              tACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGtt                                                                                                                              <  1:111356/250‑1 (MQ=255)
                                                                                                               aCGATTGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTc                                                                                                                             <  1:141407/250‑1 (MQ=255)
                                                                                                               aCGAATGTTTGCTGGGTTTCTGGTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTc                                                                                                                             <  1:423709/250‑1 (MQ=255)
                                                                                                                  aaggTTTGTTGGGTTTCTGTTTT‑AATAACATTGTCTGCACCGGCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGc                                                                                                                           <  1:208762/246‑1 (MQ=255)
                                                                                                                  aaTGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGc                                                                                                                           >  1:240812/1‑249 (MQ=255)
                                                                                                                   aTGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGACACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGGTACTTTTTGAAGTTCGc                                                                                                                           >  2:1758/1‑248 (MQ=255)
                                                                                                                      tttGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                     >  1:184441/1‑251 (MQ=255)
                                                                                                                       tttctgagTTTCTATTTT‑AACATCATTTTCGCCGCCGTCAAAAATTTTGGCTGCATCGACAGTTTTTTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTTGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                     <  2:372076/243‑1 (MQ=255)
                                                                                                                       ttGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCg                                                                                                                    >  1:192530/1‑251 (MQ=255)
                                                                                                                       ttGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                     <  1:451266/250‑1 (MQ=255)
                                                                                                                        tGCTGGGTT‑CTGTTTTGAAATACATTTTCTGGGCCGCCACAAATTTTTGCGGTTTCTACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTTCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGTGCCAGCGCAGTATCGAGTAGCATGTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCg                                                                                                                    <  2:330704/250‑1 (MQ=255)
                                                                                                                         gCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGt                                                                                                                   >  1:75668/1‑250 (MQ=255)
                                                                                                                          cTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACACATTTTGGCGGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGACGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGAAAACGTCTGTTGAGCACATCCTGAAATAAGCAGGGCCAGCGCAGTAGCGCGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATctcat                                                                                                                 >  2:65189/1‑247 (MQ=255)
                                                                                                                              gTTTCTGTTTT‑AACAACATTTTCTGTGCCTCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATTGGTGATGGGTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATgtgt                                                                                                             <  2:402955/251‑1 (MQ=255)
                                                                                                                                        t‑AACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGATTTCTTCTGCCCAATTCCAGAAACGCAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAACATTTTTTTCATGGTGTTATTACCGATGCTTTTTGAAGTTTGCAGATTCGTATGTGGAGAAAAttc                                                                                                    >  2:349739/1‑249 (MQ=255)
                                                                                                                                                aTTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAAcc                                                                                            >  1:259837/1‑251 (MQ=255)
                                                                                                                                                  tttCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCt                                                                                          >  1:219781/1‑251 (MQ=255)
                                                                                                                                                    tcgtcgccgcCACAAATTTTAGCTGCATCGACAGTTTTCGTCTGCCCAAGTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTTGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                         <  2:15799/246‑1 (MQ=255)
                                                                                                                                                    tCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                         >  1:64446/1‑250 (MQ=255)
                                                                                                                                                    tCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACAGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                         <  1:6145/250‑1 (MQ=255)
                                                                                                                                                       gcgcCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTaa                                                                                        <  1:103334/248‑1 (MQ=255)
                                                                                                                                                        cgccgcCACAAATTTTGGCTGCATCGACAGTTTTCGTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGTTGATGGTTTCCTTTGGTGCTACTGCTG‑CTGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAACCCCTAAACaa                                                                                    <  1:73381/251‑1 (MQ=255)
                                                                                                                                                           cgccACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGa                                                                                 <  2:190103/251‑1 (MQ=255)
                                                                                                                                                              cgcaAAGTTTGGCTGCATGGACATTTTCCTTCTGTCCAATTCCAGAACAGAAGAAAT‑GATGGGTGATGGTTTCATTTGGGGCTACTGCTGCCGG‑TTTGGTTTGAAAAGTAAACGTCTCTTGAGCACATCCTGTAATTAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGTTGATATTCCCGATGCTTTTTGAAGATCGCAGAATCGTATGTGTAGAAAATTAAACAAACTCTAAACAATGAGtt                                                                              <  2:435644/249‑1 (MQ=255)
                                                                                                                                                                 aaatttcggCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGaa                                                                           <  1:224594/243‑1 (MQ=255)
                                                                                                                                                                 aaaTTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGa                                                                            >  1:58768/1‑250 (MQ=255)
                                                                                                                                                                 aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGa                                                                            >  2:108993/1‑250 (MQ=255)
                                                                                                                                                                  aaTTTTAGCTGCATCGACAGTTTTATTCTGCACCATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGATGTCTG‑TTTGTTTTGAACAGGAAATGTCTGTTGCGCACATCATGTAATACGCAGGGCCAGCGCAGTAGAGAGTAGCATTTTTTTCATGGTGTTCTTCACGATGCTTTTTGAAGTTCGCAGATTCGTATGTGTAGAAAATTACACAAACCCTAAACCATGAGTTGaa                                                                           >  2:157268/1‑250 (MQ=255)
                                                                                                                                                                  aaTTTTAGCTGAATCGACAGTTTTCTTCTGCCCCATTCCAGAATCGACGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACCTCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCACGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGaa                                                                           >  2:222385/1‑250 (MQ=255)
                                                                                                                                                                        ggcggCATCGACAGTTTTCTTCGGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTGTGAACAGTAAACGGCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTGGAAATTTCa                                                                     <  2:220426/246‑1 (MQ=255)
                                                                                                                                                                        aGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGGTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCa                                                                     <  1:21844/250‑1 (MQ=255)
                                                                                                                                                                         gctgctgcggcAGTTTTCTTCAGCCCAAGTCCAGAGACTAAGTAAT‑GATGGGTGATTGTTGCCTTTTGTGCTACTGCTGCCTG‑TTTGTTTTGAACCGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGCAGCATTTTTTTCATGTTGTTATTCCCGAAGCTTTTTGAAGTGCGCATAATCGTATGTGTAGAAAATTGAACAGACCCTAAACAATGCGTTGAAATTTCata                                                                   <  2:358875/241‑1 (MQ=255)
                                                                                                                                                                             cATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGCTGGTTTCCTTTGGTGCTACTGCTGCCGG‑ATTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCGGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCACGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAACCCCTAATCAATGAGTTGAATTTTCATATTGt                                                               >  2:396127/1‑251 (MQ=255)
                                                                                                                                                                                    aGTTGTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAtatt                                                        <  1:197282/251‑1 (MQ=255)
                                                                                                                                                                                     gttctcttcgGCCCAATTCCTGCAACGAAGAAAT‑GGTGGGTGATTGTTTACTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGCACCGTAAACGTCTGTTGAGCACATTCTGTAATAAGCAGTGCCAGCGCATTAGCGAGTAGCATTTGTTGCATTGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCTTATGTGTATTAAATTAAACAAACCCTAAAAAATGAGTTTAAATTTCATATTTTTAAtattt                                                       <  2:206638/247‑1 (MQ=255)
                                                                                                                                                                                               gCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCGTTGGTGCTACTGCTGCCGT‑TTTGTGTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGGTTAAATTTCATATTGTTAATATTTATTAATGt                                               <  1:334949/249‑1 (MQ=255)
                                                                                                                                                                                                   aatgccAGAAACGAAGTAAT‑GTTTGGTGACGGTTTTCTCTGGGGCTACTGCTGGCTG‑TTTGTTGTGAACAGTTAATGTCGTTTGAGCATTTCCTGTAATAAGCATGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATTGTGTTATTACCTATGCTTTTTGAATTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACTATGAGTTTAAATTTCATATTGTTAATATTGATTAATGTATGcta                                         <  2:93911/247‑3 (MQ=255)
                                                                                                                                                                                                    aTTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTGTGTTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTGCATGGTGTTATTCCCGATGCTTTTTGAAGGTCGCAGAATCGTATGTGTAGAAAATTAAACAGACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAg                                        <  1:24583/251‑1 (MQ=255)
                                                                                                                                                                                                               gaagaaAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATg                                            <  1:18953/236‑1 (MQ=255)
                                                                                                                                                                                                               gaagaaAT‑GATGGGTGATGGATTCCTTTGGTCCTACTGCTGCCGG‑TTTGTTTTGAACAGGAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCAAGCGCAGTAGAGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCTGATTCGTATGTGTAGAAAATTAAACAAACCCTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATg                                            >  2:18953/1‑236 (MQ=255)
                                                                                                                                                                                                                 ataaaT‑GATGTGTGATTGTTTCCTTTGATGATACTGCTGCCGG‑TTTGTTTTGAACAGTAAGCTTCTTTTGAGCACATCCTGTAATATGCAGGGCAAGCGCATTAGCGAGTAGCATTTTTTTCATTGTGATATTCACGATGCTTTTTGAAGTTCGCAGATGCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATTTATGTCAGGTGCTATGAATCg                           <  2:307831/249‑1 (MQ=255)
                                                                                                                                                                                                                 agaaAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCg                           <  1:4792/251‑1 (MQ=255)
                                                                                                                                                                                                                 agaaAA‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTGGAACAGTAAATGTCTGTTTAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCg                           <  2:339694/251‑1 (MQ=255)
                                                                                                                                                                                                                  gaaAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCTTATGTGTAGAAAATTAAACAAACACTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGt                          <  1:438454/251‑1 (MQ=255)
                                                                                                                                                                                                                           gggTGATGGTTTCCTTTTGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCTATGAATCGTCATTGTAt                  <  2:317661/251‑1 (MQ=255)
                                                                                                                                                                                                                           gggTGATGGTTTCCTTTGGTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTAt                  >  1:164347/1‑251 (MQ=255)
                                                                                                                                                                                                                               gattgTTTCCTTTGGTGCTTCTGCTGCCGG‑TTTGTTTTTAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTcc               <  1:14814/246‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGAt          <  1:233165/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ttcttttGGTTCTACTGCTGCCTG‑TTTGGTTTGAACAGTAAACGTCTGTTGAGCACATCCTTTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAATTTCGCAGAATCGTATGTGTAGAAAATTTAACAATCCCTAAACAATGAGTTGAAATTTCATTTTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGAtt         <  1:37049/246‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTTTGGTGCTTCTGCCGCCGG‑TTTGGTTTGAACAGTAAACGTCTGGTGAGCACTTCCTTTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATTGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAAc      <  2:250253/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggTGCTACTGCTGCCGG‑TTTGTTTTGAACAGTAATCGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAAGCGTCATTGTATTCCCGGATTAACTATg  <  1:336681/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
CAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT‑AACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAT‑GATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTG‑TTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATG  >  NC_000913/578562‑579045

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