Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,820 C→T A36A (GCG→GCA)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,8200CT95.3% 266.1 / 1.5 87A36A (GCG→GCA) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (36/45);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

CCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATC  >  NC_000913/72574‑73063
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ccccTTCTGGCGGGCCATAGCGTGCATTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTGAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGAGAGCAAGTTCAGCAGGGTACTTTTACCGg‑cgc                                                                                                                                                                                                                                                 <  2:244066/251‑1 (MQ=255)
 cccATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑cgcc                                                                                                                                                                                                                                                >  2:328623/1‑251 (MQ=255)
    agcTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCt                                                                                                                                                                                                                                             <  1:637950/249‑1 (MQ=255)
       tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTgg                                                                                                                                                                                                                                          >  1:647474/1‑251 (MQ=255)
       tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAAAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTCCCTG‑CGCCGCTTgg                                                                                                                                                                                                                                          >  2:664217/1‑251 (MQ=255)
        ggcgggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTg                                                                                                                                                                                                                                           <  1:87633/249‑1 (MQ=255)
               gATAGCGTGCAGTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGa                                                                                                                                                                                                                                  <  1:103855/251‑1 (MQ=255)
                     gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGcc                                                                                                                                                                                                                             >  1:164242/1‑250 (MQ=255)
                         aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTg                                                                                                                                                                                                                        <  1:161604/251‑1 (MQ=255)
                         aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTg                                                                                                                                                                                                                        <  1:401387/251‑1 (MQ=255)
                               cccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTcgcc                                                                                                                                                                                                                  <  2:647474/251‑1 (MQ=255)
                                ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCCAGATCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCCCCTGCTcgccg                                                                                                                                                                                                                 >  2:73713/1‑251 (MQ=255)
                                ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTcgccg                                                                                                                                                                                                                 >  1:26406/1‑251 (MQ=255)
                                      gtgccGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGCGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTct                                                                                                                                                                                                           >  1:447914/4‑250 (MQ=255)
                                            cgTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAAGGCGTAACGGCGGCGTCGTGGTGTGATAAACGCCGTCGAGAGTCAGCGGACCGCTGGCTGTCGTCAGAAAACCGACGATCAAATTCAGAAGGGTACTTTTACCGG‑CGCCTCTTAGCCCGAGGATCGCCACCTGCACGCCGCGTTCTACCGtt                                                                                                                                                                                                     >  2:467399/1‑251 (MQ=255)
                                             gTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCATGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCCAGATCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGTGCGCCACCTGCTCGCCGCGTTCTACCggt                                                                                                                                                                                                     >  2:532623/1‑248 (MQ=255)
                                               tCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTa                                                                                                                                                                                                    >  1:391171/1‑249 (MQ=255)
                                                 cgttacaaccccGGATTTAGCCCCAGCCCGAGGGTCTGTGCGACCGTCAGGTGGCTGACCAGGTTGTCCTCCTGAAACATCATCGACACCTGACGGCGGTACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCTCGAGGATCGCCTCCTGCTCGCCGCGTTCCACCTTTAAGCt                                                                                                                                                                                                <  2:55711/242‑1 (MQ=255)
                                                        aTCCCGGATTTAGCCCCAGCCCGATTTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGc                                                                                                                                                                                         <  2:618288/251‑1 (MQ=255)
                                                          cccGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCAt                                                                                                                                                                                       >  1:120129/1‑251 (MQ=255)
                                                           ccGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATc                                                                                                                                                                                      <  1:344071/251‑1 (MQ=255)
                                                            cGGATTTAGCCCCATCCCGATGTTCCGCGCGCCCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGTCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTGGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGGTAATCTAAAACGCATCg                                                                                                                                                                                     <  2:311618/251‑1 (MQ=255)
                                                             ggATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCgg                                                                                                                                                                                    <  1:607253/251‑1 (MQ=255)
                                                                 ttAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCaa                                                                                                                                                                                 <  1:418436/250‑1 (MQ=255)
                                                                   aGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCGGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtg                                                                                                                                                                              <  2:140477/251‑1 (MQ=255)
                                                                      cccATCCCGATGTTCTGCGCGACCGTCAGGTGAATGAACAGGTTTTTTTCCTGCAACAGCATCGACACCGGTCGGCGTGACGGCGGCGTCGTGTTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGGCGATCAATTTCAGCAGGGTACTTTTTCCTG‑CGCCGTTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttg                                                                                                                                                                           <  2:566792/251‑4 (MQ=255)
                                                                          gCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTaa                                                                                                                                                                       <  1:328623/251‑1 (MQ=255)
                                                                              gATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGcc                                                                                                                                                                   >  1:322304/1‑251 (MQ=255)
                                                                                   tCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGt                                                                                                                                                               <  2:164242/250‑1 (MQ=255)
                                                                                        gCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAgggggtatt                                                                                                                                                         >  2:177651/1‑244 (MQ=255)
                                                                                         cGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtca                                                                                                                                                        <  1:274140/251‑1 (MQ=255)
                                                                                         cGACCGTCAGGTGACTGAACAGGTTGTTTTCCGGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCTGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtca                                                                                                                                                        <  1:664217/251‑1 (MQ=255)
                                                                                              gTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcagtcag                                                                                                                                                   <  2:337658/251‑1 (MQ=255)
                                                                                                   gtgacTGAACGGGTTGTTTTCCTGAATCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGTACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGGTCAGCAGGGTACTTTTACCTG‑CGCCGCTGGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTa                                                                                                                                              <  2:588556/247‑1 (MQ=255)
                                                                                                    tgacTGAACAGGTTGTTTTCCTGAAACAGCATCTACACCTGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                                          <  2:320787/219‑1 (MQ=255)
                                                                                                    tgacTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                                          >  1:320787/4‑222 (MQ=255)
                                                                                                      acTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAAc                                                                                                                                            >  1:562290/2‑250 (MQ=255)
                                                                                                       ctgtaCAGGTTGTTTTCCTGAAACAGCATCGACACAGGACGGCGTGAAGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAACAGGCTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTACCCGGGATCTCCCGTGAAACACAAACTCATCGCCAAATGTTGGTAAAGCCGGGTGATATACGCAAGTTTTAACAt                                                                                                                                          >  2:511308/5‑251 (MQ=255)
                                                                                                         gAACAGGTTTTTTTCCTGAAACAGCAGCGACACCGGACGGCGTGACGGCGGCGGCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAtc                                                                                                                                         <  2:391171/250‑2 (MQ=255)
                                                                                                              ggTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGAGTGATCACCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCCTCGGCAAATGTTGGTAAAGCCCAGTGGTATCAGTCAGTTTTAACATTTCGccc                                                                                                                                   >  2:425865/1‑251 (MQ=255)
                                                                                                              ggTTGTTGTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGccc                                                                                                                                   <  1:291200/251‑1 (MQ=255)
                                                                                                               gttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                                                                  <  1:44122/251‑1 (MQ=255)
                                                                                                                    tttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTa                                                                                                                              <  1:207164/248‑1 (MQ=255)
                                                                                                                     ttcCTGAAACATCATCGACACCGGACGGAGTGACGGCGGGGTCGTGGTGTGAACAACGCCGTCGATAACCAGCGATCCGCTGGCTGGCGCCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGttg                                                                                                                             <  2:524112/249‑3 (MQ=255)
                                                                                                                     ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAg                                                                                                                             <  1:83965/249‑1 (MQ=255)
                                                                                                                             ccAGCATCGACACAGGACGGCGTGCCGGCTGCGTCGTGGTGTGGTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGAACAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTTAGCTAAAACGCATCGGTCAATGTTGGTAAAGCCAGGAGCTATCAGTCAGTTTTAACCTTTCGCCCCGGTAATTTTTccac                                                                                                                    >  2:301675/2‑248 (MQ=255)
                                                                                                                             aCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAAt                                                                                                                    >  1:699840/1‑251 (MQ=255)
                                                                                                                              cAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGAAATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATc                                                                                                                   >  1:234547/1‑251 (MQ=255)
                                                                                                                               aGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAATAAGCGAACCGCTGGCTGGCGTCAGAAAACAAGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGATCGACGCGTTCTACCGTTAAGCTAAAACGCATCGGCAACTGTTGGTAAAGACAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTTTCAAta                                                                                                                   >  1:703972/1‑249 (MQ=255)
                                                                                                                                  atcgaacccgtccGGCGTGACGGCGCGATCGTGGTGTGATCAAGGCCGTCGATAGTCAGCGAGCCGCGGGCTGGCGTCATAAACCCGTCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCGTGCCGCGTGCTACCGTTAAGCGAAAACGCAGCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTGTTCAATCACgg                                                                                                               <  2:491846/239‑1 (MQ=255)
                                                                                                                                        cccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAca                                                                                                         <  1:530579/250‑1 (MQ=255)
                                                                                                                                         ccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcag                                                                                                        <  2:454654/251‑1 (MQ=255)
                                                                                                                                         ccGGACGGCGTGACGGCGGCGTCGTGGTGTGAGCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTGCAATCACGGTAAAcag                                                                                                        <  1:219768/251‑1 (MQ=255)
                                                                                                                                             aCGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                                                       >  1:17259/1‑248 (MQ=255)
                                                                                                                                                               cgtggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGCTCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcagc                                                                                     >  2:366371/5‑248 (MQ=255)
                                                                                                                                                                   gtgtgATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                               >  2:220108/1‑250 (MQ=255)
                                                                                                                                                                   gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACCGAGTAGCAGCAGAAtt                                                                               >  2:140136/1‑250 (MQ=255)
                                                                                                                                                                   gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCCGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                               >  1:653924/1‑250 (MQ=255)
                                                                                                                                                                    tgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATAGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTa                                                                              >  1:469040/1‑250 (MQ=255)
                                                                                                                                                                     gtgATCAACGACGTCGATAGTCAGCGAACCGATGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCGCGGTAGTGTTTCAATCACGGTAAACAGTAGAAAACAGAGTAGCAGCAGAATTAAc                                                                            <  1:237695/251‑1 (MQ=255)
                                                                                                                                                                          caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                       <  1:513750/249‑1 (MQ=255)
                                                                                                                                                                          caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTGACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                       <  2:83218/249‑1 (MQ=255)
                                                                                                                                                                           aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                                      >  1:702875/2‑251 (MQ=255)
                                                                                                                                                                            acgccggcGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGa                                                                     <  1:563089/244‑1 (MQ=255)
                                                                                                                                                                                     aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGt                                                            >  1:30549/1‑251 (MQ=255)
                                                                                                                                                                                       aGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                                          >  1:209714/1‑139 (MQ=255)
                                                                                                                                                                                       aGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgttgg                                                                                                                                                                          <  2:209714/139‑1 (MQ=255)
                                                                                                                                                                                        gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCt                                                         <  1:486554/251‑1 (MQ=255)
                                                                                                                                                                                            gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt                                                     <  2:653924/251‑1 (MQ=255)
                                                                                                                                                                                                  cGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGAt                                                >  1:67130/1‑250 (MQ=255)
                                                                                                                                                                                                     gggctggcGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGTCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCTGTCCTGGCTGCGATAGGa                                            <  1:156338/250‑1 (MQ=255)
                                                                                                                                                                                                      ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                                           <  1:73713/251‑1 (MQ=255)
                                                                                                                                                                                                              tCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAtt                                     <  1:309548/249‑1 (MQ=255)
                                                                                                                                                                                                                 gAAAACCGGCTATCAAGTTCAGCAGGGTACTTG‑ACCCGGCGCCGCTTGGCCCGAGGATCGCTACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                <  2:26406/251‑1 (MQ=255)
                                                                                                                                                                                                                     aCCGGCGAGCAAGTTACGCAGGGTACTTTTCCCTG‑CGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTag                            <  2:93357/251‑1 (MQ=255)
                                                                                                                                                                                                                           gATCAAGTTCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCAACTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACAGCGCCGTACTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATa                       >  1:219478/1‑250 (MQ=255)
                                                                                                                                                                                                                                agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                 >  1:365887/3‑219 (MQ=255)
                                                                                                                                                                                                                                agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                 <  2:365887/217‑1 (MQ=255)
                                                                                                                                                                                                                                 gttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACACGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGCGATAAAACGGc                >  1:232299/2‑251 (MQ=255)
                                                                                                                                                                                                                                  ttCAGCGGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                >  2:600197/1‑250 (MQ=255)
                                                                                                                                                                                                                                   tCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                >  1:628535/1‑249 (MQ=255)
                                                                                                                                                                                                                                      gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt           <  1:425865/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt           <  1:274904/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gggTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg        <  2:30549/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                              <  1:344441/223‑1 (MQ=255)
                                                                                                                                                                                                                                               ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                              >  2:344441/1‑223 (MQ=255)
                                                                                                                                                                                                                                               ttttACCTG‑CGCCGCTTGGCACGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGAGCGGAAatc  >  1:564100/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATC  >  NC_000913/72574‑73063

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