Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,959 C→T P253P (CCG→CCA fecR ← ferric citrate regulator FecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,9590CT97.2% 315.9 / NA 110P253P (CCG→CCAfecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (53/52);  minor base G (3/0);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516714‑4517204
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGttt                                                                                                                                                                                                                                                  >  1:475052/1‑251 (MQ=255)
gAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGtt                                                                                                                                                                                                                                                   <  2:65886/250‑1 (MQ=255)
  gTCGAACGAGACGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTATCCCCTCACCCCGTGGTTTAt                                                                                                                                                                                                                                                >  2:317963/1‑251 (MQ=255)
  gTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTTCAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCATCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTAt                                                                                                                                                                                                                                                <  1:25243/251‑1 (MQ=255)
        cGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTg                                                                                                                                                                                                                                           >  1:568885/1‑250 (MQ=255)
        cGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGa                                                                                                                                                                                                                                          <  2:76302/251‑1 (MQ=255)
         gagaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGCGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGCATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGCACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGCTCGCCGCGCAGCACGGCGTTGCGGGCACGGGATAGCGGGGCTATCACCTCCCCCCGTGGTTTATCGCTgg                                                                                                                                                                                                                                          >  1:49689/1‑249 (MQ=255)
         gagaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACGTGGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCAGCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGaa                                                                                                                                                                                                                                         <  2:369819/251‑1 (MQ=255)
          agaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCCGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGGTTATCGCTGa                                                                                                                                                                                                                                          >  1:268579/1‑249 (MQ=255)
                  aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGa                                                                                                                                                                                                                                 <  2:585251/250‑1 (MQ=255)
                              cGAAATAAGAATTATTTTCCTTACAGTGGTGCAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAAAGGGAAGCGTTTGCGCGATCACGTTCAGGATCGCATCGGTATTTTACAGCGGGAACGTCCCGCTCAGGAGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGTGGGAACGGGTTAGCGGGGCTATCACCTCACCCCGTGGTTTATCGCTGCAGCTCAGGATGTCCTTCGTccc                                                                                                                                                                                                                    >  1:567307/1‑250 (MQ=255)
                              cGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGGTTATCGCTGAAGCTCAGGATGTCCTTCGTccc                                                                                                                                                                                                                    >  1:351541/1‑250 (MQ=255)
                               gAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAg                                                                                                                                                                                                                   <  2:510253/251‑1 (MQ=255)
                                  aTAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCtt                                                                                                                                                                                                                <  1:9593/251‑1 (MQ=255)
                                     aGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGGTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGGTGTCCTTCGTCCCGCTTGt                                                                                                                                                                                                              >  2:506689/1‑250 (MQ=255)
                                           attTTCCTTACAGTGGTGAGATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGAATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACTGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTCTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCg                                                                                                                                                                                                       >  1:425778/1‑251 (MQ=255)
                                             tttCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTACTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCCGCACGCCGTTGCGGGAACGGGGTAGCGGGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTc                                                                                                                                                                                                     >  2:423974/1‑251 (MQ=255)
                                              ttCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTc                                                                                                                                                                                                     >  1:593816/1‑250 (MQ=255)
                                              ttCCTTACAGTGGTGAAATGTTTATCCAGCACCGCGTAATAGCCTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTc                                                                                                                                                                                                     <  2:516392/250‑1 (MQ=255)
                                               tCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCACCACGCCGTTGCGGTAACGGGCTAACGTGGCTATCAACTCACCCAGTGGCTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCAt                                                                                                                                                                                                   >  1:628997/1‑251 (MQ=255)
                                                ccTTACAGTGGTGAAGTGTTTATACAGTCCAGCGTATTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCAACGGAATTTTTCAGCAGGACCGGCCCGCTCAGGCGCAGCCCGGCAACGGCGGCAATGCAGCGCAGGACGCCGTCGCGGTAAAGGGTTATGGCGGCTTACGCCTCACCCAGTGGTTTATGGCTGAGGCTCAGGGTGTCGTTCGTCCAGccagtacttggttcgtt                                                                                                                                                                                                  >  2:61618/1‑235 (MQ=255)
                                                   tACAGTGGTGAAATGTTTAACCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGAGAACGTTCAGGATCGCATCGGTATTTGTCAGCGTGAACGTCCCGCTCAGGCGCAGCGCGGCAACGGCGGGATCGCAGCCCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAGCACCTCACCCAGTGGGTTATCTCTGAAGCTCAGCATGTCCTTCGTCCGGCGTGTACTCTCGTCATcc                                                                                                                                                                                                 <  2:49689/249‑1 (MQ=255)
                                                    aCAGTGGTGCACTGTTTCTCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGACGCGCCGCACGCCGTTTCGGGAACGGGTTAGCGTGGCTAACCCCTCACCCAGTGGTTTATCGCTGAAACTAAGGATGTCCCTCGTCCCGCTTGTTATCTCCTCATccgac                                                                                                                                                                                              >  2:314522/1‑248 (MQ=255)
                                                       gtggtgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCTCTGAAGCTCAGGAGGTCCTTCGTCCAGCTTGTACTCTCGTCATCCCGCGGt                                                                                                                                                                                           >  1:236158/1‑251 (MQ=255)
                                                                aTGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGGTCAGGAGCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGGACTCTCGTCATCCAGCGGTTTCACTgc                                                                                                                                                                                   <  2:236377/251‑1 (MQ=255)
                                                                       tCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCCCAGCGCAGCACGCCGGTGCGGTAACGGGTTAGCGTGGCTATCACCTCCCCCCGTGGTTTTTCCCTGGAGCTCAGGGTGTCCTTCGTCCAGCTTGTACTCTCGCCATCCAACGGTTTCACTGCGGCAAAc                                                                                                                                                                            >  2:566871/1‑251 (MQ=255)
                                                                        ccAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCTCTGGAGCTCa                                                                                                                                                                                                                                    >  2:226338/1‑194 (MQ=255)
                                                                        ccAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCa                                                                                                                                                                                                                                    <  1:226338/194‑1 (MQ=255)
                                                                         cAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTc                                                                                                                                                                          >  1:312320/1‑251 (MQ=255)
                                                                             aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTCAGCGGGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGCTTCACTGCGCCCACCTCag                                                                                                                                                                        >  2:340931/1‑249 (MQ=255)
                                                                               cgcgTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTTCGGTAACGGGTTAGCGTGGCTATCACCTCAACCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCACACTCAGAggg                                                                                                                                                                    >  2:144869/1‑250 (MQ=255)
                                                                               cgcgTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGgc                                                                                                                                                                    <  2:416216/251‑1 (MQ=255)
                                                                                       aGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACt                                                                                                                                                            >  1:587761/1‑251 (MQ=255)
                                                                                       aGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACt                                                                                                                                                            >  2:196480/1‑251 (MQ=255)
                                                                                          cTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAACGTGGCTATCACCTCACCCAGTGGGTTATCTCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCCGCGGGTTCACTGCGCCCAACTCAGAGGCGCTGAACTGc                                                                                                                                                          >  2:183647/1‑250 (MQ=255)
                                                                                               ttttAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCCGCACGCCGTTGCGGTAACGGGTTAGCGTGGATATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCCTCCACCCCCGGTTTTACTGCGCCAAACCCCGAGGCGCTGAACCGCAGGCtt                                                                                                                                                    >  2:525241/1‑251 (MQ=255)
                                                                                                   aaCGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCCGGCGCAGCCCGGCAACGGCGGGGTCGCAGCGCAGCACGCCGTTGAGGTAACGGGTTAGCGTGGCTATTCCCCAACCCAGTGGTTTATCCCTGAAGATCCGGGTGTAATTCGTCCAGCTTGAACTCTCGTCATCCCGCGGGTTCACTGCGCCAAACCCAGAGGCCCTGAAACGCAGGCCTTCAc                                                                                                                                                >  2:541742/1‑251 (MQ=255)
                                                                                                     cGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGGTCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCCAACTCAGAGGCGCTGAACTGCCGGCTTTCAccc                                                                                                                                              >  2:209077/1‑250 (MQ=255)
                                                                                                      gggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAAAGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAACGGTTTCACTGCGCCCAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                                                                                                               >  1:498300/1‑249 (MQ=255)
                                                                                                            cGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCa                                                                                                                                       <  1:116438/251‑1 (MQ=255)
                                                                                                             gTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                      <  1:143767/251‑1 (MQ=255)
                                                                                                                    cgATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTAACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGt                                                                                                                               >  1:50686/1‑251 (MQ=255)
                                                                                                                       tAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGGTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGgttt                                                                                                                            >  2:205456/1‑251 (MQ=255)
                                                                                                                        aaCGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCACGGCAACGTGGGGATAGCAGCGCAGCACGCCGGTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCGTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           <  1:21923/251‑1 (MQ=255)
                                                                                                                                 gATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACGGGTTAGCGTGGCTATCCCCTCACCCAGGGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCCAACTCCGAGGCGCTGAAATCCAGGCTTTCACCAGCGTTCACGATGCGTTTTTTCGCggggg                                                                                                                  >  2:453059/1‑251 (MQ=255)
                                                                                                                                     gCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACt                                                                                                              >  1:148278/1‑251 (MQ=255)
                                                                                                                                      cATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTg                                                                                                             <  1:501837/251‑1 (MQ=255)
                                                                                                                                       aTCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCATCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCCGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTgg                                                                                                            >  1:350342/1‑251 (MQ=255)
                                                                                                                                         cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGc                                                                                                           <  1:252780/250‑1 (MQ=255)
                                                                                                                                         cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGc                                                                                                           <  1:144869/250‑1 (MQ=255)
                                                                                                                                         cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCg                                                                                                          <  2:155109/251‑1 (MQ=255)
                                                                                                                                         cGGTATTTTTCAGCCGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCg                                                                                                          <  2:190500/251‑1 (MQ=255)
                                                                                                                                           gTATGTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCgag                                                                                                        <  2:243236/251‑1 (MQ=255)
                                                                                                                                            tATTTTTCAGCGGGAGCGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGGGGCTATCACCTCACCCCGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTTACTGAGCCAAACTCAGAGGCGCTGAACTGCACGCTTTCACCAGCGTTCACGATGCGTTTTTGAGCGGGGGCACTGGCg                                                                                                          >  2:519523/1‑248 (MQ=255)
                                                                                                                                             aTTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGaa                                                                                                      <  2:587761/251‑1 (MQ=255)
                                                                                                                                                 ttCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCa                                                                                                   <  2:468578/250‑1 (MQ=255)
                                                                                                                                                   cAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCACAGCACGCCGTTGCGGTAAAGGGTTAGCGTGGATATCACCCCACCCAGTGGTATATCGCTGAAGCTCCGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGAGCTTAACTGCAGGCTTTCACCAGCGTTCCCGATGCGTTTTTGCGGGGGGCACAGGCGAGCAGCACTTc                                                                                               >  1:397112/1‑251 (MQ=255)
                                                                                                                                                    aGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGGAACGGGGTAGCGTGGCTATCACCTCACCCCGGGGTTTATCGCTGAAGCTCAGGAGGTCCTTCGTCCAGCTTGTAATCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAgagcttc                                                                                               >  2:187393/1‑245 (MQ=255)
                                                                                                                                                     gCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGGTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTcc                                                                                              <  1:365011/251‑1 (MQ=255)
                                                                                                                                                       gggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCCGCCAGCCGTTGCGGTAACGGGTTAGCGTGGCTAACACCTCCCCCCGTGGTTTATCGCTGAAGCTCAAGAGTTCCTTCGTCCCGCTTGTAATCTCGTCCTCCAGCGGTTTCACTGCGCCAAACTCAGAAGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGACGCACTTCcac                                                                                            >  1:568113/1‑251 (MQ=255)
                                                                                                                                                       gggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTTAAGCTCAGGATGTCCTTCGTCCAGCTTGTTCTCTCGTCATCCATCGGTTTCACTGCGCCAAACTCCGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCcac                                                                                            >  1:192444/1‑251 (MQ=255)
                                                                                                                                                        ggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGAGCAGACCGCCGTTGCGGTAAACGGTTAGCCTGGCTATCACCTCACCCTGTGGCTTATCGCTGACGCTCAGGGTGTCCTTCGTCCCGCTTGTACTCTCATCACCCAGCCGCTTCCCTACGCCAAACTCAGAGGCGCTGAACTGCAGCCTTTCCCCAGCTTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCcac                                                                                            >  1:297504/1‑250 (MQ=255)
                                                                                                                                                             gTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGGAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAACAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtg                                                                                      >  1:57242/1‑251 (MQ=255)
                                                                                                                                                             gTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGCACTGGCGAGAAGCACTTCCACAGCGtgc                                                                                     >  2:141430/1‑251 (MQ=255)
                                                                                                                                                               cccGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGCACTGGCGAGAAGCACTTCCACAGCGtgct                                                                                    <  1:141430/250‑1 (MQ=255)
                                                                                                                                                                  gCTCAGGCGCATACCTGCATCCGCGGGATCGCAGCGTAGCACGTCGGGGCGGTAACGGGTTAGCGTTGCTATTACCTCACCCATTGTTTTATCGCTGAATCATAGGATGTCCTTCGTCCAGCTTGTTCTCTCGTCAGCCATCGGTTTCACTGCGCTAAACTCATAGGCGCTGAACTGCAGGCTTGCAACAGCGTTCACTATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctt                                                                                   <  2:497394/249‑2 (MQ=255)
                                                                                                                                                                     cAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctgctgaa                                                                              <  1:434488/251‑3 (MQ=255)
                                                                                                                                                                     cAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGAACTGGCGAGAAGCACTTCCACAGCGtgctgctgaa                                                                              <  1:561414/251‑3 (MQ=255)
                                                                                                                                                                       ggCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCGCAGGATGTCCTTCGTCCAGCTTGTACGCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACt                                                                                                                                                            <  1:407205/171‑1 (MQ=255)
                                                                                                                                                                       ggCGCAGCCCCGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGCAACGGGTTAGCGTGGCTATCACCTCACCCAGGGGTTTATCGCTGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTAACTGCGCCAACCTCAGAGGGGCTGAACt                                                                                                                                                            >  2:407205/1‑171 (MQ=255)
                                                                                                                                                                          gcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTGTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCa                                                                         <  1:117553/251‑1 (MQ=255)
                                                                                                                                                                           cAGCCCGGCAACGGCGGGATCGCAGCGCCGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCCCCTCACCCCGTGGTTTATCGCTGAAGCTCAGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAAGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTcc                                                                         >  1:125989/1‑249 (MQ=255)
                                                                                                                                                                            aGCCCGGAAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCaat                                                                       <  2:1298/251‑2 (MQ=255)
                                                                                                                                                                               ccGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGc                                                                   <  2:177002/251‑1 (MQ=255)
                                                                                                                                                                                cGGCAACGGCGGGCTCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGGAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCCGGTTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGCGAAGCACTTCCACAGCGGGCTGCTGAACGTCAAGCTGc                                                                   >  2:612746/1‑251 (MQ=255)
                                                                                                                                                                                cGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGc                                                                   <  2:64617/251‑1 (MQ=255)
                                                                                                                                                                                        gcggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCGCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTGTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAAt                                                            <  2:425778/250‑1 (MQ=255)
                                                                                                                                                                                         cggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCGTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTa                                                          <  1:263131/251‑1 (MQ=255)
                                                                                                                                                                                               cgcagcgcagCACGCCGTTGCGGTAACGGGTTAGCGTTGCTATCACCTCCCCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctg                                                     <  1:340931/250‑1 (MQ=255)
                                                                                                                                                                                                 cagcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGcgcg                                                                                                                      >  2:127483/1‑183 (MQ=255)
                                                                                                                                                                                                 cagcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGcgcg                                                                                                                      <  1:127483/183‑1 (MQ=255)
                                                                                                                                                                                                   gcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGAGGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGc                                                <  2:4240/251‑1 (MQ=255)
                                                                                                                                                                                                    cgccgcaCGCCGTTGCGGTAACGGGTGAGCGTGTCTATCACCTCCCCCTGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGGCCTGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCTCTAAACTCAGATGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGGGAAATTATCCTGCTGGCg                                               <  2:424629/247‑1 (MQ=255)
                                                                                                                                                                                                       agcaCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGAc                                            <  1:407953/251‑1 (MQ=255)
                                                                                                                                                                                                             ccGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAAt                                     <  1:317963/251‑1 (MQ=255)
                                                                                                                                                                                                              cGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGCAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAAt                                     >  1:429162/1‑251 (MQ=255)
                                                                                                                                                                                                                ttGCGGTATCGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAAtt                                    <  1:132810/250‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCt                                  >  1:507736/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTg                                 >  1:492389/1‑251 (MQ=255)
                                                                                                                                                                                                                          gggTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGACCGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTaa                          >  2:83814/1‑250 (MQ=255)
                                                                                                                                                                                                                                gCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  2:312320/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTGGCTCTCACCTCAACCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCAACCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGCGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGCTGCGTTTTTGCGAGGGGGCACTGGCGAGAAGCACTTCCACCGCGTGCTGCTGAACGTCAAGCTGCGTGTACTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTACAGCAGTg                   >  2:629351/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cGTGGCTATCCCCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCATCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  2:350342/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTg              <  1:471336/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGc            >  1:139157/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTgg             <  2:106561/249‑1 (MQ=255)
                                                                                                                                                                                                                                          cacctcaccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGccc          <  1:196480/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cctcaccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGa       <  1:172550/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ctcacaCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGCGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGAc      >  2:407572/1‑251 (MQ=255)
                                                                                                                                                                                                                                               caccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGAc      <  2:475052/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516714‑4517204

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