Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,939 A→G E235G (GAA→GGA)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,9390AG99.2% 399.9 / ‑5.5 122E235G (GAA→GGA) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (60/61);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 122 total aligned reads displayed.

GGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGAAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGA  >  NC_000913/3152694‑3153184
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ggAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAAt                                                                                                                                                                                                                                                   <  1:252889/250‑1 (MQ=255)
        cAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTg                                                                                                                                                                                                                                          <  1:438278/251‑1 (MQ=255)
           gCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATg                                                                                                                                                                                                                                       <  1:264661/251‑1 (MQ=255)
                     cATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATgg                                                                                                                                                                                                                             >  1:161164/1‑251 (MQ=255)
                      aTGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGt                                                                                                                                                                                                                            >  1:663502/1‑251 (MQ=255)
                      aTGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGt                                                                                                                                                                                                                            <  1:243545/251‑1 (MQ=255)
                      aTGGAAGGCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGt                                                                                                                                                                                                                            <  2:161164/251‑1 (MQ=255)
                            gTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGc                                                                                                                                                                                                                      <  1:643773/251‑1 (MQ=255)
                             tCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGcc                                                                                                                                                                                                                     >  2:138706/1‑251 (MQ=255)
                                  gacgTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACCCCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATAc                                                                                                                                                                                                                <  2:590403/251‑1 (MQ=255)
                                    cgTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCg                                                                                                                                                                                                              >  1:629659/1‑251 (MQ=255)
                                    cgTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCg                                                                                                                                                                                                              <  1:164421/251‑1 (MQ=255)
                                             gtttgttgCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCTTTCAGATGCGATGATTGGCACTGCCGTGTGAAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCATATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATaa                                                                                                                                                                                                     <  2:211224/249‑1 (MQ=255)
                                             gATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATAGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCtatata                                                                                                                                                                                                      >  2:177381/1‑250 (MQ=255)
                                              aTTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGGGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAAc                                                                                                                                                                                                    >  1:287052/1‑251 (MQ=255)
                                                                 gtgtGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGc                                                                                                                                                                                                             <  1:271159/223‑1 (MQ=255)
                                                                  tgtgGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACAt                                                                                                                                                                                >  2:602449/1‑251 (MQ=255)
                                                                    tggtgCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACAtt                                                                                                                                                                               >  2:332129/1‑250 (MQ=255)
                                                                               gccttcatTATGATCGACAACCCCTGGGCAGCCGGTGGGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGttt                                                                                                                                                                   <  2:264772/247‑1 (MQ=255)
                                                                                    catTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTc                                                                                                                                                              >  2:294540/1‑251 (MQ=255)
                                                                                     atTATGATCGACAACCCCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCa                                                                                                                                                             <  1:363101/251‑1 (MQ=255)
                                                                                      ttATGATAGACAACACCTGGCCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGAGGTTTCTATTCAAGCCGCCCCCAAATATCTGGTGGGGAATTCAGATGCGATGATTGGCACTGCCTTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTc                                                                                                                                                              <  2:377483/249‑1 (MQ=255)
                                                                                       tATGATCGACAACACCTGGGCAGTCGGTGTGCTGTTTACGGCGCTGGATTTGGGCATCGATGTTTCTATTCAAGCAGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGAAGCTACGGGGAAATGCCTAGCTGATGGTCCAGATGGTCGATGCCGATACCCCCTATATAACCAGACGTGCCCTGCGCACAGTAGGTGTGCGTTTGCGTTAAc                                                                                                                                                           >  2:298721/1‑251 (MQ=255)
                                                                                                 aacaCCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGaa                                                                                                                                                   >  1:261324/1‑249 (MQ=255)
                                                                                                  acacCGGGGCAGCCGGTGTGCTGTTTAAGGCTTGGGATTTTGGCACCGCTGTTGCTATTCAAGCCGCCCCCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGACTTGTGCAATGCCAGTTGCTGGGAGCAGTTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATTCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGc                                                                                                                                                <  2:126878/251‑1 (MQ=255)
                                                                                                       ggggCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCt                                                                                                                                           <  1:138706/250‑1 (MQ=255)
                                                                                                         ggCAGCGGGTGTGCGTTTTAAGGCGCTGGTTTTTGGCAGCGATGTTTCTATGCAAGCCGCCACCAAATTTCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGGGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGa                                                                                                                                         <  2:11406/251‑1 (MQ=255)
                                                                                                             gCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGt                                                                                                                                     >  2:244798/1‑251 (MQ=255)
                                                                                                               cGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg                                                                                                                                   >  2:412927/1‑251 (MQ=255)
                                                                                                                ggTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCCCTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGc                                                                                                                                  <  1:20466/251‑1 (MQ=255)
                                                                                                                  tgtgCTGTTTAAGGCGCTGGCTTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTg                                                                                                                                >  2:518828/1‑251 (MQ=255)
                                                                                                                   gtgCTGTTTAAGGCGTTGGTTTTTGGCATCGTTGTTTCTATTCAAGCCGCCCCCAAATATCGGGTTGGGAATTCAGAGGCGATGATTGGCCCTGCCGTGTGCAATGCCCGTGGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGa                                                                                                                               <  2:506526/251‑1 (MQ=255)
                                                                                                                     gCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGaa                                                                                                                              <  1:332129/250‑1 (MQ=255)
                                                                                                                      cTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGAt                                                                                                                                                                                                                  <  2:53727/165‑1 (MQ=255)
                                                                                                                      cTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGAt                                                                                                                                                                                                                  >  1:53727/1‑165 (MQ=255)
                                                                                                                      cTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtg                                                                                                                            >  2:265403/1‑251 (MQ=255)
                                                                                                                        gttaaaGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCGGGTTGGGAATTAAGATGCGATGATTGGCACTGCCGTGTGAAATGCCCGTGGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtgg                                                                                                                           <  1:354704/246‑1 (MQ=255)
                                                                                                                        gTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCCCCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTGGATGCCGGTCCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtg                                                                                                                            <  2:254804/249‑1 (MQ=255)
                                                                                                                        gTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                                                                                                                          >  2:151907/1‑251 (MQ=255)
                                                                                                                           aaaGGCGCTGGATTTTGGCATCGATGTTTCTATTCAGGCCGCCACCAAATATCTGGTTGGGCATTCAGAGGCGAGGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggctgg                                                                                                                       <  1:630269/250‑1 (MQ=255)
                                                                                                                            aaGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggctgg                                                                                                                       <  1:177381/250‑1 (MQ=255)
                                                                                                                                     gATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                                                                                                             >  2:486533/1‑251 (MQ=255)
                                                                                                                                         ttGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCTATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa                                                                                                          >  1:358784/1‑250 (MQ=255)
                                                                                                                                         ttGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa                                                                                                          >  1:416477/1‑250 (MQ=255)
                                                                                                                                              aTCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTgc                                                                                                    >  1:34205/1‑251 (MQ=255)
                                                                                                                                              aTCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTg                                                                                                     >  1:350653/1‑250 (MQ=255)
                                                                                                                                                 gATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                                                                 >  1:660245/1‑251 (MQ=255)
                                                                                                                                                 gATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCAGCAAGTTGcgcg                                                                                                 >  2:362178/1‑251 (MQ=255)
                                                                                                                                                    gtttttttttaggcccccccaaaaTTTCTGGTTGGGCATTCAGAGGCGATGATTGGCACTGCGGTGTGCAATGCCGGTTGCTGGGAGCGGCTACGGGAAAGTGCCTATCTGATGGGCCAGATGGGCGATGCCGATTCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAg                                                                                               <  2:287052/229‑1 (MQ=255)
                                                                                                                                                        cTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACCTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                                                          >  2:506155/1‑251 (MQ=255)
                                                                                                                                                          aTTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                                                          <  1:137105/249‑1 (MQ=255)
                                                                                                                                                              aggccgccACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCGACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                                                                    <  2:46506/249‑1 (MQ=255)
                                                                                                                                                              aaGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                                                                    <  1:486533/251‑1 (MQ=255)
                                                                                                                                                              aaGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                                                                    <  1:265403/251‑1 (MQ=255)
                                                                                                                                                              aaGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                                                                    <  1:151907/251‑1 (MQ=255)
                                                                                                                                                                      accaAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGcc                                                                            >  1:264067/1‑251 (MQ=255)
                                                                                                                                                                         aaaTATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTgg                                                                         >  2:685945/1‑251 (MQ=255)
                                                                                                                                                                          aaTATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggct                                                                                                                         >  1:354501/1‑202 (MQ=255)
                                                                                                                                                                           atatCTGGTTGGGCATTCATATGCGATGTTTGGCACTGCGGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGctcn                                                                               <  2:328360/243‑2 (MQ=255)
                                                                                                                                                                           atatCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGctct                                                                               >  1:328360/1‑243 (MQ=255)
                                                                                                                                                                           atatCTGGATGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGc                                                                        >  2:472485/1‑250 (MQ=255)
                                                                                                                                                                                  gTTGGGCATTCAGAGGCGATACTTGGCACTGCGGTGGACAATGCCGGTTGCTGGGGGCAGCTCCGGGGAAATGCCTATCGGAGGGGCAAGATGGTCGATGCATATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAgg                                                                <  2:350653/251‑1 (MQ=255)
                                                                                                                                                                                    gggtcATTCAGATGCGATGATTGGCACTGCCGTGGGAAATGCCCGTGGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTc                                                              <  2:660245/247‑1 (MQ=255)
                                                                                                                                                                                      ggCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCAc                                                            >  1:78349/1‑251 (MQ=255)
                                                                                                                                                                                                   aatgatTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAACGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGGCACGAATTCTGGAAAc                                               >  2:359728/2‑251 (MQ=255)
                                                                                                                                                                                                     tgatTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTGGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACg                                              >  1:412328/1‑250 (MQ=255)
                                                                                                                                                                                                      gatTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACgag                                            <  1:294540/251‑1 (MQ=255)
                                                                                                                                                                                                      gatTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACg                                              >  1:557632/1‑249 (MQ=255)
                                                                                                                                                                                                      gatTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACg                                              >  1:87709/1‑249 (MQ=255)
                                                                                                                                                                                                       atTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTgc                                                                                                    <  1:658845/194‑1 (MQ=255)
                                                                                                                                                                                                       atTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACg                                              >  1:207044/1‑248 (MQ=255)
                                                                                                                                                                                                         tGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACt                                         >  2:582465/1‑251 (MQ=255)
                                                                                                                                                                                                          ggcccTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACtt                                        <  1:315453/247‑1 (MQ=255)
                                                                                                                                                                                                            cACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGACTTCTGGAAACGAGACttt                                       >  1:531124/1‑250 (MQ=255)
                                                                                                                                                                                                               tGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTAc                                     >  2:357909/1‑249 (MQ=255)
                                                                                                                                                                                                                  ngtgtGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATTCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGgc                                 >  2:317814/2‑250 (MQ=255)
                                                                                                                                                                                                                   gtgtGCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCGCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGgcag                               <  1:671269/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggCAATGCCCGTTGCGGGGAGCAGTTACGGGGAAAGGCCTCTCTGAGGGGCCGGATGGTCGATGCCGATACCGGCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCGGGAAACGAGACTTTACAGGcagcag                            <  1:244798/250‑1 (MQ=255)
                                                                                                                                                                                                                       gCAATGCCCGTTGCTGGGAGCAGCTCCGGGGAAATGCCTATCTGATGGGCCAGATGGGCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCagcagc                           <  2:160132/251‑1 (MQ=255)
                                                                                                                                                                                                                       gCAATGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCagcagc                           >  2:604271/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGgcagca                             >  2:146744/1‑246 (MQ=255)
                                                                                                                                                                                                                          aTGCCCGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGAGGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGgcagca                             <  1:146744/246‑1 (MQ=255)
                                                                                                                                                                                                                            gCCCGTGGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGTAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCt                      <  2:348380/251‑1 (MQ=255)
                                                                                                                                                                                                                             cccGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGTTGGTCGATGCCGATCCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGAAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTa                     <  1:685945/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGt                                                                                              >  2:295900/1‑176 (MQ=255)
                                                                                                                                                                                                                               cGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTAtt                   <  1:113632/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGTTGCTGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTAtt                   >  2:666501/1‑251 (MQ=255)
                                                                                                                                                                                                                               cGTTGCGGGGAGCGGTTACGGGGAAATGCCTATCTGATGGGCCAGAGGTTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGt                                                                                              <  1:295900/176‑1 (MQ=255)
                                                                                                                                                                                                                                     tGGGAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCtt             >  2:642282/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggggAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCtt             <  1:602449/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gggAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCttt            <  1:642282/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gggAGCAGCTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGACAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCt              <  1:361872/249‑1 (MQ=255)
                                                                                                                                                                                                                                         ggcagcTAGGGGGAAATGCCTATCTGATGGGCCAGATGGTCGAGGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTtgtg         <  2:663502/250‑1 (MQ=255)
                                                                                                                                                                                                                                         agcagcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc                                                                                                                                                           <  1:236625/105‑1 (MQ=255)
                                                                                                                                                                                                                                         agcagcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc                                                                                                                                                           >  2:236625/1‑105 (MQ=255)
                                                                                                                                                                                                                                         agcagcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGCAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTtgtg         >  2:379474/1‑251 (MQ=255)
                                                                                                                                                                                                                                         agcagcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTtgtg         >  1:673474/1‑251 (MQ=255)
                                                                                                                                                                                                                                            agcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATTCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGcgt      <  1:317814/251‑3 (MQ=255)
                                                                                                                                                                                                                                             gcTACGGGGAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGCGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCtt      >  1:290084/1‑250 (MQ=255)
                                                                                                                                                                                                                                                aCGGGGAAATGCATATTTGATGGGCCAGATGGTCGTTCCCCATACCGCCTATATAACCAGACGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAg   <  1:598343/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGATTGTTGCCGCCGTACGCAGTGTGGTGCCGGATGCCATCATTATGATCGACAACACCTGGGCAGCCGGTGTGCTGTTTAAGGCGCTGGATTTTGGCATCGATGTTTCTATTCAAGCCGCCACCAAATATCTGGTTGGGCATTCAGATGCGATGATTGGCACTGCCGTGTGCAATGCCCGTTGCTGGGAGCAGCTACGGGAAAATGCCTATCTGATGGGCCAGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGA  >  NC_000913/3152694‑3153184

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