Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,305 A→C 11.5% H130P (CAC→CCC)  yfjX → CP4‑57 prophage, putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,3050AC11.5% 194.7 / 2.1 80H130P (CAC→CCC) yfjXCP4‑57 prophage, putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (40/29);  new base C (2/9);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGG  >  NC_000913/2776067‑2776543
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
tttgCCTGGGTGGACCGTTTTTGGGATGGGTACTGGGGGGGTG‑TCTGGTCCTTTTACCCGCTCAGCAATGGCGGAGGTTTTATGGCTCCGGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCa                                                                                                                                                                                                                                                <  2:539488/236‑1 (MQ=255)
ttCGCCTGGATGGACCGCTTCTGGGATGAGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCa                                                                                                                                                                                                                                                <  1:661691/239‑1 (MQ=255)
    ccTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCCCAAAGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGGAGCCGCAGGA‑AGCGC‑TGTCTGCCTGATTGAATATAGCCAGCCAGCCTGCCTCACCGAAAGCGATGCCATGAC‑GGAGCACTAt                                                                                                                                                                                                                                            >  1:1080696/1‑239 (MQ=255)
             aCCGCTTCTGTGATGAGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATAGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATTGTGCGGGAATG‑AGCGCGGAAGCCGCAGGG‑ATCGC‑GTTCTGCCTGGTTGAAAATAGCCATCACGCCGGCCCCACCCAATGCGATTCCCTTAC‑GGGAAACCATTAccccgt                                                                                                                                                                                                                                    >  2:419749/1‑234 (MQ=255)
              agttttgtataagggccggcgggggggTATTTTGGCCCTTTTCCCCGCCCAAGAAGGGCGGGGCGTGTGTGGCCCCTGGTGTTGA‑GGGTGACGATAAAGGGCGTTGGCTCAGCGGCATGAATGGCACTGGTGCGGAAAGG‑AGCGGGGAAGCCGCAGGA‑ATCGT‑TGTCTGCCTGATGGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTg                                                                                                                                                                                                                                   <  2:572364/217‑1 (MQ=255)
                       ttgtgtgttcggggggggTT‑TTTGGTCCTTTTACCCGCTCAGAAATGGGGGGGGTTTTAGGGCTCCGGATGCTGG‑CGGGGACGAAAAATGGCTTTTGCTCAACGGCATGAATGGCAATGGTGCGAAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACg                                                                                                                                                                                                                         <  2:664269/226‑1 (MQ=255)
                          ggatgccggggggtgTC‑TCTGGTCCTTTTACCCGCCAAGGAAGGGCGGGGGGTTTAGGGTCCCTGATGCTGA‑CGGTGGCGATAAGTGGCTTTGGCTCAACGGCATGGAGGGAAATGGGTCGGAAATG‑GGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCACCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGct                                                                                                                                                                                                                       <  2:519643/228‑1 (MQ=255)
                           gaGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTAAGc                                                                                                                                                                                                                        >  1:852022/1‑236 (MQ=255)
                            aGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCCTGAATGGGAATGGTGCGGAAATG‑GGCCCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGCTTGAATATAGCCATACAGCCCGCCTCACCGAATTCGAGGCCATGAC‑GGGGACCTATTACGCCCTCCGGGGTTACGccctc                                                                                                                                                                                                                    >  2:230645/1‑238 (MQ=255)
                                 tGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCCATGGTGCGGAAATG‑AGGGCGGAAGCCGCAGGG‑ATCGC‑TGTCTGCCTGATTGAATAAAGCCATCACGCCTCCCCCACCGAAGGCGATGCCCTGCG‑CGAGCACACTTACccccccggggttaaccccctgcaca                                                                                                                                                                                                                >  2:891437/1‑217 (MQ=255)
                                   cGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGGAATG‑AGCGCGGAAGCCGCAGGG‑ATCGC‑TGTCTGCCTGATTGAAAATAGCCATACAGCCTGCCCCACCGCATGCGATGCCATTAC‑GGGGCACCATTACCCCCCGCGGGAATACGCCCCGCAacac                                                                                                                                                                                                              >  2:901875/1‑238 (MQ=255)
                                       gtgtttttGGTTCTTTTACCCGCTCAGAAAGGGGGGGGGGTTTAGGGCTCCTGGTGCGGG‑GGGTGGAGAAAAAGGGCATTTGCCAAACGGAATGAATGGCAATGGGGGGGAAATG‑GGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCt                                                                                                                                                                                                           <  2:352663/231‑1 (MQ=255)
                                        gTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGG‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACCATTACCCCCTGCGGGATTACGCTCTGCACAACCCCGaa                                                                                                                                                                                                        >  2:980815/1‑239 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGGACTATTTCCCCCTGCGGGATTACg                                                                                                                                                                                                                         >  1:572364/1‑219 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGGTTACGCTCTGCAGCACCCTGAATc                                                                                                                                                                                                      >  1:253788/1‑238 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGG‑AACGC‑TGTCTGCCTGATTGAATATAGCCCTCACGCCTGCCTCACCGAATGCCATGCCATGAC‑GGGGCCCCTTTACCCCCTCCGGGATTCCGCTCTGCACCACCCCGAAAcc                                                                                                                                                                                                     >  2:27541/1‑239 (MQ=255)
                                              ttgTCCTTTATCCCGCTCAGCATGGGGGGGGGTTTTAGGGCTCATGGTG‑TGAGGGGGGCCGAAAAAGGGCTTTGCCTAAACGGAAGGAATGGAAAGGGGGGGGAAAGG‑ACCGGGGAAGCCGGAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAg                                                                                                                                                                                                   <  2:991986/237‑1 (MQ=255)
                                                gTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAgcg                                                                                                                                                                                                 <  1:25143/239‑1 (MQ=255)
                                                 tCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGAGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCCACGGGACTACGCCTTGCAGCACCCTGGATCCAgcg                                                                                                                                                                                                 >  1:202743/1‑238 (MQ=255)
                                                     tttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGCTTACGCTCTGCAGCACCCTGAATCCAACGCCCtt                                                                                                                                                                                            >  1:38061/1‑239 (MQ=255)
                                                      ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGAAATG‑AGCGCGGAAGCCGCAGGA‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACCCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGAACTATTAcc                                                                                                                                                                                                                                        >  1:913817/1‑194 (MQ=255)
                                                              tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGGTGCCATGAC‑GGAGAACTATTTCCGCCTGCGGGCTTAACCCCTCCACCCCCCCGAAACCGACGCCCTTATggcccta                                                                                                                                                                                   >  2:857273/1‑233 (MQ=255)
                                                               cAGCAATGGCGGGGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGAATGAATGGCAATGGTGCGGAAAGG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatta                                                                                                                                                                                  <  1:118004/239‑1 (MQ=255)
                                                               aaggaaTGGCGGGGCTTTTAGGGCTCCTGATGCTGA‑CGGGGACGATAAAGGGCATTTGCTCAACGGCATGAATGGAAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatta                                                                                                                                                                                  <  1:729691/235‑1 (MQ=255)
                                                                     tGGGGGGGCGTTTATGGCCCCTGATGCGGA‑CGGTAACGATAAATGGCATTTGCTCAAGGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAtt                                                                                                                                                                            <  1:857273/239‑1 (MQ=255)
                                                                          gAGCGTTTATGGCTCCTGATGCTGA‑CGGGGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGACTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGAACTATAACCCCCTGCGGGATTACGCCCTGCAGCCCCCCGAATCCAGCCCCATTATGGCCattatttgataaggg                                                                                                                                                                       >  2:758399/1‑230 (MQ=255)
                                                                               tttATGGCTCCGTGTGCTGC‑CGGTGACGATAAAGGGCATTTGCTCAACGGCATGAATGGAAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTTTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:672129/238‑1 (MQ=255)
                                                                               tttATGGCTCCGGATGCTGA‑CGGGTACGATAAAGGGCTTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:891437/238‑1 (MQ=255)
                                                                                   tGGCTCTTGAGGGTTA‑GGGTGACGATAAATGGCTTTTGCTAAACGGCATGAATGGCAATGGTGCGGAAATG‑GGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTTATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatga                                                                                                                                                              <  1:336831/239‑1 (MQ=255)
                                                                                        gcTGATGCTGA‑GGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAAGG‑GGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGaa                                                                                                                                                         <  1:351001/238‑1 (MQ=255)
                                                                                         cTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTCCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTCTGCGCATTATTTATTAAGGGAACCCAAGATGGaaa                                                                                                                                                        >  2:3627/1‑238 (MQ=255)
                                                                                         cTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTCCGGGATTACGCCTTGCAGCCCCCTGAATCCAGCCCCATTTTGCGCCTTTTTGATTAAGGaaaacactaatgaaaa                                                                                                                                                        >  2:68548/1‑224 (MQ=255)
                                                                                              gCTGA‑CGGTGACGAAAAATGGCATTGGCCCAACGGCATGAAGGGCAATGGTGCGGAAATG‑GGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:601224/223‑1 (MQ=255)
                                                                                               cTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCGGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAA‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGAAGCCCCCTGAATCCAGCGCCCTTATCCGCCTTATTGATTAAGGAAACTCATGATGGAACAGtaa                                                                                                                                                   >  2:24632/1‑236 (MQ=255)
                                                                                                  gagggTGACGATAAATGGCATTTGCTAAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTc                                                                                                                                                <  1:27541/236‑1 (MQ=255)
                                                                                                     ggTGACGATAATTGGTATGTGCCCAACGGGTGAAATGGCAAGGGTGGGAAAATG‑AGCCCGGACGCCGGCGGT‑AACGC‑TGTCGGCCTGATTGAATATAGCCAACCCGCCTGCCTCCCCGAATGCGATGCCATGGC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCAt                                                                                                                                              <  2:38061/238‑1 (MQ=255)
                                                                                                         aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCAATATGCGCATTTTTGATTAAGGAAACTCATGATGGAAccggtctcttataaac                                                                                                                                         >  2:601224/1‑224 (MQ=255)
                                                                                                             aaaaTGGCATTTGCCCAAGGGCATGAATGGCAATGGGGCGGAAATG‑AGCGGGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      <  1:758399/237‑1 (MQ=255)
                                                                                                                aTGGCATTGGCCCAACGGCTGGACGGGCAAGGGTGCTGAAATG‑AGCGCGGAAGCCGAAGGT‑ATGCC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCAGGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACAcc                                                                                                                                  <  1:1002310/239‑1 (MQ=255)
                                                                                                                     ttttGCTCACGGGTATGAAT‑GCAATGGGGCCGAAAAGTGGCGGGGAAGCC‑CGGGTGATCGC‑TGTCTGCCGGTTTGATTATGGCCATCCCGCCAGCCTCACCGAATGCGGTGCCCGGGC‑GGAGCACTTTTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTAATGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGCCACCGGTa                                                                                                                              <  2:182166/237‑1 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAt                                                                                                                                                                             <  1:407891/191‑1 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGt                                                                                                                               >  1:292103/1‑237 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTAATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTTTTGAt                                                                                                                                                                             >  2:407891/1‑191 (MQ=255)
                                                                                                                      tttGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACGGAACCTCATCCCACAGACACCGGAAc                                                                                                                             >  1:1031109/1‑238 (MQ=255)
                                                                                                                              aCGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgc                                                                                                                                                                                        >  2:415947/1‑171 (MQ=255)
                                                                                                                              aCGGCATGAATGGAAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgc                                                                                                                                                                                        <  1:415947/171‑1 (MQ=255)
                                                                                                                              aCGGAATGAATGGCAAGGGGGCGGAAATG‑AGCGCGGAAGCCGAAGGT‑ATCGC‑TGTCTGCCTGTTTGAATAGAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCTCTGa                                                                                                                    <  1:121851/239‑1 (MQ=255)
                                                                                                                               cGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      <  2:1041596/220‑1 (MQ=255)
                                                                                                                               cGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      >  1:1041596/1‑220 (MQ=255)
                                                                                                                                      aagggcaggggTGCGAAAATG‑GCCCGGGAGGCCGCAGGT‑‑TCGCGTGTCGGCCTGATTGAATAAAGCCATCACCCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGcaca                                                                                                             <  1:368903/229‑1 (MQ=255)
                                                                                                                                             aTGGTGCGGAAAGG‑AGCGCGGAAGCCCCAGGT‑ATCGG‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACgg                                                                                                     <  1:901875/239‑1 (MQ=255)
                                                                                                                                                gTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTAACTCATCCCACAGCAACCGTTACTTCCACTGACCCCACAAAGCACGGGaa                                                                                                  >  2:708043/1‑239 (MQ=255)
                                                                                                                                                    ggAAATG‑AGCGCGGAAGCCGCAGGT‑ACCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGTCACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:418554/239‑1 (MQ=255)
                                                                                                                                                     gAAATG‑GGCGCGGAAGCCGCAGGT‑ATCGG‑TGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCCCCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:68548/238‑1 (MQ=255)
                                                                                                                                                     gAAATG‑GGCGCGGAAGCCGCAGGT‑A‑CGCTTGTCGTCCTGAATGAATATAGCCAGCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:153985/238‑1 (MQ=255)
                                                                                                                                                      gaatccgcccgccaaGCCGCGGGT‑ACCCC‑TGTCTGCCTGATTGAATAGAGCCATCCCGCCTGCCCTCCCGAATGCGATGCCAGGCC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  2:1080696/225‑1 (MQ=255)
                                                                                                                                                         tG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCCTGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGAACAACGCACGGTAAAACGCGCCt                                                                                         >  1:410728/1‑239 (MQ=255)
                                                                                                                                                                   aaGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCCCTCATCCCACAGACACCGGTa                                                                                                                              >  2:1028689/1‑193 (MQ=255)
                                                                                                                                                                   aaGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTa                                                                                                                              <  1:1028689/193‑1 (MQ=255)
                                                                                                                                                                    acccGCAGTT‑TGCGC‑TGTCTGCTTGATTAAATATCGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgc                                                                                <  1:1058759/236‑1 (MQ=255)
                                                                                                                                                                       cccgcgggttttcctttttGCCTGTTTGAATATGCCCACCCCCCCTGCCCCCCTGAATGGCTTCACATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                               <  2:1031109/220‑1 (MQ=255)
                                                                                                                                                                                         tGCCTGATTGAATATAGCCATCACGCCTGCCTCACGGAATGCGATGCCATGAC‑GGAGCACTATTACAGCATGAGGGATTACGCTCTGCAGAACCCTCAATCCAGCGCCATTATGAGCATTATTGATTAAGGATACTCATGATGGAACAGTCAATCACCACACAGACACAGGTACTTCAACCGACCGCACAAGGCACGGTAAAACGAGCCTTAACACTGCTTGACGTACACCTGCCCGaa                                                            >  1:78435/1‑239 (MQ=255)
                                                                                                                                                                                         tGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGaa                                                            >  1:964827/1‑239 (MQ=255)
                                                                                                                                                                                                     tatgGCCATCCCCCCTGCCCCACCGAATGGGTTGCCA‑GACGGGAGCCCTTTTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCa                                                <  2:292103/239‑1 (MQ=255)
                                                                                                                                                                                                       tggCCATCACGCCTGCCTCCCCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                                              <  1:980815/237‑1 (MQ=255)
                                                                                                                                                                                                       aaGCCACCCCGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                                              <  1:230645/238‑1 (MQ=255)
                                                                                                                                                                                                           cctcccgtctgcctCACCGAAGGCGGTGCCATGAC‑GGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                                          <  1:21875/231‑1 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCCCTCATCCCACAGACACCGGTACTTCCACTGACCGCCCAAGCCACGGTACAACCCGCCTTAACGCTCCTTACCCGACACCGTCGCGAAAAAGGCGGGGCATTCACCTc                                        >  2:443715/1‑239 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAttaccct                                         >  1:871763/1‑233 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCt                                         >  1:711513/1‑238 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACCCTCCTTGACCCACACCT‑CCCCAAAAAAGCGGGGGATTCACcccc                                       >  2:459908/1‑239 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTAATTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTCCGCGAAACAGGCGTGGCATTCACCt                                         >  1:1043283/1‑238 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAACGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACACACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCCGCGCGAAAAAGGAGGGGCATGCAccc                                         >  2:463724/1‑237 (MQ=255)
                                                                                                                                                                                                              cACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCCCGGGAAAACCAGCCTTACCCCGGGTTACCCGAAACCTGCGGCAAACAGACCGGGCTTTCCCCTc                                        >  2:606303/1‑238 (MQ=255)
                                                                                                                                                                                                                      cctCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCCCGAAACAGGCGTGGCATTCACCCCCACCCAg                                 >  2:428596/1‑237 (MQ=255)
                                                                                                                                                                                                                      cctCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACGGACCGCCCAACGCACGGTCAAACCCGCCTTAACCCTGCTTGACCGAAACCTCCGCGACAAAGACTGGCCATTAACCTCCACCCAggt                               >  1:1055030/1‑238 (MQ=255)
                                                                                                                                                                                                                                tGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAAAGCTGCTTGACCGACCCCTGCGCGAAACAGGCGGGGCATTCACCTCCACTCAGGCCGCCCGGGAc                     >  2:506109/1‑239 (MQ=255)
                                                                                                                                                                                                                                       cccttgcgggAGCACTTTTTCCGCCTGCGGGATTTCGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACCCCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgctg               <  2:253788/231‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcgcg                                                                                            <  1:784525/151‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcgcg                                                                                            >  2:784525/1‑151 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGGGGCATTCCCCCCCCCCCAGCCTGCCCGGGACTGgctagaacggaaaaa    >  2:294101/1‑227 (MQ=255)
                                                                                                                                                                                                                                                   cACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                               >  1:215885/1‑163 (MQ=255)
                                                                                                                                                                                                                                                    aCTTTTACCGGCTGCGGGATTACCCTCTGCAGCACCCTGAATCCCGCCCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATgg  <  1:599883/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
TTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTA‑TCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGA‑CGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATG‑AGCGCGGAAGCCGCAGGT‑ATCGC‑TGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC‑GGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGG  >  NC_000913/2776067‑2776543

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