Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,943 T→G 22.9% S138S (TCA→TCC)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9430TG22.9% 149.2 / 5.1 79S138S (TCA→TCC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/37);  new base G (12/10);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTg                                                                                                                                                                                                                                             <  1:52679/238‑1 (MQ=255)
gtTCCCAGCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTg                                                                                                                                                                                                                                             <  2:1099742/238‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATTGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  2:110450/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCTCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:950384/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:611882/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:529084/239‑1 (MQ=255)
    ccATCCGCTGGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTTCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  2:961014/239‑1 (MQ=255)
     cATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgat                                                                                                                                                                                                                                       <  2:635802/239‑1 (MQ=255)
       tCGGCGTGAGGGAAAGCCGACGCACCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgatag                                                                                                                                                                                                                                     <  1:249339/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  1:1075988/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:211184/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:1220634/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:1023546/239‑1 (MQ=255)
           cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGAGAAAGTGACCAATTTTCTTTCCAACGGGTC‑GGGGGCGTTGGC‑AAACTC‑GGCCTGCGGTGagcggtagcgat                                                                                                                                                                                                                                 >  1:414676/1‑239 (MQ=255)
           cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGCTAACACCCAGAAGACGAGTAAATTGCTTTTCAGCGCGGA‑GGGGGCGTGCGC‑ACCGCC‑GCGCCGCCGGGatcggtagcgag                                                                                                                                                                                                                                 >  2:548900/1‑238 (MQ=255)
            gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                                <  1:182974/239‑1 (MQ=255)
                gggAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTAc                                                                                                                                                                                                                             <  1:992390/238‑1 (MQ=255)
                     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACCAATTTGCTTTCCAGCGTGTCGGAGGGCTGGGCC‑AAC‑TC‑GGCCTGCGGTGAGCGCTAGCGCTGTACCActt                                                                                                                                                                                                                        >  2:897698/1‑238 (MQ=255)
                     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTAGACTAATTTGCTTTCCAGCGTGTC‑GAGGGGTTGCGC‑CAAGTC‑GGCCTGCGGGGAGCGAGAGCGATGTAcccattt                                                                                                                                                                                                                       >  1:148373/1‑234 (MQ=255)
                                    tCGTATTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGA‑CGTTTAACTATGGTTTCTTTTACTACGAGAAAAAAGCGGCGACGGGTGCCAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAGGGCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGTTGGCGC‑AACGCC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCG‑cc                                                                                                                                                                                                         <  2:270343/237‑1 (MQ=255)
                                            gcggtggCCCACAAACCAGGATTCCGGTTGCAGA‑CGTTTAACGATGTTTTCTTTTACTTCGCGAAAAAAAGGGCGACGGGTCCCAAAATCCACCAG‑ATGGAAGCCCGACATACAAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGaa                                                                                                                                                                                               <  2:148373/233‑1 (MQ=255)
                                             cTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCCGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGGGATGGACCATTTTACCGATAACCG‑CCAGCAAAAGa                                                                                                                                                                                                >  1:1209956/1‑237 (MQ=255)
                                              tGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GGGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTGATAACCC‑CCCGCACAAGAACt                                                                                                                                                                                             >  2:428772/1‑239 (MQ=255)
                                               ggTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTGCTTTCCAGCGTGTC‑GGGGGGGTGCGC‑AAACTC‑GGCCTGCGGGGGGCGATAGCGCTGTCGCAATTTACTGATAAACGCCCACCAAGGGAACt                                                                                                                                                                                             >  2:44251/1‑239 (MQ=255)
                                                  gCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGGACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGACAAATCGACTAATTTTCTTTCCCGCCTGGT‑GAGGGCTTAAGC‑AAAGTC‑CGCCTGCGGTGAGCGAAAAGGGTTGACCCTTTTCCTGAGACACG‑CCAGCAAAGGGAATTcgc                                                                                                                                                                                         >  2:649585/1‑237 (MQ=255)
                                                  gCCCACAAACCAGGATTCCTGTGGCAGA‑CGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTcc                                                                                                                                                                                          <  1:1212552/238‑1 (MQ=255)
                                                       caAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTTCAACGGGGC‑GAGGGCCTGCGC‑AACGCC‑GGCCCGCGGTGAGCGATAACGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCCGAgg                                                                                                                                                                                    >  1:136552/1‑239 (MQ=255)
                                                         aaCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCAAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTACCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                  <  1:428772/239‑1 (MQ=255)
                                                          aCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGTGACTAATTTGCTTTTCAACCGGTC‑GAGGGGCTGCGCAAACGTG‑GGCC‑GCGGGGGGTGGAAGAGGTGTACCATTTAACTGATAACCC‑CCAGGAAAGAAACATCCCAGAGGaag                                                                                                                                                                                 >  2:934169/1‑237 (MQ=255)
                                                               gATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTTCAACGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGGCCGCGGGGGGCGATAGCGGTTTACCATTTCACCGATAAACG‑CCCGCAAAGGAAATTTCCCGAGGGTGACtt                                                                                                                                                                             >  2:519363/1‑238 (MQ=255)
                                                                aTTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                            <  1:484931/238‑1 (MQ=255)
                                                                 tttccgTTGCAGGCCGCTTAAGGATGGTTTCTTCTACGTCGCGAAAAAACCGGCGACGGGTGCAAAAAGCCATAAG‑ATGGAAGGCCGAAATATCAAGACCGGGGTTTAACGCCGAAAAGTCGATTAATTTGTTTCCCAGCGTGTC‑GAGGGCTTGCGC‑ACCGTC‑GGCCTGCGGTGGGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                            <  2:784158/233‑1 (MQ=255)
                                                                    cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTTCTTTTCAGCGGGTC‑GAGGGCTTGCGC‑AACGGC‑GGCCTGCGGTGAGCGGTAACGATGTAACATTTTACTTATAAACG‑CCCGCAAAGGAACTTCCCCGAGGGTGACTTTACGc                                                                                                                                                                        >  1:532900/1‑238 (MQ=255)
                                                                    cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAGTCGACTCATTTTCTTTCCCGCGGGTC‑GGGGGGTTGCGC‑AACGCC‑CGCCCGCGGGGGGGGGTAGCGGTGTACCATTTTCATGATAAACC‑CCCGCAAAAGGACATCCCAGGGGGAGACtttccgc                                                                                                                                                                        >  2:182688/1‑233 (MQ=255)
                                                                                           ttttCTTGTACTGCGGAAGAAAAAGGGCGAGGGGTGCCAAAATCCA‑CAGAAGGGAAGCGCGACATATCAAGACGGGCGTTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGA‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAatat                                                                                                                                                  <  1:967773/237‑1 (MQ=255)
                                                                                                 tGTCCTTGGGGAAAAAAGCGGGGACGGGGCCAAAAACCCATCAG‑ATGGAAGCCGGACATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTGGCTTTCCAGCGTGTC‑GAGGGCTTGCGA‑ACCGTC‑GGCCTGCGGGGAGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                                                                                                           <  2:961553/239‑1 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCAATGATAACCG‑CCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGc                                                                                                                                                                 >  1:1235500/1‑210 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCCGCGGTGGGCGATAGCGATGTAACATTTCACTGCTAAACG‑CCAGCAAAAGAACTTCCCCCAGGATGACTTCACGCCACCGGCCGCTTAAACGCATATCCAACTTTCCCttt                                                                                                                                    >  1:1001106/1‑238 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGACAACAAATTGCTTTACAGAGGGTCGGAGGGC‑TGCGC‑GAAATC‑GGGCTTCAGGGGGCGATAGACATTTCACATTTCCTCGCAAAACG‑CCCGCCAAGGGACTTACCAGAAGAAGGGTTCCAGCCACCGGGCGCCTAAACCGATATTCCACTTGCCACTa                                                                                                                                    >  2:210005/1‑239 (MQ=255)
                                                                                                                aaccGGCGAGGGGGGCAAAAACCCATCAG‑ATGGAAGCGCGACATACCAAGACCGGCGGTAAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                            <  1:296234/236‑1 (MQ=255)
                                                                                                                              gcaaaaacacATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCGGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                             <  2:532900/230‑1 (MQ=255)
                                                                                                                              gacAAAACCCATAAA‑GTGGAAGGCCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTCTTTTCCAGCGTGTC‑GGGGGTTGGCGC‑ACCGCC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                             <  2:849192/237‑1 (MQ=255)
                                                                                                                               ccAAAATCCATCAG‑AGAAAAGCGGAAAATATAAGGACCGGCTGTTAACGCGAATAAGTCGACTATTTGGCTTCCCGGCGTGTG‑TGGGGCTTGCGC‑ACCGTG‑CGCCGGCGGTGGTCGAGAGCGAGTTCCCATTTCACTGATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                             <  2:943213/238‑1 (MQ=255)
                                                                                                                               ccAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGGCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                             <  1:519363/238‑1 (MQ=255)
                                                                                                                                         tCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTGATAACCG‑CCCGCCAAGGAACTTCCCAAGAGATGACTTCACGCCAACGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCATTGGACACAGGGGCTTG‑TTCCGGGGTaaa                                                                                                   >  2:572585/1‑238 (MQ=255)
                                                                                                                                              aTGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTg                                                                                                                                                                                                                                             <  2:950375/98‑1 (MQ=255)
                                                                                                                                                ggAAGCGCGAAAGATCAAGACCGGGGGTTAACGCCGAGAAGCCGACTAATTTGCTTCCCAGCGTGTC‑GAGGTCTTGCCC‑AACGTC‑GGCCGGGGGTGAGCAATAGCGATGTACCTTTCCCCTGAAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGaa                                                                                             <  2:1209956/238‑1 (MQ=255)
                                                                                                                                                        aaCATATCAAGCCCGGGGGTTAACGCCGGGAGGTCGACTATTTTGTTTCCCAGCGTGTC‑GGGGGTTTGCGC‑AACGGC‑GGCCTGCGGGGGGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  2:289057/238‑1 (MQ=255)
                                                                                                                                                         aCATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  1:44251/238‑1 (MQ=255)
                                                                                                                                                            tatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTTCTTTCCACCGTGTC‑GAGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTGATAACCG‑CCCGCAAACGCACTTCCCCGAGGATGACTTCCCCCCACCGGCCGCTTAAACGAATATCCCGCCTGCCCTTATTCTTGTACACGGTGGCTTGTTTCCGGGTTAAA‑CGGACAACGCGGAACCac                                                                                >  2:689178/1‑238 (MQ=255)
                                                                                                                                                               cAAGCCCGGGGTTTAACGTCGAAACGTCGGCTAGTTGGCTTCCCAGGGTGTG‑GAGGGTTTGCCC‑AAGGTC‑GGCGCGGGGGGACGGATAGCGTGGTCCCATTTCTCGAATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                                                             <  2:1226803/239‑1 (MQ=255)
                                                                                                                                                                      gggggTAAACGCCGGGAAGGCGAATAATTTGCTTTCCAGCTGGTC‑GAGGGCATGCCA‑ACCGTC‑GGCTGGCGGTGACGAATAGCGATGTACAATTGCCCTGAGAACCG‑CCAGCAAAGGAACTGCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  2:157560/236‑1 (MQ=255)
                                                                                                                                                                       gcggTTTAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCGTTGGACACGGGGGATTG‑TTCCGGGGTAAAGCGGAACTCCGGGAAACAGTTaaaaaaaac                                                                     >  2:686493/1‑232 (MQ=255)
                                                                                                                                                                       gcggTTAACGCCGAGAAGTCGACTAATTTTCTTTCCAGCGTGTC‑GGGGGCTTGCGC‑AACGGC‑GGCCTGCGGGGAGCGAGAGCGGTGTACCATTTTACTGATAAACG‑CCCGCCAAGGAACTTCCCCGAGGGTGACTTTCCCCCACCGGCCGCTTAAACGAAAATACACCTTGCCCTTATCGTTTGACACGGTGACCTG‑TTTCGGGTTCAAGCGGAACTCCCGTAAACCGTTACGATAGt                                                                      >  2:631283/1‑238 (MQ=255)
                                                                                                                                                                        cggTTACCGCCGAGAAGTCGACTAATTTGCTTCCCAGGGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGAGATTTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                     <  1:740722/238‑1 (MQ=255)
                                                                                                                                                                        cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCGGCGTTGC‑GTGTGCTTGCGA‑ACCGGC‑TGCGTGCGGGGATCGATAGCGATGTCCCATTTCCCTAAAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                     <  2:136552/238‑1 (MQ=255)
                                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                        >  2:379034/1‑239 (MQ=255)
                                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                        >  2:116645/1‑239 (MQ=255)
                                                                                                                                                                                      aaGTCGACTAATTGGCTTTCCAGCGTGTC‑GGGGGCTTGCCA‑ACCGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                      <  1:901034/239‑1 (MQ=255)
                                                                                                                                                                                       agtcatagaaTTTGTTTCCCACTGTGCC‑GAGGGTTTGCGC‑CCCGTG‑GGCCTGCGGGGAGAGTTAGTGATGTACCATTGCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                      <  1:1158402/230‑1 (MQ=255)
                                                                                                                                                                                        gTCGGCTAATTTGCTTTCGTTGGTGGG‑GAGGTCTTGCGC‑ATGGCC‑GCCCTGCGTGGAGGGAGAGCGATGTCCTTTCTCATTGATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                     <  2:379445/238‑1 (MQ=255)
                                                                                                                                                                                         tCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                   >  2:1218794/1‑239 (MQ=255)
                                                                                                                                                                                            aCTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                                                >  1:204556/1‑239 (MQ=255)
                                                                                                                                                                                                   tgcttccgggCGTGTG‑GGGGGCTTGCG‑‑AACGGCGGGCGAGCGGTGGGGGATGGCGGTGTTCCATTTCGCTAATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                         <  2:20286/230‑1 (MQ=255)
                                                                                                                                                                                                       ttCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                        >  1:707889/1‑236 (MQ=255)
                                                                                                                                                                                                          cAGTGTGTC‑GGGGGTTTCGGC‑AATCGC‑GGCCTGCGTGGGCCAATAGCTATGTACCATTTCACTTATATCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                  <  1:865714/239‑1 (MQ=255)
                                                                                                                                                                                                              gtgtC‑GAGGGTTTGCGC‑AACGCC‑GGCCGGCGGGGAGCGATAGCGATGTACAATTTCACAGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                               <  1:416752/238‑1 (MQ=255)
                                                                                                                                                                                                              gtggC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                               <  1:1218794/238‑1 (MQ=255)
                                                                                                                                                                                                               tgtC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                             >  1:819315/1‑239 (MQ=255)
                                                                                                                                                                                                                           tGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  1:731040/1‑239 (MQ=255)
                                                                                                                                                                                                                            gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  1:1091819/1‑238 (MQ=255)
                                                                                                                                                                                                                            gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  2:249316/1‑238 (MQ=255)
                                                                                                                                                                                                                                         gCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct       >  2:334793/1‑238 (MQ=255)
                                                                                                                                                                                                                                           cggCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:1152536/237‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGGGGTGAGCGATACCGTTGTGCCATTTCACTAATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:680435/239‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCGGTGGGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:182688/239‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCGGTGAGGGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc     <  1:1117046/238‑1 (MQ=255)
                                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  2:678394/1‑238 (MQ=255)
                                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  1:467246/1‑238 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177

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