Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,8310TC73.4% 129.9 / 54.3 95T21T (ACA→ACGborDprophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (32/37);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTG  >  NC_000913/578587‑579081
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          
gCAACTCATGCAATTATTGTGAACAATACACACGCGCTTCCAGCGGCGTATAAATACCTAAAGTAATAAAACCGCGCAATCCATTTACGACTGTTTTCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTTCCCAATTCCAGAAACGAAGAAATGATGGGTGTTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAGCAGTAAACGTCttt                                                                                                                                                                                                                                                      >  2:428243/1‑249 (MQ=255)
  aaCTCATGCCATTATTTTGCGAAATACACACGCGCCTCCAGCGGAGTCTAACCGCCTAAAGTAATAACCCCGAGCAATCCATTTACGAGTGTTTGCGGGGTTTCTGTTTTACCAACATTTTCTGAGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGCAGAAATGATGGGTGTTGGTTTCCTTTGGTGCTTCTGCTGTCTGTTTGTTTTGAACAGTAATTGTCTGttt                                                                                                                                                                                                                                                    >  2:451266/1‑250 (MQ=255)
    cTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTg                                                                                                                                                                                                                                                    >  1:418547/1‑249 (MQ=255)
           aaTTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCACTTCCAGAAACGAAGAAATGATGGTTGATGGTTTCCTTTGGTGCTACTGCTGCCGTTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATcc                                                                                                                                                                                                                                           >  2:4792/1‑251 (MQ=255)
                 tgtGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGACTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGAGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACCTACTGtaat                                                                                                                                                                                                                                     >  2:312373/3‑251 (MQ=255)
                    gAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAg                                                                                                                                                                                                                                  <  1:117971/251‑1 (MQ=255)
                     aGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGc                                                                                                                                                                                                                                 <  1:184126/251‑1 (MQ=255)
                            cccactcgtttCCAGCGGGTTATAAATGCCTAATGTAATAAAACCGACGACTCCTTTTACGAATGTTTGCTGGGTTTCTTTTTTAACATCATTTTCTGCGCCGCTACAAATTTTAGCTGCATCGTCAGTTTTCTTCTGCCCAATTCCAGATACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCCGGGCCa                                                                                                                                                                                                                          <  2:390001/242‑1 (MQ=255)
                            cacccgctcttccatggtaGTAAAAATTCCTAAAGTAATAAAACCGTGCAATCCATTTACGAATGTTTGATGGGTTTCTGTTTTAAAACCATTTTCAGCGCCGCCACAAATGTTATCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAGACGAAGAAATGATGGGTGATGGTTGACTTTGGTGCTGCTGCTGTCTGTTTGTTATGAACAGTAAAGGTCTGTTGAGCACTTCCTGTAATAAGCAGGGCCa                                                                                                                                                                                                                          <  2:323021/233‑1 (MQ=255)
                             acacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGAGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGcc                                                                                                                                                                                                                           >  1:229617/1‑249 (MQ=255)
                              cacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCCGGGCCAgc                                                                                                                                                                                                                        >  2:110172/1‑251 (MQ=255)
                                      tCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                                                                >  1:250253/1‑251 (MQ=255)
                                       ccAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                                                                <  1:318392/250‑1 (MQ=255)
                                        taGCGGAGAATTAATGACGAAAGTAATAAAACCGAGCAATCCATTTACGAACGTTTGCTGGGTTTCTGTTTTAACAAAATTTTCTGCGCCGACACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGCAAAGAAGAAATGATGGGAGGTGGTTTCTTTTGTTGCGACTGCTGCCGGTTGGTTTTGAGTAGTAAACGGCTGTTGTGCACATCCTGTAATAAGCTGGGCCAGAGCAGTAgat                                                                                                                                                                                                               >  2:14814/2‑248 (MQ=255)
                                        aagctgaatataAATGCCTAAAGTAATAAAACCGAGCAATCAATTTACGAACGATTGATGGGATTCTATATCAACAACAATTCCTGCGCCGCCACAAATTTTGGCAGCCTCGACAGTCATCTTCTGCCCCATTCCAGCAACGATGAAATTGACAATGATGATATCATTTGTTGCTACTGCTGACGATTTGTTTGGAACACTAAACTTCTTTTGAGCACAACCTGTCATAATCAGGCCCAACGCacacgatg                                                                                                                                                                                                              >  2:129546/9‑244 (MQ=255)
                                         aGCGGAGTAAAAATGCCTAAAGTAATAAAACCGTGCAATCCATTTACGAATGTTTGCTGGGTGTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTGTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAg                                                                                                                                                                                                             <  2:110212/251‑1 (MQ=255)
                                             gAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                                                         >  1:183170/1‑251 (MQ=255)
                                             gAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTGTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                                                         <  1:300043/251‑1 (MQ=255)
                                              tgTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTACTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGAATCGACAGTTTTCTTCTGCCAAATTCAAGAAACGAAGAAATGATGGGTGATGGTTTACTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCATGGAATAAGCAGGGCCAGCGCAGTAGCGAGTAGc                                                                                                                                                                                                         >  2:120353/2‑250 (MQ=255)
                                               gTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAt                                                                                                                                                                                                       <  1:446150/251‑1 (MQ=255)
                                                tataAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAACTTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGCTGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACATCTGTTGAGCACATCCTGTAATAAGCAGGGCCCGCTCAGTAGCGAGTAGCAtt                                                                                                                                                                                                      >  2:398027/1‑251 (MQ=255)
                                                tataAATGACTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGCTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGCTTTCCTTTGGTGCTACTGCTGACTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCACGCGCAGTAGCGAGTAGCAtt                                                                                                                                                                                                      >  2:73562/1‑251 (MQ=255)
                                                          tAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATgg                                                                                                                                                                                            >  1:190103/1‑251 (MQ=255)
                                                           aaaGTAATAAACCCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTGCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGgt                                                                                                                                                                                           <  2:439443/251‑1 (MQ=255)
                                                             aGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTACTGCGCCGCCACAAATTTTGGATGCATCGACAATTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGAAGGGCCAGCGCAGTAGCGTGTAGCATTTTTTTCATGgtgt                                                                                                                                                                                         >  2:106712/1‑251 (MQ=255)
                                                              gTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGtt                                                                                                                                                                                        <  1:398027/251‑1 (MQ=255)
                                                                 ataaAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTAtt                                                                                                                                                                                     >  2:229662/1‑251 (MQ=255)
                                                                  taaAACCGAGCATTCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAAAAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATACTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTc                                                                                                                                                                                    >  1:181737/1‑251 (MQ=255)
                                                                   aaaaCCGAGCACTCCATTTACGCCTGTTTGCTGGGGTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCACATCCAGAAACGACGAAATGCTGGGAGATGGTTTCCTTTGATGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTTTTGAGAACATCCTATAATAAGCAGGGCCAGCGAAGTAGCGAGTAGAATTGTTTTCATGGTGTTATTcc                                                                                                                                                                                   >  2:447769/1‑251 (MQ=255)
                                                                       ccGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGAt                                                                                                                                                                               <  1:184910/251‑1 (MQ=255)
                                                                        cGAGCAATCCATTTCAGAATAGTTGCTTGGTGTCTGTTTCAACAACATTTTCTGCTCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGTAGTAGCGAGTAGCATTTTTTTCATTGTGTTATTCCCGATg                                                                                                                                                                              <  1:220455/251‑1 (MQ=255)
                                                                             aacccttttACAAATGTTTGCGGGGGTTCTGTTTTAACAAAATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCtttt                                                                                                                                                                         <  2:23962/245‑1 (MQ=255)
                                                                                ccATTTACGAATGTGTGCTGGGTTTCTGTTTTAACAACATTTTCTGCAACGCCACCAATTTTGGCTGCTTCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAAAGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCAGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATACTGTAATAAGCTGGGCCAGCGCCGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCttttttt                                                                                                                                                                      >  2:275235/1‑249 (MQ=255)
                                                                                 aatttacgaattctggttGGGCTTCGGTTTTAATCACATTTTCTGCGCCGCTACAATGTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCTTTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                                                                                                                                     <  1:11834/234‑1 (MQ=255)
                                                                                    ttACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTCGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGtt                                                                                                                                                                  >  2:224594/1‑251 (MQ=255)
                                                                                     tACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACACATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCCGAAACGACGACATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGCACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCCTGGTGTTATTCCCGATGCTTTTTGCAGTTc                                                                                                                                                                 >  1:163670/1‑251 (MQ=255)
                                                                                     tACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGtt                                                                                                                                                                  <  1:111356/250‑1 (MQ=255)
                                                                                      aCGATTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTc                                                                                                                                                                 <  1:141407/250‑1 (MQ=255)
                                                                                      aCGAATGTTTGCTGGGTTTCTGGTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTc                                                                                                                                                                 <  1:423709/250‑1 (MQ=255)
                                                                                         aaggTTTGTTGGGTTTCTGTTTTAATAACATTGTCTGCACCGGCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGc                                                                                                                                                               <  1:208762/246‑1 (MQ=255)
                                                                                         aaTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGc                                                                                                                                                               >  1:240812/1‑249 (MQ=255)
                                                                                          aTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGACACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGGTACTTTTTGAAGTTCGc                                                                                                                                                               >  2:1758/1‑248 (MQ=255)
                                                                                             tttGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                                                         >  1:184441/1‑251 (MQ=255)
                                                                                              tttctgagTTTCTATTTTAACATCATTTTCGCCGCCGTCAAAAATTTTGGCTGCATCGACAGTTTTTTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTTGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                                                         <  2:372076/243‑1 (MQ=255)
                                                                                              ttGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCg                                                                                                                                                        >  1:192530/1‑251 (MQ=255)
                                                                                              ttGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                                                         <  1:451266/250‑1 (MQ=255)
                                                                                               tGCTGGGTTCTGTTTTGAAATACATTTTCTGGGCCGCCACAAATTTTTGCGGTTTCTACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTTCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGTGCCAGCGCAGTATCGAGTAGCATGTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCg                                                                                                                                                        <  2:330704/250‑1 (MQ=255)
                                                                                                gCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGt                                                                                                                                                       >  1:75668/1‑250 (MQ=255)
                                                                                                 cTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACACATTTTGGCGGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGACGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGAAAACGTCTGTTGAGCACATCCTGAAATAAGCAGGGCCAGCGCAGTAGCGCGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATctcat                                                                                                                                                     >  2:65189/1‑247 (MQ=255)
                                                                                                     gTTTCTGTTTTAACAACATTTTCTGTGCCTCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATTGGTGATGGGTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATgtgt                                                                                                                                                 <  2:402955/251‑1 (MQ=255)
                                                                                                               tAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGATTTCTTCTGCCCAATTCCAGAAACGCAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAACATTTTTTTCATGGTGTTATTACCGATGCTTTTTGAAGTTTGCAGATTCGTATGTGGAGAAAAttc                                                                                                                                        >  2:349739/1‑249 (MQ=255)
                                                                                                                       aTTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAAcc                                                                                                                                >  1:259837/1‑251 (MQ=255)
                                                                                                                         tttCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCt                                                                                                                              >  1:219781/1‑251 (MQ=255)
                                                                                                                           tcgtcgccgcCACAAATTTTAGCTGCATCGACAGTTTTCGTCTGCCCAAGTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTTGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                                                             <  2:15799/246‑1 (MQ=255)
                                                                                                                           tCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                                                             >  1:64446/1‑250 (MQ=255)
                                                                                                                           tCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACAGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTa                                                                                                                             <  1:6145/250‑1 (MQ=255)
                                                                                                                              gcgcCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTaa                                                                                                                            <  1:103334/248‑1 (MQ=255)
                                                                                                                               cgccgcCACAAATTTTGGCTGCATCGACAGTTTTCGTCTGCCCAATTCCAGAAACGAAGAAATGATGGTTGATGGTTTCCTTTGGTGCTACTGCTGCTGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAACCCCTAAACaa                                                                                                                        <  1:73381/251‑1 (MQ=255)
                                                                                                                                  cgccACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGa                                                                                                                     <  2:190103/251‑1 (MQ=255)
                                                                                                                                     cgcaAAGTTTGGCTGCATGGACATTTTCCTTCTGTCCAATTCCAGAACAGAAGAAATGATGGGTGATGGTTTCATTTGGGGCTACTGCTGCCGGTTTGGTTTGAAAAGTAAACGTCTCTTGAGCACATCCTGTAATTAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGTTGATATTCCCGATGCTTTTTGAAGATCGCAGAATCGTATGTGTAGAAAATTAAACAAACTCTAAACAATGAGtt                                                                                                                  <  2:435644/249‑1 (MQ=255)
                                                                                                                                        aaatttcggCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGaa                                                                                                               <  1:224594/243‑1 (MQ=255)
                                                                                                                                        aaaTTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGa                                                                                                                >  1:58768/1‑250 (MQ=255)
                                                                                                                                        aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGa                                                                                                                >  2:108993/1‑250 (MQ=255)
                                                                                                                                         aaTTTTAGCTGCATCGACAGTTTTATTCTGCACCATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGATGTCTGTTTGTTTTGAACAGGAAATGTCTGTTGCGCACATCATGTAATACGCAGGGCCAGCGCAGTAGAGAGTAGCATTTTTTTCATGGTGTTCTTCACGATGCTTTTTGAAGTTCGCAGATTCGTATGTGTAGAAAATTACACAAACCCTAAACCATGAGTTGaa                                                                                                               >  2:157268/1‑250 (MQ=255)
                                                                                                                                         aaTTTTAGCTGAATCGACAGTTTTCTTCTGCCCCATTCCAGAATCGACGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACCTCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCACGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGaa                                                                                                               >  2:222385/1‑250 (MQ=255)
                                                                                                                                               ggcggCATCGACAGTTTTCTTCGGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTGTGAACAGTAAACGGCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTGGAAATTTCa                                                                                                         <  2:220426/246‑1 (MQ=255)
                                                                                                                                               aGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGGTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCa                                                                                                         <  1:21844/250‑1 (MQ=255)
                                                                                                                                                gctgctgcggcAGTTTTCTTCAGCCCAAGTCCAGAGACTAAGTAATGATGGGTGATTGTTGCCTTTTGTGCTACTGCTGCCTGTTTGTTTTGAACCGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGCAGCATTTTTTTCATGTTGTTATTCCCGAAGCTTTTTGAAGTGCGCATAATCGTATGTGTAGAAAATTGAACAGACCCTAAACAATGCGTTGAAATTTCata                                                                                                       <  2:358875/241‑1 (MQ=255)
                                                                                                                                                    cATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGCTGGTTTCCTTTGGTGCTACTGCTGCCGGATTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCGGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCACGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAACCCCTAATCAATGAGTTGAATTTTCATATTGt                                                                                                   >  2:396127/1‑251 (MQ=255)
                                                                                                                                                           aGTTGTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAtatt                                                                                            <  1:197282/251‑1 (MQ=255)
                                                                                                                                                            gttctcttcgGCCCAATTCCTGCAACGAAGAAATGGTGGGTGATTGTTTACTTTGGTGCTACTGCTGCCGGTTTGTTTTGCACCGTAAACGTCTGTTGAGCACATTCTGTAATAAGCAGTGCCAGCGCATTAGCGAGTAGCATTTGTTGCATTGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCTTATGTGTATTAAATTAAACAAACCCTAAAAAATGAGTTTAAATTTCATATTTTTAAtattt                                                                                           <  2:206638/247‑1 (MQ=255)
                                                                                                                                                                      gCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCGTTGGTGCTACTGCTGCCGTTTTGTGTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGGTTAAATTTCATATTGTTAATATTTATTAATGt                                                                                   <  1:334949/249‑1 (MQ=255)
                                                                                                                                                                          aatgccAGAAACGAAGTAATGTTTGGTGACGGTTTTCTCTGGGGCTACTGCTGGCTGTTTGTTGTGAACAGTTAATGTCGTTTGAGCATTTCCTGTAATAAGCATGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATTGTGTTATTACCTATGCTTTTTGAATTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACTATGAGTTTAAATTTCATATTGTTAATATTGATTAATGTATGcta                                                                             <  2:93911/247‑3 (MQ=255)
                                                                                                                                                                           aTTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTGTGTTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTGCATGGTGTTATTCCCGATGCTTTTTGAAGGTCGCAGAATCGTATGTGTAGAAAATTAAACAGACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAg                                                                            <  1:24583/251‑1 (MQ=255)
                                                                                                                                                                                      gaagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATg                                                                                <  1:18953/236‑1 (MQ=255)
                                                                                                                                                                                      gaagaaATGATGGGTGATGGATTCCTTTGGTCCTACTGCTGCCGGTTTGTTTTGAACAGGAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCAAGCGCAGTAGAGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCTGATTCGTATGTGTAGAAAATTAAACAAACCCTAAACACTGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATg                                                                                >  2:18953/1‑236 (MQ=255)
                                                                                                                                                                                        ataaaTGATGTGTGATTGTTTCCTTTGATGATACTGCTGCCGGTTTGTTTTGAACAGTAAGCTTCTTTTGAGCACATCCTGTAATATGCAGGGCAAGCGCATTAGCGAGTAGCATTTTTTTCATTGTGATATTCACGATGCTTTTTGAAGTTCGCAGATGCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATTTATGTCAGGTGCTATGAATCg                                                               <  2:307831/249‑1 (MQ=255)
                                                                                                                                                                                        agaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCg                                                               <  1:4792/251‑1 (MQ=255)
                                                                                                                                                                                        agaaAAGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTGGAACAGTAAATGTCTGTTTAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTATCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCg                                                               <  2:339694/251‑1 (MQ=255)
                                                                                                                                                                                         gaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCTTATGTGTAGAAAATTAAACAAACACTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGt                                                              <  1:438454/251‑1 (MQ=255)
                                                                                                                                                                                                  gggTGATGGTTTCCTTTTGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCTATGAATCGTCATTGTAt                                                      <  2:317661/251‑1 (MQ=255)
                                                                                                                                                                                                  gggTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTAt                                                      >  1:164347/1‑251 (MQ=255)
                                                                                                                                                                                                      gattgTTTCCTTTGGTGCTTCTGCTGCCGGTTTGTTTTTAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTcc                                                   <  1:14814/246‑1 (MQ=255)
                                                                                                                                                                                                          gTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGAt                                              <  1:233165/251‑1 (MQ=255)
                                                                                                                                                                                                            ttcttttGGTTCTACTGCTGCCTGTTTGGTTTGAACAGTAAACGTCTGTTGAGCACATCCTTTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAATTTCGCAGAATCGTATGTGTAGAAAATTTAACAATCCCTAAACAATGAGTTGAAATTTCATTTTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGAtt                                             <  1:37049/246‑1 (MQ=255)
                                                                                                                                                                                                              ccTTTGGTGCTTCTGCCGCCGGTTTGGTTTGAACAGTAAACGTCTGGTGAGCACTTCCTTTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATTGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAAc                                          <  2:250253/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggTGCTACTGCTGCCGGTTTGTTTTGAACAGTAATCGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAAGCGTCATTGTATTCCCGGATTAACTATg                                      <  1:336681/250‑1 (MQ=255)
                                                                                                                                                                                                                            gctgccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGACGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGccc                            >  1:383257/8‑251 (MQ=255)
                                                                                                                                                                                                                            gctgccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGcc                             <  1:229662/243‑1 (MQ=255)
                                                                                                                                                                                                                                  ggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACgg                      >  1:195279/2‑251 (MQ=255)
                                                                                                                                                                                                                                     tttgtttTGAACAGAAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGAGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACgggg                    >  1:83240/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttttGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACt                >  1:250936/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tttGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACtt               >  1:445303/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACtt               <  1:73562/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTCATGGTGTTATTCACGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAAAAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAgtg  >  1:186969/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          
GCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTG  >  NC_000913/578587‑579081

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