Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,563,2370TC70.3% 33.5 / 46.9 68V173A (GTT→GCT) LY180_07445transposase IS609
Reads supporting (aligned to +/- strand):  ref base T (7/12);  new base C (19/26);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGT  >  CP006584/1563086‑1563531
                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       
gAAAACAGTCGGTTGTTTCTGGCGAAATTGGGCGGGATGCGGTTCCGGAATAGCCGTCAGGTCACGGGGGTTGTAAAAAATGTCCCTGTACGCCAGGTCTGCGGTAAGTGGTACATCATTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTGCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                       <  2:9200/298‑1 (MQ=11)
gAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGTGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTAAATCAGTATTCAGACAGAAAGTGAAGTATCCACTACGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGAGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTATTTGAGCCTGTACACAGTTTTCAGAAAAACAAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                       >  1:369776/1‑298 (MQ=255)
gAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                       <  1:543988/298‑1 (MQ=14)
 aaaaCCGCCGTTTTTTTCGGCCGAAACTGGGCTGGAGGCGCTACCGGAACAGCCTGAAGGTCACGGGTGTGGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCATACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGt                                                                                                                                                      <  2:397505/298‑1 (MQ=11)
 aaaaCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCTCTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTTAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGt                                                                                                                                                      <  1:395867/298‑1 (MQ=14)
 aaaaCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGt                                                                                                                                                      <  1:93296/298‑1 (MQ=31)
 aaaaCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGt                                                                                                                                                      <  1:2525/298‑1 (MQ=14)
   aaCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCa                                                                                                                                                    >  1:602500/1‑298 (MQ=31)
     cgggcGTTTTTTTCAGCTGGAGCGGGGCTGGATGCGTTAACGGAAGAGCCGCCAGGTCACGGGTGTGGTGAAAAATGTGAGTGTCAGGCAGGCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGGGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTGTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCaaa                                                                                                                                                  <  2:898232/294‑1 (MQ=11)
     cAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCaaa                                                                                                                                                  >  1:515278/1‑298 (MQ=31)
       gCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCCATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACCCTGCCCGATGGCACAGTCTCTGAGCCTGTAAACAGTCTCAAGAAAAACCCAAAGACCCCGGCCAAACTCCACACACAGTTCAGCCGCAAGGTCAAAtt                                                                                                                                                >  2:235619/1‑298 (MQ=11)
         cGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCa                                                                                                                                              <  1:972315/298‑1 (MQ=31)
         cGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCa                                                                                                                                              <  1:637436/298‑1 (MQ=31)
         cGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTAAGCGAGTCCTGCGGTAAGTGGTACATCAGTATTAAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGATTTCAGAAAAACCAGAAGACGCTGACGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAAttaa                                                                                                                                              >  1:613727/1‑296 (MQ=9)
          gTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAg                                                                                                                                             <  1:977666/298‑1 (MQ=31)
                  cTGCCGAAACTGGGTTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACt                                                                                                                                     <  1:176053/298‑1 (MQ=14)
                         aagtGGGCTGGATGCGCTACCGGAACGGCCGTTAGGTCAGGGGGGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGgcagaa                                                                                                                              <  2:677563/295‑1 (MQ=9)
                         aaCTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAGCAACTGgcagaa                                                                                                                              >  1:553186/1‑298 (MQ=14)
                              ggCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTCTCCACTCCGGTTCACCCTTCAGCATCAATGGTCCGGCTGGCTCCTTGCGTGGCTACACTCGCCAGGCTGTCAGATGGCACATCCTTTGAACCCGTCAACATTTTTCATAAAAACCAGAAGACCACGGCGAGACTTCAGCGAACGTTAAGCCGCCAGGTCAAATTAAGCAACACCTGgcagaaaccaa                                                                                                                         >  2:244295/1‑293 (MQ=11)
                                cTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGCCAAATTCAGCAACAACTGGCAGAAGCAGaaa                                                                                                                       <  1:585032/298‑1 (MQ=14)
                                  ggATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAACCGCAAGGGCACATTCAGCAACCACCGGCAGAAGAAGaaaac                                                                                                                     >  2:67292/1‑296 (MQ=9)
                                   gATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGACGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGCAAACAGTTGTCAGAACAACCAGAAGACACTGGCGAGACGTCAGACAAATTTAAGCCGCAAGGGCAAATTCAGCAACAACGGGCCGAAGACGAAACga                                                                                                                    >  2:635333/1‑297 (MQ=11)
                                       ggggccaggAAGAGCCGTGCAGGCACGGGGGTTGTGAAAAATGGCACTGTCAGCCGGGCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCGTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCaaaa                                                                                                                <  2:108534/291‑1 (MQ=11)
                                        gcTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAAt                                                                                                               >  1:301320/1‑298 (MQ=31)
                                          tACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATAc                                                                                                             >  1:100388/1‑298 (MQ=14)
                                           aCCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACa                                                                                                           <  1:66143/298‑1 (MQ=31)
                                                 aTAGCCGTCAGGTCACTGTGGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTAAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACt                                                                                                     <  1:237493/298‑1 (MQ=14)
                                                      ccgtacggtccaggggGTTTTGAAAAATGTCACTTGTAGCCCTTACTGCGGTAAGTGGTACATCAGTATTCAGCAAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTTGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGACTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCAt                                                                                                 <  2:1000687/286‑1 (MQ=11)
                                                      cGTCAGGTCACGGGTGGTGTTAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGCTACATCAGTATTCAGCCAGAAAGAGACATATCCCCTCCGGCTCACCCTACAGCAGCAATGGTAGGGCTGGATGCTGGCGTGGCTACACGCGCCACGCTGTCAGCTGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAATACCAGAAGACACTGGCGAGACTTACGCGCCAGTTCAGCAGCAAGGTCAACTCCAGCAACCAAGGGCAGAAGCAGCAACGCCAAAGACacagccaccatc                                                                                                >  2:1020274/1‑287 (MQ=11)
                                                      cGTCAGGGTACGGGGGGTGTGAAAAATTTCACTGGCAGACAGTCCTGCGGTAAGGGGTACATCAGTATTCAGACAGAAAGTGAAGGATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCAtt                                                                                                <  2:981562/298‑1 (MQ=11)
                                                       gTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACGGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTc                                                                                               <  1:779673/298‑1 (MQ=14)
                                                        tCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTcc                                                                                              <  1:999216/298‑1 (MQ=31)
                                                        tCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTcc                                                                                              >  1:246720/1‑298 (MQ=14)
                                                          cAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCt                                                                                             <  1:998121/298‑1 (MQ=14)
                                                          cAGGTAGCGGGTGTTGTGAAAAATGTCAATGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCt                                                                                             <  2:341376/298‑1 (MQ=14)
                                                           aGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                       <  1:355508/240‑1 (MQ=14)
                                                           aGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                       >  2:355508/1‑240 (MQ=14)
                                                                 cGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCATGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATTAATGGTCGGGCAGGATGCTGGCGTGGCTAAACTCGTTACGCTGTCAGATGGCACAGTCTGTGAGCCTGTAAACAGTTATCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGCTAAGCCGCAAGGGCAAATTCAGCAACAACTGGCAGAAGCAGCAACGCAAAATACAGCGACTGCATTCCGGTATCGc                                                                                      >  1:277484/1‑298 (MQ=11)
                                                                   ggTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACt                                                                                                     >  2:1009066/1‑281 (MQ=14)
                                                                   ggTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACt                                                                                                     <  1:1009066/281‑1 (MQ=14)
                                                                     tgttgtGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAat                                                                                  >  1:879888/1‑298 (MQ=17)
                                                                            tcttccgatcttGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCa                                                                           <  1:639559/287‑1 (MQ=17)
                                                                              aaaTGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAgg                                                                         <  1:745592/298‑1 (MQ=17)
                                                                              aaaTGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAgg                                                                         <  1:838014/298‑1 (MQ=17)
                                                                                     aCTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTAcc                                                                  <  2:440571/298‑1 (MQ=17)
                                                                                      cgggcggccatgccTTCGGTAAGGCGTACATCAGGATGCAGAAAGAAAGTGAAGTATCCACTCCGGGTCAACCTTCAGCATCAATGGGCGGGCTTGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCt                                                                 <  2:192528/286‑1 (MQ=11)
                                                                                       tGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCtt                                                                <  1:859429/298‑1 (MQ=17)
                                                                                         gcAGCCAGTCCGGCGGGAAGTGGTACATCAGTAGTCAGACAGAAAGTGAAGTATCCCCTCCGGCGCACCCTTCAGCATCAATGGTCGGGCTGGAGGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTca                                                              <  2:930100/297‑1 (MQ=12)
                                                                                            gCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCAGCAAGGAcgaacctcccaa                                                           >  2:204234/1‑287 (MQ=12)
                                                                                             ccAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACaaa                                                          <  2:208812/298‑1 (MQ=17)
                                                                                              cAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAg                                                         >  1:861845/1‑298 (MQ=17)
                                                                                                gTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTACACTCCCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGCAAACAGTTTTCAGAAAAACAGGAAGACGCTGGCGAGACTCAACCGCCAGTTACGCCCCAAGGTAAATTACAGCACCAAATGGCAGAAGCAGAAACACAAATACAAGCACAGGCCTTCCGTTAACACAAATAACCGAAGGGACTACATACACAAAgca                                                       >  2:384531/1‑296 (MQ=11)
                                                                                                      cGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTAACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACCCGCCACGATGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTAAGAAAAACAAGAAGACGCTGGCGAGAATTCAGCGACAGTTAAGACGCAAAGTCAAATTCAGCAACAACTGGCAGACGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAACATCCGCAGGGACTAAATTCACAAAGTCACAACg                                                 >  1:308695/1‑298 (MQ=11)
                                                                                                        gTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGgcagaa                                                                                                                              <  2:252034/220‑1 (MQ=14)
                                                                                                        gTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGgcagaa                                                                                                                              >  1:252034/1‑220 (MQ=14)
                                                                                                          aaGTGGTACACCAGTATGCAGACAGAAAGTGAAGTGTCCACTCCGGGTCCCCCTTCAGCATCAATGGTCGGGCGTGATGCTGGCGCGGCTAAACTCGCCACGCTGTAAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATGCAGCAACAACTGGCAGAAGCAGAAACGCAAAATGCAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCg                                             <  2:629403/298‑1 (MQ=11)
                                                                                                                         aTTCGGACAGAAAGTGAGGGATCCCCCCCGGGTCACCCTTCAGCATCAATTGTCGGGCTGGATGCTGGCGTGGCTAAACGCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACg                              <  2:40685/298‑1 (MQ=255)
                                                                                                                         aTTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACg                              >  1:29373/1‑298 (MQ=17)
                                                                                                                             agacaggACGGGAAGTATCCACTCCGGCTCACCCTTCAGCTTCAATGGTCGGGCGGGAGGCTGGCGGGGCTAAACTCGCCACGCGGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGGTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAAt                          <  2:648321/298‑1 (MQ=11)
                                                                                                                              gacagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATg                         >  1:469656/1‑298 (MQ=17)
                                                                                                                                  gaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGCAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGACACGCAAAATACAGCGACTGCATTCCTGTATCGCACATAACAGCAGGGACAACCTCCACAAAGACACACAGCACGCCCGCAAACACCACGCAATGACAg                     >  2:149339/1‑298 (MQ=11)
                                                                                                                                   aaaggcatgttgccACTCGGGGGCACCCTTAGGAATCAATGGTCGGGCTGGAGGCTGGCGTGGCTAAACTCGCCAGGCGGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGt                    <  2:390642/286‑1 (MQ=11)
                                                                                                                                     aGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCTGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCa                  <  1:7582/298‑1 (MQ=17)
                                                                                                                                       tGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCAtt                <  1:302020/298‑1 (MQ=17)
                                                                                                                                             ttCCACTCCGGCTCACCCCTCGGCGTAAATGGTCGGGCTGGATGCTGGGGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAt          <  2:861845/297‑1 (MQ=12)
                                                                                                                                              tcccccccGGCTCACCCTCTCCCATCAATGGTCGGGGTGGAGGCTGGCGTGGCTAAACCCGCCAAGCTGTCAGATGGCACAGTCTTTGGGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAGCAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAtt         <  2:553186/292‑1 (MQ=11)
                                                                                                                                                  cTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGaa     <  1:593812/298‑1 (MQ=17)
                                                                                                                                                     cGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGt  >  1:66658/1‑298 (MQ=17)
                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       
GAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGT  >  CP006584/1563086‑1563531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: