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,9590CT94.6% 246.0 / ‑5.8 95P253P (CCG→CCAfecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (44/43);  minor base G (5/0);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 5.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516717‑4517204
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
gcGAACGAGAAGGGTAAAAAGGACAAGCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAGGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATc                                                                                                                                                                                                                                                     <  1:403821/250‑1 (MQ=255)
  gAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGGGGCTATCACCTCACCCAGTGGGTTATCGc                                                                                                                                                                                                                                                  >  1:285162/1‑251 (MQ=255)
     cGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGa                                                                                                                                                                                                                                               <  1:664282/251‑1 (MQ=255)
      gagaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGa                                                                                                                                                                                                                                               >  1:293106/1‑250 (MQ=255)
               aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGAATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGTTCGCAGAGAAGCACACCGTTTCCGTGACGGGCTGGCGAGGCTAACAACCCACCCAATGGGTTATTGCTGCAGCTCAggc                                                                                                                                                                                                                                      >  1:642241/1‑249 (MQ=255)
                 aaaGGACAATCGAAATAATAATTATTTTCCTTATAGTGGTGAAATGTTTATCCAGTACCGCGTAATAAACTGAATTTTAACGTGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGTTAACTTCCCTCTCAGGCTCAGCCCGGCAACGGCGGGATCGCATCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCTCTGAAGCTCAGGATGt                                                                                                                                                                                                                                   <  2:523342/251‑1 (MQ=255)
                   aGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAAGTTTAACGGGAAGCGTTTGCGCTATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTcc                                                                                                                                                                                                                                 <  1:248499/251‑1 (MQ=255)
                     gACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCtt                                                                                                                                                                                                                               <  1:500435/251‑1 (MQ=255)
                     gACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCtt                                                                                                                                                                                                                               <  1:205807/251‑1 (MQ=255)
                         aTCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTc                                                                                                                                                                                                                           <  1:174508/251‑1 (MQ=255)
                            gAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGATTGCGCGATAACGTTCAGGCTCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAACGGGGTAGCGTGGCTATCCCCTCCCCCAGTGGGTTATCGCCTCAGCTCTGGATGTCCTTCGTcccg                                                                                                                                                                                                                        >  1:86908/1‑249 (MQ=255)
                                    aaTTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTATCACCTCCCCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAct                                                                                                                                                                                                                >  1:447496/1‑251 (MQ=255)
                                              cTTACAGTTGTTAAATGTTTACCCAGTACCGCGTAATAGACTGAATTTTAACGCGAAGCGTTGGCGCGATAACGTTCAGGAGCGCATCGGTATTTTTCAGCGGGAACGGCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCTGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCTCCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGt                                                                                                                                                                                                                            <  2:568806/229‑1 (MQ=255)
                                              cTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGAGGGAACGGGTTAGCGTGGCTATCACCTCAACCAGTGGTTTATCGATGAAGCTCAGGATGTCCTTCGt                                                                                                                                                                                                                            >  1:568806/1‑229 (MQ=255)
                                                   aGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAGCgg                                                                                                                                                                                                >  1:668006/1‑251 (MQ=255)
                                                      ggtgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCCCCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAGCgggtt                                                                                                                                                                                             >  1:642709/1‑248 (MQ=255)
                                                                  aTCCAGTACCGCTGAATTGACTGAATGTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCAGCGGTATTTTGCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGGCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCaaa                                                                                                                                                                             <  2:351799/251‑1 (MQ=255)
                                                                         aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCaga                                                                                                                                                                       <  1:623328/250‑1 (MQ=255)
                                                                             cgTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGgcgc                                                                                                                                                                  <  1:572905/251‑1 (MQ=255)
                                                                               tAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGGTTGCCGCGCAGCACGCCGTTGCGGTAACGGGTTAGGGTGGCTGTCCCCCCACCCAGTGGGTTTTCGCTGAAGCTCAGGGTGTCCTTCGTCCAGCTTGTCCTCTCGTCATCCAGCGGTTTCACTGCGCCAAACCCAGAGgcgcgg                                                                                                                                                                >  2:385606/1‑249 (MQ=255)
                                                                                aaTAGGCTTAATTTTAACGGGAAGTGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTg                                                                                                                                                                <  1:124028/250‑1 (MQ=255)
                                                                                aaTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAAAGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCCGAGGCGCTg                                                                                                                                                                >  1:351246/1‑250 (MQ=255)
                                                                                aaTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACCCCGGTGCGGTAAAGGGATAGCGGGGCCATCAACCCACCCCGGGGTTTTTCGGCGAAGCTCAGGATGTCCCTTGTCCAGCTTGTACTCTCGTCATCCACAGGTTTTCATGCGCCAAAACCCGAGGCGCTg                                                                                                                                                                >  2:660739/1‑250 (MQ=255)
                                                                                   aGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGGTTCACTGCGCCAAACTCAGAGGCGCTGAACt                                                                                                                                                            >  1:367587/1‑251 (MQ=255)
                                                                                      cTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAACGCAGCGCAGCACGCCGTTGCGGGAACGGGGTAGCGTGGCTATCACCCCAACCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACAGCGCCCAACTCAGAGGCCCTGAACttc                                                                                                                                                          >  2:284191/1‑248 (MQ=255)
                                                                                       tGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCCTCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGc                                                                                                                                                          >  1:284687/1‑249 (MQ=255)
                                                                                         aaTTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACATCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCCGCGGTTTCACTGCGCCCAACTCAGAGGCGCTGCACTGCAGGc                                                                                                                                                      >  1:369704/1‑251 (MQ=255)
                                                                                             ttAACGGGAAGCGTTGGCTCGATAACTTTCAGGATCGCATCGGTATTTTGCAGCGGTAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATTTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCttt                                                                                                                                                   <  2:658978/250‑1 (MQ=255)
                                                                                              tAACGGTAAGCGTTGGCGCGATAACGTGCAGGATAGCGTCGGTATTTTTCAGCGGGAACGTCCCGCTCATGCGCAGCCCGGCACCGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGTCTATCACCTCACCCAGTGGTTTAGCGCTGAAGCGCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTc                                                                                                                                                  <  2:588174/250‑1 (MQ=255)
                                                                                               aaCGGGAAGCGTTTGCGCCATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCa                                                                                                                                                 <  2:367587/250‑1 (MQ=255)
                                                                                                aCGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAACGGGTTAGCGGAGCGATCCCCTAACCCAGATGGTTTTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGACATCCAGCGGGTTCACTGCGCCAAACTCCGAGGCGCTGAAATGCAGGCTTTCAcc                                                                                                                                               >  1:25665/1‑251 (MQ=255)
                                                                                                 cGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                                                                                                               >  1:162396/1‑250 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATGGCAACGCAGCAAGCCGTTGCGGAAACGGGTTAACGTGGCTATCAACTCACCCAGTGGTTTATCGGTGACACTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCAGGTTCAATGCGACAAACTCAGAGGCGCTGAACTGCAGGCTTTTACCAGc                                                                                                                                            >  1:107958/1‑251 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAACACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCCCCCCGGGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCCGAGGCGCTGAACTTCAGGCTTTCAcccgc                                                                                                                                            >  1:616647/1‑248 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCCCAGAGCCACACGCCGTTGCGGTAAACGGGTAGCGTGGCCATCACCTCACCCAGGGGGTTATCGCTGAAGCTCCGGCTGACCTTCCTCCACCTTGTACTCTCGTCATCCCGTCGGTTCACTGGCCCCAAACTAGAGGAGCTGAAcctcaggatttccccacc                                                                                                                                            >  2:662575/1‑234 (MQ=255)
                                                                                                     aaGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGt                                                                                                                                          >  1:387515/1‑251 (MQ=255)
                                                                                                      aGCGTTTGCGCGATAACGTTAAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACCCCACCCAGTGGTTTTACGCTGAAGCCAAGGAGGTCCTTCGCCCAGCTTGGACTCTCTCCATCCAGCGGTTTCACTGCCCCCAACTCAGAGGAGCGAAACAGCAGGCTTTTACCAGCGtt                                                                                                                                         >  2:467375/1‑251 (MQ=255)
                                                                                                          gttGCGCGATATCGTTCAGGATCGCATCGGGATGTTGCAGCGGGTACGTTCCTCTCAGGCGCAGCCCGGCAACGGCGGGACCGCAGGTCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTAGCACCGCACCCAGTGGGGTAGCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTGGTACTCGCTTCATCCAGCGGTTTCACTGCGCTAAACTCAGAGGCGCTGAACGGCAGGCTGTCACCAGCGGTCACg                                                                                                                                     <  2:86908/250‑1 (MQ=255)
                                                                                                            ggcgcGATAACGTGTAGGATCCCACGTGTGTTTTTCAGCGGAAACGTCCCGCGCAGGAGCAGCCCGGCAACGTCGGGATCGCAGCGCTGCACTCCGTTGCGGTAGCGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCGGTAGCGCAGGATGTCCTTCGTCTAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCTGCGTTCACGAt                                                                                                                                   <  2:636325/250‑1 (MQ=255)
                                                                                                             gcgcGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATg                                                                                                                                  <  1:514150/251‑1 (MQ=255)
                                                                                                               gcgATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGTAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGGCATCCAGCGGTTTCACTGCGCCAAACTCAAAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGCTGCg                                                                                                                                >  1:175727/1‑251 (MQ=255)
                                                                                                                cgATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCCGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGGACTGCAGGCTTTCACCAGCGTTCAAGATGCGt                                                                                                                               >  1:445614/1‑251 (MQ=255)
                                                                                                                cgATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGt                                                                                                                               <  1:170390/251‑1 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCCGCGCAGCACGCCGGTGCGGTAACGGGATAGCGGGGCCATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTAATCTCGTCATCCAGCGGTTTCACTGCGCCCAACTCAGAGGCGCTGAACTACAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           >  1:74631/1‑251 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           <  1:108412/251‑1 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           <  1:284191/251‑1 (MQ=255)
                                                                                                                       gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGGGGCCATCACCTCACCCAGTGGTTTATCGCTGAAGCTCCGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAACAGCGGTCACGATGCGGTTTTgcg                                                                                                                        >  1:132116/1‑251 (MQ=255)
                                                                                                                         tcatgatagCATCGGTATTTTTCAGCAGGAACGCCCCGCTCAGGCGCGGACCTGCAACTGCGGGAGCTCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGGACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCTCGATTCGTTTTTgcgc                                                                                                                       <  2:462605/242‑1 (MQ=255)
                                                                                                                             gATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGGTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggggg                                                                                                                  >  1:560022/1‑251 (MQ=255)
                                                                                                                                cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCa                                                                                                                <  1:511348/250‑1 (MQ=255)
                                                                                                                                cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAATGCAGAGCAGAAAGAGGCTGCGGTAACGGGGTAGGGCGGCCATCACCTCGACCAGTGGGTTATCAGCGGCGCTAAGGGTGTCCTTCGCCAAGCGTGTATTCTCGCCATCCAGCGGGTTCACTCGCACCAACACCAGCGCGCTACACCGCAGCCTTTCAACAGCGTTCACGATGCGgtttttagccggggcg                                                                                                                >  1:666805/1‑239 (MQ=255)
                                                                                                                                 gCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACt                                                                                                              <  1:162289/251‑1 (MQ=255)
                                                                                                                                 gCATCGGTATTTTTCAGCGGGAACGTCCCGATAAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCACCACCCCACTTCGGTAAAGGGGTAGCGAGGGTATCACCTCACCCCAGGGGTTATCGCCAAAACAAAGGAAGTACTTCCTCCAACTTGTAATACCTTCATCCACAGGTTTCCCAACTCCAAACTCAGAGGAACTTAACTGCAGGCCTTCACCCCCGTTCACGactaggttttgagacggcggccc                                                                                                              >  2:594667/1‑229 (MQ=255)
                                                                                                                                   agcGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                      <  2:400871/223‑1 (MQ=255)
                                                                                                                                   aTCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGGAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCCGGCTTTCACCAGCGTTCAc                                                                                                                                      >  1:400871/1‑225 (MQ=255)
                                                                                                                                     cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGc                                                                                                           <  1:88817/250‑1 (MQ=255)
                                                                                                                                     cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCg                                                                                                          <  1:362099/251‑1 (MQ=255)
                                                                                                                                         tATTTTTCAGCGGGAACGTACCGCTCAGGCGCAGCCCGGCAACGGCGGGCTCGCAGCGCAGCCCGCCGTTGCGGTAACGGGCTCGCGGGGCCATCACCTCACCCAGTGGGTTATCGCTCAAGGACAGGGTTGCCTTCCTCCCGCTTGTACTCTCGTCATCCAGCGGTGTCACCGGGCCAAACTCCGCGGAGCTGAACAGCAGGCTTTCACCAGCGTTACCCATACGTTTTTGCGCGGGGAGACTGGCgcg                                                                                                        >  2:435297/1‑248 (MQ=255)
                                                                                                                                                   cGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGTGAGAAGCACTTccc                                                                                             >  1:318910/1‑250 (MQ=255)
                                                                                                                                                   cGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTATCACCTCCCCCAGTGGGTTATAGCCGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGGTTCACTGCGCCAAACTCAGAGGCTCCGAACTGCAGGCTTTCACCAGCGTTCACCATACGTTTTTGTGCGGGGCGACTGGCGCGAAggccttcac                                                                                             >  2:439379/1‑243 (MQ=255)
                                                                                                                                                     ggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAGGCGCTGAACTGCCGGCTTTCACCAGCGTTCACGATGCGTTTTTTCGCGGGGGCACTGGCGAGAAGCACTTCcacc                                                                                           >  1:435973/1‑250 (MQ=255)
                                                                                                                                                          gTCCCGCTCAGGCGCAGCCCGGCACCGGCGGGATCGCAGCGCAGCACGTCGTTGCGTTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACGCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtg                                                                                      <  2:598189/251‑1 (MQ=255)
                                                                                                                                                            cccGCTCATGCGCAGCCGTGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTTACGGGTTAGCGTTGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTTCTCTCGTCATCCAGCGGTTTCTCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGGTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgct                                                                                    <  2:447496/251‑1 (MQ=255)
                                                                                                                                                             ccGCCCAGGCGCAGCCCGGCAACGGCGTGATCGCTGCGCGTCACTCCGTTGCGGTAACGTGTTAGCGTGGCTATCACCTCTCCCAGTGGTTTATCGCTGAAGCTCAGGAGGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCTCTTCCCCAGCGtgctg                                                                                   <  2:449960/251‑1 (MQ=255)
                                                                                                                                                              cGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGGTAGCGTGGCTATCACCTCACCCAGGGGTTTATCGCTGAAACTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCCAACTCAGAGGCGCTGAACTGCCGGCTTTCACCCGCGTTCCCGATGCGTTTTTGCGCGGGGGCCCTGGCGAGAAACACTTCCACAGCGggctgt                                                                                  >  2:172785/1‑250 (MQ=255)
                                                                                                                                                                   aGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctgctgaac                                                                             <  1:179737/251‑4 (MQ=255)
                                                                                                                                                                      cgcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTCGCGTGGCTATCACCTCACACAGTGGTTTTTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAGGCCCTGAACTGCAGGCTTTCACCAGCGTTCCCGATGCGTTTTTGCGCGGGGGCACTGGCGCGAAGGACTTCCACAGCGTGCTGCTGAACGTc                                                                          >  1:634668/1‑251 (MQ=255)
                                                                                                                                                                          gCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCCAGc                                                                      >  1:22214/1‑251 (MQ=255)
                                                                                                                                                                                cAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTg                                                                <  1:223241/251‑1 (MQ=255)
                                                                                                                                                                                cAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTg                                                                <  1:467011/251‑1 (MQ=255)
                                                                                                                                                                                    ggcggGATCGCAGTGCAGCACGCCGTTTCTTTATTGGGCTAGCGTGCCTTTCCCCTCCCCCAGTCGTTTATCGCTGAAGGTCAGGATGTCCTTCGTGCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCTAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTGCACGATGCGTTTTTTCGCGGGGGCACTGGCGAGAAGCACTGCCACTGCTTGCTGCTGAACGTCAAGCTGCGTGAAAt                                                            <  2:410395/251‑1 (MQ=255)
                                                                                                                                                                                     gcggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCTCAGCGTGCTGCTGAACGTCAAGCTGCGTGAAAtt                                                           <  1:176172/251‑1 (MQ=255)
                                                                                                                                                                                     gcggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAAtt                                                           >  1:589717/1‑251 (MQ=255)
                                                                                                                                                                                        tgctcgccccgcagcACGCCGTTGCGGCCACGCGATTGCGAAGCTATCCCCCCACCCCGTGCTTTCTCGCTCCAGATCATGATGTCCTTCGACAATCTTGTACTCTCGTCATCCAGCGGTTTCACCGCGCACAACTCCGACGCCCTGAACCGCAGGCTTTCAACCTAGTTCACGCTACGTTTTTACGCGGGGGCCCTGGCAAGAAGCACTTCCACAGCGTGATGCTGAAAGTCAAGCTGCggtctattacc                                                        >  2:341489/10‑241 (MQ=255)
                                                                                                                                                                                            cgcagcgcagCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGt                                                                                       >  1:517100/1‑216 (MQ=255)
                                                                                                                                                                                            cgcagcgcagCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctg                                                     <  1:45141/250‑1 (MQ=255)
                                                                                                                                                                                             gcagcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctgct                                                   >  1:48254/1‑251 (MQ=255)
                                                                                                                                                                                               agcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTgg                                                 >  1:679571/1‑251 (MQ=255)
                                                                                                                                                                                                gcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGc                                                <  1:42633/251‑1 (MQ=255)
                                                                                                                                                                                                    agcaCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGAc                                            <  1:583683/251‑1 (MQ=255)
                                                                                                                                                                                                       ccGCCGTTGCGGTAACGGGTTAGCGTGACTATCACCTCTACCCGGGGTTTATAGCTGAAGCTCAGGATGACCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCGGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAGGCTCTTCCCAAGCGTGCGGCTGATCGTCCAACTGCGTAAAATTCTCCTTCGGGGGGACggg                                         >  2:368197/2‑250 (MQ=255)
                                                                                                                                                                                                                             aGCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGCCGGTAAATTCTGTCCCTAAAGCAgg                    >  1:131453/1‑250 (MQ=255)
                                                                                                                                                                                                                              gCGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg                   <  1:385606/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggtggcttgcacctctccCGGTGGTTGATCGCTGAAGCTCAGGAGTTCCTTCGTCCAGCTTGTACTCTCGGCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  2:15954/242‑1 (MQ=255)
                                                                                                                                                                                                                               cGTGGCTCTCACCTCACCCAGTGGTTTCTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCGGGCGGACGGGAAATTCTGCCCCTAAAGCAGTgg                  >  2:350731/1‑250 (MQ=255)
                                                                                                                                                                                                                               cGTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCGCAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  1:18559/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAACGCAGTg                   >  1:448741/1‑249 (MQ=255)
                                                                                                                                                                                                                                gTGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg                   >  2:322847/1‑249 (MQ=255)
                                                                                                                                                                                                                                 tGGCTATCACCTCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGc                >  1:518223/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCTATCTCCTCGCCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGGCATCCTGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTTCACGTCAAGCTGCGTGAAATTATCCTGCTGGTGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCt               <  2:333333/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCACCTCACCCAGTTGTTTTACGCTGGAACTCCGGATGTCCTTTGTCCTGCTTGTACTCTCGTCAACCAGCCGTTTCACTGCGCCAAACTCAAAGGTGCTGGACTGCCGGCTTTCAACAGCCTTCCCGATGCGTTTTTGCTCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGcc           <  1:100337/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tcacctcaccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTAACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGccc          >  1:21824/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cccAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgct                                                                                    <  1:372135/169‑1 (MQ=255)
                                                                                                                                                                                                                                                cccAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTACATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt  >  2:107852/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cccAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGAGtgct                                                                                    >  2:372135/1‑169 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
TCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516717‑4517204

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