Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 433,649 T→G 21.5% G65G (GGT→GGG)  ribD → fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG21.5% 111.5 / 7.4 92G65G (GGT→GGG) ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (24/41);  new base G (14/7);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 2.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

AAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433414‑433883
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
aaGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCg                                                                                                                                                                                                                                            <  1:1110967/239‑1 (MQ=255)
     tttGGCGAAGAGATCGCGCGCCTGGAGGCCTAAGCCGTGCAGACCGAGTATGCCGTGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAAGGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGc                                                                                                                                                                                                                                       <  2:372900/239‑1 (MQ=255)
     tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCg                                                                                                                                                                                                                                        >  1:1255105/1‑238 (MQ=255)
     tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTACTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGGGCGACCg                                                                                                                                                                                                                                        >  1:837203/1‑238 (MQ=255)
       tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACC‑CAAAGTGCGAGCGcc                                                                                                                                                                                                                                      >  2:1037516/1‑238 (MQ=255)
       gggggAAGGGATCGGGCGGCTGGAGGACTAGGCCGTGCAGGACGATGATTCAATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGCCGTTTTACCACGCATACCAGCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGcc                                                                                                                                                                                                                                      <  2:1347344/234‑1 (MQ=255)
        ggCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGGAAAGCC‑CAAGGTGGGACCGCCTa                                                                                                                                                                                                                                    >  2:395726/1‑239 (MQ=255)
            aaGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                                 <  2:643917/238‑1 (MQ=255)
             aggggccgcgcgCCTGGAGCACTAAGCCGTGCAGGAGGAGTATTACATGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGGCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCGCCAACGTGCGGGTAAACCACATGCCGAAGTCCACGCTTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  2:664418/233‑1 (MQ=255)
             agagtgcgcgcgGCTGGAGTACTAAGCCGAGCGGGAACAGAATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACTGTTCACCACGCAACCCAACCCGAATGGCGGGTGGCTCATTGTCAAAGATGGCGAAATTGTCGTTGAGGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCTTTGCGTAGGGC‑GGGTAAAAAAGCA‑AAAGGTGCGACCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  2:200359/233‑1 (MQ=255)
             agagATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  1:552655/239‑1 (MQ=255)
             agagATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  1:1299057/239‑1 (MQ=255)
             agagATCGCGCGCCTGGAGGACTAAGCCGGGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  1:768998/239‑1 (MQ=255)
             agagATCGCGCGACTAGATGAATGAGCCGTGCCGGATGAGTATTCCACGGCGCGTGCGACAAAGCTGCCGCAACGAGGACGTTTTACCACGCATCCCAACCTGAATGTCGGGTGCGTCATTGTCAAAGAGGCCGAAATTGTCGGTGAAGTTTACCACCAACTTAGGGGTAAGCCACATGCCGAAGTCCACGCGTTGCTTATGGC‑GGGTAAAAAAGCC‑AAAGGTGCGCCCGCCTATGT‑c‑a                                                                                                                                                                                                                             <  2:255908/239‑1 (MQ=255)
                   cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCg                                                                                                                                                                                                                       >  1:56810/1‑239 (MQ=255)
                      gcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACAGGGCGCGGGCGCTAAAGTTGGCGCAACGAGGACGTTTTACCACGCACCCCACCCCGAATGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTTTGGG‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGaa                                                                                                                                                                                                                     <  2:809054/238‑1 (MQ=255)
                                  cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGATGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACCCTCGACACCTGTAGccct                                                                                                                                                                                                         >  1:1062342/1‑236 (MQ=255)
                                   tAAGCCGGGCAGGACAAGTTTGACATGGCGCGGGCGAGAAGTGTGGCGCAACGGGGACGTTTTACCACGCATCCCAGCCCGAATGCCGGGTGCGGTATGGCCAAAGAGGGCGAAAGTGTCGGTGAGGTTGACCCCCAGGGGGCGGTTGACCCACATGCCAAACTACCCGCTTTGCTTATGCC‑GGGT‑AAAAAGCCAAAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCcatca                                                                                                                                                                                                       <  2:17498/239‑1 (MQ=255)
                                     aGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATg                                                                                                                                                                                                     >  1:333977/1‑239 (MQ=255)
                                        cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                                  <  2:485248/239‑1 (MQ=255)
                                            cAGGACGAGTATTACAGGGCGGGGGCGCTAAGGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGGCAAAGATGGCAAAATTGTCGGTGAAGGTTCCCACAAACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                               <  2:310427/238‑1 (MQ=255)
                                              ggACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                                                             <  1:426203/238‑1 (MQ=255)
                                                 ggAGTATTAAAGGGCGGGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTGGGTGAAGGTTCCCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                          <  2:289904/237‑1 (MQ=255)
                                                  gggTATTACATGGCGCGGGCGCTAAAGCGGGCGCAAGGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGCCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc                                                                                                                                                                                        <  1:221695/237‑1 (MQ=255)
                                                  gaaagttACAGGGCGCGGGCGCGAAAGCGGGCGCAACGAGGACGTTTTACCACGCATCCAAACCGGAAGGTCGGGTGCGTCAGTGTCAAAGAGGGCGAAATGGTCGGTGAAGGTTACCACCAACGGGGGGGTGAACCACAGGCCGAAGTACACGGGTGGCGTAGGGC‑GGGTGAAAAAGCA‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                          <  1:73938/232‑1 (MQ=255)
                                                   attATTACATGGCGCGGGCGCTAAAGCGGGCGCACCGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                    <  1:1132474/224‑1 (MQ=255)
                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTCTGT‑C‑ACACTCGAAACCTGGAGCCATACTGGGCGTACGCCACCgc                                                                                                                                                                                      >  1:517212/1‑238 (MQ=255)
                                                       ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTCGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGGGCGACCGCCCATTT‑C‑ACACTCCCACCCTGTGACCAACACGGTCGTAACCCACCgcgtc                                                                                                                                                                                   >  2:806366/1‑235 (MQ=255)
                                                            tGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGAC‑AAAGAGGCGGCCGCCTATGTCC‑ACAC‑CGAACCCTGGAGCCCACTTGGTCGAACGCCCCCCTggcgtggc                                                                                                                                                                              >  2:844166/1‑232 (MQ=255)
                                                             ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCC‑AAAGG‑GCGAACGCCCATGT‑CAACAC‑CGAACCCCGTAACCATCAAGGGCGTCAGCCAcccctctcttata                                                                                                                                                                             >  2:1132474/1‑228 (MQ=255)
                                                             cgcgggcgcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGGGGGACCGCCCATGT‑C‑ACACCCGCACCCCGTAGCCATCATTGTCGTACCCCCCCGGGCTGGGAc                                                                                                                                                                              >  2:1207346/8‑238 (MQ=255)
                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACA‑AAGGGGGGGCCCGCTTATGC‑G‑ACATACGACCCTGGTAG‑CAACCTGGGCGGACGCCCCCGTGCtgacgagc                                                                                                                                                                            >  2:1280436/1‑232 (MQ=255)
                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGc                                                                                                                                                                            >  1:676894/1‑238 (MQ=255)
                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGc                                                                                                                                                                            >  1:707308/1‑238 (MQ=255)
                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGGGCGACCGCCCATGT‑C‑ACACTCGAACCCTGGAGCCATCCTTGGCGTACGCCACCGTGCTGTGACGc                                                                                                                                                                            >  1:321368/1‑238 (MQ=255)
                                                                gcgggggtaaaaGCGGGGGAAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTAAAAGATGGCGAAATTGGCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGG‑GGGTGAAAAGGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                                                           <  2:529480/229‑1 (MQ=255)
                                                                  gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGCC‑AAAGGGGCGACCGCCCCCGT‑C‑CCCCCCCCACCCTCTAGCCAACATCGTCGTACGCCACCgcgcggggcccgtcc                                                                                                                                                                         >  2:73938/1‑227 (MQ=255)
                                                                     cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACC‑AAAGGTGCGACCGCCTCTGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTCTGACGCACTCAt                                                                                                                                                                      >  1:12905/1‑238 (MQ=255)
                                                                     cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCAt                                                                                                                                                                      <  1:225989/238‑1 (MQ=255)
                                                                         aaaGCTGGCGCAACGAGGACGTTTTCCCACGCACCCCACCCCGAATGGCGGGTGCGTCATTGACAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTAAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccg                                                                                                                                                                 <  2:78611/239‑1 (MQ=255)
                                                                                 cgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAGCC‑AAAAGGGCGACCGCCCATGT‑C‑ACAATCGCACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGTCGCACTCATCCCCCCCGGGGTa                                                                                                                                                         >  2:381696/1‑239 (MQ=255)
                                                                                              ttttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGtt                                                                                                                                             >  1:369817/1‑238 (MQ=255)
                                                                                                ttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGTAAAAGCC‑AAAGGCGCGGCCCCCCCTTT‑C‑GCCCTCGTAACCTGTAGCCATCATCGCCGTACCCCACCCGGCTGTGACGCAGTCCTCCCCCCTGGCGCATCGCGCGTCGGTGc                                                                                                                                           >  2:104903/1‑238 (MQ=255)
                                                                                                      cGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAAGAGCGACCGCCTATTG‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGCGCTGGTACGCCCCCATCCCCGCTGTCGTAGCGCGCCGGGGTGcccctgtg                                                                                                                                    >  2:293739/1‑233 (MQ=255)
                                                                                                        cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGGCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCCCCATGCa                                                                                                                                  >  1:250912/1‑239 (MQ=255)
                                                                                                        cATCCCAACCCGAAGGTGGGGTGCGTCATTGTCAAAGATGGCGAAATGGCCGGGGAAGGTTCCCACCAACGTGGGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                  <  2:56810/239‑1 (MQ=255)
                                                                                                         aTCCCACCCCGAATGTCGGGTGCGTCATTGTAAAAGATGGGGAAATTGTGGGGGAAGGTTACCACAAACGTGCGGGTGAACCACATGCGGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                  <  2:125393/238‑1 (MQ=255)
                                                                                                           cccAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAACC‑AAAGGTGCGGACGCCTATGT‑C‑ACCCTCGAACCCTGTAGCCATCATGGTCGGAAGCCACCGCGCTGTGGCGCCCTCACCCCCCCTGGCGTGTGGCGCGTGGGTGCCTCGATggcc                                                                                                                                 >  2:374901/1‑234 (MQ=255)
                                                                                                               aCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACC‑AAAGGTGCGGCCGCCTATGT‑C‑CCCCCCCGCCCCCTTAGCCATCATGGTCGTACCCCCCCGTGCTTTGACCCACTCCTCTCCCCTGGCGGAGCGCGCGCGGTTGCCTCCCTGCCAGATccc                                                                                                                           >  2:120947/1‑238 (MQ=255)
                                                                                                                cccGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTa                                                                                                                          >  1:1300201/1‑239 (MQ=255)
                                                                                                                cccGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAAGGGGGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCCTACTGGTCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCCGGCGTAGCGCCCGGTGGTGCCTCCATGCAAGATCCTa                                                                                                                          >  2:1032343/1‑239 (MQ=255)
                                                                                                                 ccGAATGTCGGGTGCTTCATTGTCAAAGATGGCAAAATTGGCGTGGAAGGTTACCCCCAACGGGCGGGTGAACCAAATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                     <  2:356622/179‑1 (MQ=255)
                                                                                                                   gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAAc                                                                                                                        >  1:719907/1‑238 (MQ=255)
                                                                                                                     aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGGGCGACCGCCCATGT‑T‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGGCGCCCTCATCGCCGCTGGCGTAGCGCGCGGGGGTGGCTCGGTGCACGATCCTAAccc                                                                                                                      >  2:692678/1‑238 (MQ=255)
                                                                                                                           ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGGCATCATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                                                           >  1:356622/1‑179 (MQ=255)
                                                                                                                           ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTc                                                                                                               >  1:472585/1‑239 (MQ=255)
                                                                                                                              gCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                             <  1:1032343/238‑1 (MQ=255)
                                                                                                                                  aaTTGTCAAAGATGGCGAAATTGGCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                        <  1:1348047/238‑1 (MQ=255)
                                                                                                                                       taaaaGAGGGCGAAATTGTCGGTGAAGGTAACCGCCAACGTGGGGGTGAACCACAGCCCGAAGCACACGCGTTGCGTATGGC‑GGGTGAAAAAGCA‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                             <  1:1194411/227‑1 (MQ=255)
                                                                                                                                         aaaGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACCCTCGAAACCTGTAGCCATCACGGTCGGACGCCACCGGGCTGTGACGCACTCATCGCCGCCGGGCGAGCGCGCGTGGGTGCCTCGCTGCAAGATCCCAACCCCCACGGCGCGCGGCGGGGACt                                                                                                 >  2:87965/1‑239 (MQ=255)
                                                                                                                                              tGGCAAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                              <  1:1302762/237‑1 (MQ=255)
                                                                                                                                                 cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAACC‑AAAAGCGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGACAACATGGTCGTACGCCACCGGGCGGTGACGCACCCATCGCCGCTGGCGTAGCTCGCGGTGTTGCCTCGATGCAAGATCCTAACCCGCACGTCCCTTGGCGGGGaactgtccctt                                                                                         >  2:1194411/1‑229 (MQ=255)
                                                                                                                                                 cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGTGCGACCGCCCATGT‑C‑ACACCCGAACCCTGTTGCCATCATGGTCGTACGCCACCGTCCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCggggggcttatcg                                                                                           >  2:1302762/1‑228 (MQ=255)
                                                                                                                                                 cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAAGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCCTCATGGGCGTACGCCCCCCTGCTGTGACGCACTCACCGCCCCTGGCGGAGCGCGCGTGGTTGCCCCCATGCAAGATCCTAAACCCCCGGTCGCCGGGCGCGGGCTTTccccc                                                                                          >  2:970611/1‑236 (MQ=255)
                                                                                                                                                                   ttcccaccaCCGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCAt                                                                        <  1:898030/234‑1 (MQ=255)
                                                                                                                                                                       caccaACGTGGTGGTAAACCACATGCCGAAGTACACGTGTTGCTTATGGG‑GGGTGAAAAGGCC‑AAAGGTGCGGCCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGt                                                                  <  1:816540/239‑1 (MQ=255)
                                                                                                                                                                           accGTGCGGGTGACCCCCATGCCGAGGTACCCGCGTTGCTTTGGGG‑GGGGGAAAAAGCC‑AAGGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGc                                                              <  2:116566/237‑1 (MQ=255)
                                                                                                                                                                             cgggggggTAAACCACAGGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCA‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCa                                                            <  2:584103/234‑1 (MQ=255)
                                                                                                                                                                             cgggCGGGTGAACCCCAGCCGGAAGTCCCCGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCa                                                            <  2:1300201/236‑1 (MQ=255)
                                                                                                                                                                             cGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCa                                                            <  1:87965/239‑1 (MQ=255)
                                                                                                                                                                                 cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATCGCCCCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGGCTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACgg                                                         >  1:531646/1‑238 (MQ=255)
                                                                                                                                                                                  gggtacaCCACTGGCCGAATTGCACGCTTTGGGTATGGG‑GGGTAAAAAGGCA‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGc                                                        <  2:517212/232‑1 (MQ=255)
                                                                                                                                                                                         ccACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAAACGAGAACGCCTATTT‑T‑ACACTCGAACCCTGGAGCCATCATGGTCGTAAGCCACCGAGCCGGTTCGCCCTCAACGCCGCTGGGGTAACGCTCCGGGTTTCCACCCTGCAAGATCCTACCCCGCAGGGCGCTTGGCGGGGGCTTTTCCGTCTGCAACAGGCTAGGTGTTG‑CGTCAGGCACGGCCtgctgat                                                >  2:134069/1‑239 (MQ=255)
                                                                                                                                                                                              tGCCGAAGTCCACGCGTTGCGGATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                      <  1:663153/181‑1 (MQ=255)
                                                                                                                                                                                                ccGAAGTCCCCGCGTTGCGTATGGG‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGa                                          <  1:806366/238‑1 (MQ=255)
                                                                                                                                                                                                 ggagctACACCCGGTCCGTATGGG‑GGGTGAAAAAGCC‑AAGGGTGCGGCCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg                                        <  2:1255105/234‑1 (MQ=255)
                                                                                                                                                                                                    aGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                                     <  1:491760/239‑1 (MQ=255)
                                                                                                                                                                                                        cacGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGAGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGCGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATACTAACCCGCAGGGCGCTGGGCGCGGGCCTTAcc                                                                                            >  2:663153/1‑181 (MQ=255)
                                                                                                                                                                                                        cacGCGTTGCGTATGGC‑GGGGGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCAACATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTTCCTCGATGCAAGATCCTAAACCCCAGGTCCCCGGGGGTGGGCTTTTCCGGCCGCAACAGGCT‑GGCATTTGCGTCAGACACGGCCTGGTGATGAGGGAAGCCGGGCa                                 >  2:98318/1‑239 (MQ=255)
                                                                                                                                                                                                         acGGGTTGCGTAGGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                <  1:395726/239‑1 (MQ=255)
                                                                                                                                                                                                         acGGGTGGCGTATGGG‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                <  1:104903/239‑1 (MQ=255)
                                                                                                                                                                                                          cgcgTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt                               <  1:1037516/239‑1 (MQ=255)
                                                                                                                                                                                                             gTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAAGCACCATTGTCGTACGCCACCGTGCTGTGACGCACGCATCGCCCCCGGCGTAGCGCGCGTGGGTGCCTCCATGCCAGATCCTAACCCGCAGGTCCCTGGGCGTGGACTTTACCGCCTGCAACAGGCT‑GGCATTTACGTCCGCCACGGCCTGATGTAGAGTGAAGGCGGGCAATTgg                            >  2:466806/1‑238 (MQ=255)
                                                                                                                                                                                                              ttGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGaa                           >  1:148470/1‑239 (MQ=255)
                                                                                                                                                                                                                gCGTATGGC‑GGGGGAAAAGGCC‑AAGGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt                          <  2:973757/238‑1 (MQ=255)
                                                                                                                                                                                                                 cGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                                     >  2:1309842/1‑226 (MQ=255)
                                                                                                                                                                                                                 cGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                                     <  1:1309842/226‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGG‑GGGTAAAAAAGCC‑AAAGGGGCGCCCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAg                      <  2:684928/238‑1 (MQ=255)
                                                                                                                                                                                                                            gTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtc                >  1:973757/1‑237 (MQ=255)
                                                                                                                                                                                                                            gTGAAAAAGCC‑AAAGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAACCGCCGGTCGCCGGGCGGGGGCTTTACCGGCTGCAACAGGCT‑GGCCTTGACGTCCGCCCCGGCCTTGTGGTGGGTGAAGCCGAGCAATTCAATAAAGACtttat               >  2:1082169/1‑236 (MQ=255)
                                                                                                                                                                                                                               aaaaaGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg           <  1:381696/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aaGGGGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgcgc  <  2:472585/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
AAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCC‑AAAGGTGCGACCGCCTATGT‑C‑ACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCT‑GGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433414‑433883

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: ‑