Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 67,596 T→C K248E (AAA→GAA)  araA ← L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091367,5960TC97.8% 281.8 / ‑2.3 89K248E (AAA→GAA) araAL‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (41/46);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

GAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTT  >  NC_000913/67349‑67830
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
gaggTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGcc                                                                                                                                                                                                                                         >  2:164334/1‑251 (MQ=255)
        gccCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGAttt                                                                                                                                                                                                                                 <  1:360152/251‑1 (MQ=255)
             gCAGACCGGTTGAAAACACCTTCCAGATCGGAAGCAGGGCGGCGTTTTTCCAGTCGCCTGCGCCCGCAAAGCCGTAACCCTGCGGCCTCAGACGCTGTAGGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGT‑AAGCTACCTTGTTCCA‑GGAATCGCTTCATCCCCAGCTC‑AATACGCGCCGCTGCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTtgtgtt                                                                                                                                                                                                                            <  2:263108/251‑2 (MQ=255)
                        gACATCACCTTAATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGTTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCAt                                                                                                                                                                                                                 <  2:238008/251‑1 (MQ=255)
                         ccATCTCCTTCATGAAGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTGCGCCGTGGATTTGTGTTGCAGGCGTCATg                                                                                                                                                                                                                <  2:399304/250‑1 (MQ=255)
                             cACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCGTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATTCGCGCCGCTTCCAGCACGTTCTGTCTTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGgtg                                                                                                                                                                                                             <  1:587519/250‑1 (MQ=255)
                             cACCGTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCCGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGACCGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGgtgt                                                                                                                                                                                                            <  1:699692/251‑1 (MQ=255)
                                     tgatgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCtt                                                                                                                                                                                                    >  1:460072/1‑251 (MQ=255)
                                       atgagAAGCATGGCTGCTGTTTTCCCTTCGCCTGCTGCCGCATAGCAGTAACCCTGCTGCATCACACGCTGTACTGCCTGACCTGGAAGCTGTTTCATACCG‑TGTAAATCTTCAAAGGGTGTGGTGACCGCGTGG‑AAGCGTCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCTTTGTTCGCCGTTGATTTGTGTGTCAGGcttcttggggtatcacctgta                                                                                                                                                                                                  <  2:479236/247‑21 (MQ=255)
                                       atgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACATGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGAGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGAGTAGCAGCTTTc                                                                                                                                                                                                  >  1:79539/1‑251 (MQ=255)
                                            aaGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACt                                                                                                                                                                                             <  1:146945/251‑1 (MQ=255)
                                             aGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc                                                                                                                                                                                            <  1:236986/251‑1 (MQ=255)
                                             aGCAGGGCGGCGGGTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGTCCTGGCAGCTGTTTAAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGGGG‑AAGCCACCTGGTTCCA‑GGAAACGCTTCATACCCAGCTC‑AATACGCGCCGCTTCCAGCCAGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTTCAGTCGTCATGGTGTAGCCTCTTTCGTACTc                                                                                                                                                                                            >  2:161559/1‑251 (MQ=255)
                                                 gggcggcggtgttccagacggtgtCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCCCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCGTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCg                                                                                                                                                                                        <  2:385415/228‑1 (MQ=255)
                                                  ggcggcGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCg                                                                                                                                                                                        >  1:245108/1‑250 (MQ=255)
                                                          ttttccagcggcttccccccGCAAAGCGTTAACCCGCCTCCACCAGCCGCTGTCCGGCAGGACCGGGAGGTTGTTTCAAACCG‑TGCAAACCTCAAAAGGTGGGGTTGAACGCGTGG‑AAGCCCCCTT‑TTCCAAGGAACCGTTTCTTCCCCAGCTA‑AATCCGCGCGGCTTC‑AGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGcgcg                                                                                                                                                                               <  2:139648/233‑1 (MQ=255)
                                                             tCCAGTCGCCGTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGTCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGGTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGtt                                                                                                                                                                             <  2:258273/250‑1 (MQ=255)
                                                              ccAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTaa                                                                                                                                                                           <  1:37222/251‑1 (MQ=255)
                                                              ccAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTa                                                                                                                                                                            >  1:415795/1‑250 (MQ=255)
                                                                 gTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACAt                                                                                                                                                                        >  1:620651/1‑251 (MQ=255)
                                                                   cGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAACGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCATCTTTCGTACTCATCGACCAGCGCGTTAACAtgg                                                                                                                                                                      >  2:502182/1‑249 (MQ=255)
                                                                   cGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCg                                                                                                                                                                      >  1:413884/1‑251 (MQ=255)
                                                                   cGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTAAGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCg                                                                                                                                                                      >  1:588747/1‑251 (MQ=255)
                                                                    ccctgcGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑ATGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGGTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGc                                                                                                                                                                     <  2:464917/246‑1 (MQ=255)
                                                                     ccTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCa                                                                                                                                                                                                                         >  1:434410/1‑198 (MQ=255)
                                                                     ccTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGAGTTGTGTTGCa                                                                                                                                                                                                                         <  2:434410/198‑1 (MQ=255)
                                                                         cGCCCGCAAGGTCGTAACCCTGCGGCATCAGCCGCTGTACGGCCAGTCCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAATGTGGGGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTGCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTc                                                                                                                                                                 <  2:710408/250‑1 (MQ=255)
                                                                           cccGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCt                                                                                                                                                              <  2:607654/251‑1 (MQ=255)
                                                                                     cGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTc                                                                                                                                                    >  1:152808/1‑251 (MQ=255)
                                                                                      gTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGAGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGCTTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCCCGATAACATCGCCGTCGCTGATgcattta                                                                                                                                                    >  2:20964/1‑244 (MQ=255)
                                                                                         aCCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTcacca                                                                                                                                                <  1:62477/251‑1 (MQ=255)
                                                                                          cccTGCTGCATCAGACGCTGTACGGACAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAAAGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTcaccac                                                                                                                                               >  1:228523/1‑251 (MQ=255)
                                                                                             tgctgcATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTg                                                                                                                                            >  1:482550/1‑251 (MQ=255)
                                                                                              tcgtcatcaggtgcgtgaCGGCCAGACCGTGAATCTGGTTCAGACTT‑TGCAAATCTGCAAAGGTGGTGGTGAACGCGTGG‑AAGCCGTCTTTTTTCT‑GGAAACTCTTCATGCCCAGCTC‑AATACGCGCCGTTTCCAGCACGTTCTGTCGGTTGTCGCCGTGGATTTGGGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGGTAACAGCGCCGTCGCTTATGGAGCTcacctccg                                                                                                                                             <  2:704344/232‑5 (MQ=255)
                                                                                                 gcATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCAc                                                                                                                                         <  1:560468/250‑1 (MQ=255)
                                                                                                   aTCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCAcc                                                                                                                                        >  1:161711/1‑249 (MQ=255)
                                                                                                       gACGCTGTACGGCCAGACCTGGAAGCTGGTTCAGACCG‑TGCATATCTACAAAGGTGGTGGTGAAAGCGTGG‑AAGCCGCCTGGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTGCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCg                                                                                                                                  <  2:269347/251‑1 (MQ=255)
                                                                                                        aCGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGAGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCCGATCGc                                                                                                                                 >  2:708316/1‑251 (MQ=255)
                                                                                                         cGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGcc                                                                                                                                >  2:556401/1‑251 (MQ=255)
                                                                                                           cTGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCa                                                                                                                               >  1:8114/1‑250 (MQ=255)
                                                                                                            tGTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAAc                                                                                                                             >  1:248847/1‑251 (MQ=255)
                                                                                                             gTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCaa                                                                                                                              <  1:18783/249‑1 (MQ=255)
                                                                                                             gTACGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAAcc                                                                                                                            >  1:310712/1‑251 (MQ=255)
                                                                                                                cGGCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGcc                                                                                                                         >  1:656795/1‑251 (MQ=255)
                                                                                                                  gCCAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGGCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCa                                                                                                                       <  2:142712/251‑1 (MQ=255)
                                                                                                                   ccAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGccc                                                                                                                        >  1:539311/1‑249 (MQ=255)
                                                                                                                    cAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCCCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGGTCTGTCGTTTTGCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAgg                                                                                                                     <  2:660688/251‑1 (MQ=255)
                                                                                                                    cAGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAgg                                                                                                                     >  1:367349/1‑251 (MQ=255)
                                                                                                                     aGACCTGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGt                                                                                                                    <  1:604908/251‑1 (MQ=255)
                                                                                                                           gtaaGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGCTTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGATTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGAc                                                                                                              >  1:394786/3‑251 (MQ=255)
                                                                                                                            gAAGCTGTTTCAGACCGATGCAA‑TCTTCAAAGGTGGTGGTGAACGCGTTG‑AAGCTCCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGATTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTTATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGaag                                                                                                             <  2:152808/251‑3 (MQ=255)
                                                                                                                               gCTGTTTCAGCCCG‑GGCAAATCCTCAACGGGGGGTGTGAACGC‑TGGTAAGCCACCTATTTTCA‑GGAAACGCTTCGTCCCCAG‑TCTAATACGCGCCGCTGCCAGTCCGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACAGCTCCGTCGCTGATGTAGTTCACCCCCTGCACCAGATCGCCAACCGACTAGGTATTGACGgag                                                                                                          <  2:421592/251‑1 (MQ=255)
                                                                                                                                   tttCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCa                                                                                                                                                                                                                  >  1:110888/1‑143 (MQ=255)
                                                                                                                                   tttCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCa                                                                                                                                                                                                                  <  2:110888/143‑1 (MQ=255)
                                                                                                                                     tCAGACCG‑TGCAAATCTTCAAAGGGGGGGGTGAACGCGTAG‑AAGCCTGCTTGTTCCA‑GGAACCGCCTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTTTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGAGGGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAGGTATTGACGGCGAAAcc                                                                                                     <  2:336907/250‑1 (MQ=255)
                                                                                                                                      cAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCg                                                                                                                                                                <  1:24164/190‑1 (MQ=255)
                                                                                                                                      cAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGa                                                                                                   <  1:385982/251‑1 (MQ=255)
                                                                                                                                       aGACCG‑TGCAAATCTTCAAAGGAGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATCGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAGGGATTGACGGAGAAACCGa                                                                                                   >  2:433701/1‑250 (MQ=255)
                                                                                                                                         aCCG‑TGCAAATCTTCAAATGTGGTGGTGAACGCGTGG‑AAGCCCCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACt                                                                                                <  2:430798/251‑1 (MQ=255)
                                                                                                                                          ccG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACtt                                                                                               <  2:595258/251‑1 (MQ=255)
                                                                                                                                            g‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTCGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGAAGTTCACCACCTGCACCAGATCGCCAACCGCCCCGGTATTGACGGCGAAACCGAACTTg                                                                                              >  2:606337/1‑250 (MQ=255)
                                                                                                                                               gCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACAGCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATc                                                                                           <  2:645998/251‑1 (MQ=255)
                                                                                                                                                 aaaTCTTCAACGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCC‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCTCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAAAAACCGAACTTGctctt                                                                                         >  2:634976/1‑250 (MQ=255)
                                                                                                                                                  aaTCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCtgt                                                                                        <  1:281478/251‑1 (MQ=255)
                                                                                                                                                  aaTCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCtgt                                                                                        <  1:560212/251‑1 (MQ=255)
                                                                                                                                                       tCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGATACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGc                                                                                      <  2:165335/248‑1 (MQ=255)
                                                                                                                                                            ggtggtggtGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACtt                                                                              >  1:512341/1‑251 (MQ=255)
                                                                                                                                                               ggtggggAACGCGTGT‑AAGCCCCCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGAGCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTAt                                                                           <  2:656795/251‑1 (MQ=255)
                                                                                                                                                                gtggtgAGCGCGTGG‑AAGCCACCTTGGTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATGCTCGCCGCTTCCAGCTCGTTTTTACGTTTGTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCTACCAGCGCGTTAACATCGCCGTCGCGGATGGAGTTCACCACCTGCACTAGATCGCCATCCGCCCAGGTATTGACGGAGAAACCGAACTTGTTCTGTGCGGCAACTTTATc                                                                          <  2:14933/251‑1 (MQ=255)
                                                                                                                                                                 tggtgAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGCGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATAGACCAGCGCGTTAACATCGCCGTCGCTGATGGGGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGGGAAACCGAACTTGAGCTGTGCGGCAACTTTATCg                                                                         >  2:277001/1‑251 (MQ=255)
                                                                                                                                                                  ggtgAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATca                                                                         <  1:165345/250‑2 (MQ=255)
                                                                                                                                                                   gtgAACTCGTGG‑AAGCCACCTTGTCCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATcacc                                                                       <  1:502182/251‑4 (MQ=255)
                                                                                                                                                                      aaCGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATc                                                                    <  1:498733/251‑1 (MQ=255)
                                                                                                                                                                       aCGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCg                                                                   <  1:164334/251‑1 (MQ=255)
                                                                                                                                                                            tGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGAc                                                              >  1:55350/1‑251 (MQ=255)
                                                                                                                                                                             gg‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGAcc                                                             <  1:520251/251‑1 (MQ=255)
                                                                                                                                                                                     cccTTTTTCCA‑GGAAACGCTTCACCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCACCGGTGACCGCCACtt                                                      <  1:600691/250‑1 (MQ=255)
                                                                                                                                                                                         tGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCTGCAACTTTATCAACATCGGGGACCGCCACTTCACg                                                  >  2:197128/1‑251 (MQ=255)
                                                                                                                                                                                         tGTTACA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCCGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCACGTTAACATCGCCGTCGCTGCTGGAGTTCACCACCTGCACCAGATCGCCACCCGCCCAGGTATTGGCGGAGAAACCGAACTTTATCTGTGCGGCACCTTTATAACAAACGATGACCGCCACTTCACg                                                  >  2:106805/1‑251 (MQ=255)
                                                                                                                                                                                          gTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCTGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGCTCGCCAACCGCCCAGGTATTGACGGAGAAACCGACCTTGATCTGTGCGGCAACTTTATCCACATCGATGACCGCCACttaccgc                                                 >  2:261346/1‑246 (MQ=255)
                                                                                                                                                                                              cA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGAGTTGACGGAGAAACCGCACTTGATCTGTGCGGCAACTTTCTCACCATCGGTGACCGCCGCTTCACGCATGt                                             >  2:265365/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaCGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATc                                         >  1:416179/1‑251 (MQ=255)
                                                                                                                                                                                                    aaccgcttcaccccaagtcaaatccgcgCCGTTTCCAGCATTTTCTGTCTTTTTTCGCCGTGGATTTTTGTTGCAGCTGTCAGGGTGTAGCAGTTTTCGTCCTCACCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTCCCCCACCTGCCCCAGACCCCCACCCGCCAAGGTATTGACGGAGAAACCTAACTTGATCTGGGCGGCAACTTTATCCCCATCGGTGACCGCCACTTCACGCATGTTATc                                         <  1:421967/227‑1 (MQ=255)
                                                                                                                                                                                                           cATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGCACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCaaaac                                 >  2:270389/1‑249 (MQ=255)
                                                                                                                                                                                                                ccAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGc                              >  1:109411/1‑249 (MQ=255)
                                                                                                                                                                                                                        aaTACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGAAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGCCCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTc                      >  2:701615/1‑250 (MQ=255)
                                                                                                                                                                                                                            cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCaaaac                                 >  1:84604/1‑233 (MQ=255)
                                                                                                                                                                                                                                  gcTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCGAACCGCCCAGGt                                                                                                                    <  1:217278/146‑1 (MQ=255)
                                                                                                                                                                                                                                            cGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTAACCACCTGCACCAGATAGCCAACCGCCCAGGTATTGCCGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACAGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGttt  >  1:468504/1‑250 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
GAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCG‑TGCAAATCTTCAAAGGTGGTGGTGAACGCGTGG‑AAGCCACCTTGTTCCA‑GGAAACGCTTCATCCCCAGCTC‑AATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTT  >  NC_000913/67349‑67830

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