Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,107 C→T G207D (GGT→GAT)  insH5 ← IS5 family transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,1070CT98.9% 260.3 / ‑4.0 92G207D (GGT→GAT) insH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (48/43);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTG  >  NC_000913/1427860‑1428352
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
cACGAAGCCGAACTGCCGCGTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGGTTTGTTCGTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCa                                                                                                                                                                                                                                                    <  1:338744/251‑1 (MQ=255)
   gAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGTTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCt                                                                                                                                                                                                                                                 <  1:372155/251‑1 (MQ=255)
    aaGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAg                                                                                                                                                                                                                                                   >  1:407886/1‑248 (MQ=255)
    aaGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAg                                                                                                                                                                                                                                                   <  2:407886/248‑1 (MQ=255)
    aaGCCGAACTGCCGCTTGATGATGCGAAACGGGCGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTg                                                                                                                                                                                                                                                >  2:382771/1‑251 (MQ=255)
      gCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGTATTCTGCGTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGAt                                                                                                                                                                                                                                              <  1:571051/251‑1 (MQ=255)
         gAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTg                                                                                                                                                                                                                                            >  2:194924/1‑250 (MQ=255)
                    tgatgaTGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGCCAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCg                                                                                                                                                                                                                                >  2:274165/1‑251 (MQ=255)
                     gatgatGCGAAACGGGTGCTCCACCCTGGCACGGCTGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCCAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCAACTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGTTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGATGCAGATTATCCAGCTGATTGAGGTCATGCCCGt                                                                                                                                                                                                                               >  2:384856/1‑251 (MQ=14)
                         atgCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGc                                                                                                                                                                                                                           >  1:286807/1‑251 (MQ=255)
                           gctaaaCGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTGGCGCGGATTCTGCTGCAAGGTTTTTCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCAGGCTCGTTGGcc                                                                                                                                                                                                                          <  2:337061/247‑1 (MQ=255)
                            cGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCcgc                                                                                                                                                                                                                        >  1:439588/1‑251 (MQ=255)
                            cGAAACGGGTGCGCCACCCTGGCCCGGATGCTGGCTTTCATGTATTCGAGGTTGATGGCCGTTTTGTTCTGGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCTATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCGCCTTGGTAGCCGGCATCGGCTGAGTCAAATTGCTCCTCTCCATGAAGCAGAGTATCCAGCTGATTGAGTGCGGGCTCGTTGGCcgc                                                                                                                                                                                                                        <  2:264448/251‑1 (MQ=14)
                             gAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCcgcg                                                                                                                                                                                                                       >  2:291598/1‑251 (MQ=255)
                                aCGGTTGCTCCACCCTGGCACGGATGCTGGCTTTCAGTTATTCGATGTTGATGGCCGTTTTTTTCGTGCGCGGATTCTGCGTCAAGTTGTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCTGCTGATTTAGGTCATGCTCGTTGGCCGCggtg                                                                                                                                                                                                                    <  2:361115/251‑1 (MQ=14)
                                   ggggCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGgtggtg                                                                                                                                                                                                                 >  2:542287/4‑251 (MQ=255)
                                         ccACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCa                                                                                                                                                                                                             <  2:222668/249‑1 (MQ=255)
                                                          cTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTGCCTTGCCGGGACGCTCGCCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTAAGGCCAct                                                                                                                                                                                          >  2:541621/1‑251 (MQ=255)
                                                            ggCTTTCATGTATTCGATGTTGATGGCCTTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCGCCTCTCCATGAATCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCAct                                                                                                                                                                                          <  1:542827/249‑1 (MQ=255)
                                                                ttCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTgg                                                                                                                                                                                     <  1:98253/250‑1 (MQ=255)
                                                                     gTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGGCAGGCCACTCTTGGCATCg                                                                                                                                                                                >  2:407759/1‑250 (MQ=255)
                                                                       aTTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGac                                                                                                                                                                              >  1:385715/1‑250 (MQ=255)
                                                                        ttCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCTGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGaca                                                                                                                                                                             <  2:423505/250‑1 (MQ=255)
                                                                            aTGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAAt                                                                                                                                                                        >  2:102072/1‑251 (MQ=255)
                                                                                ggATGGCCGTTTTGTGCTTGCGCTGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTTCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTgg                                                                                                                                                                    <  2:315063/250‑1 (MQ=255)
                                                                                 gATGGACGGTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGTGACGCTCGGCGATCAGCCAGTCCACATCCCCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGCCAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTTAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGTCACCAATGTggg                                                                                                                                                                   >  2:128930/1‑251 (MQ=14)
                                                                                    ggCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGcc                                                                                                                                                                 <  2:123206/250‑1 (MQ=255)
                                                                                    ggCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCt                                                                                                                                                                >  1:318316/1‑251 (MQ=255)
                                                                                        gTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCAt                                                                                                                                                            >  1:298368/1‑251 (MQ=255)
                                                                                        gTTTTGTTCTTGCGCGGATTCTGATTCAAGGTTTCGCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCCCGCGCTGTGGCGCTCCTTGGTAGCCTGCATCGGCTGAGACAAATTGCTCCTCGCCATGAAGCAGATTATCCAGCTGATTGAGATCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCTCTCTTGGCATCGACACCAATGTGGGCCTTCAt                                                                                                                                                            <  1:541621/251‑1 (MQ=14)
                                                                                          tttGTTCTTGCGCGGATTCTGCTTCAAGTTTTTTCCCTTGCCGGGACGCGCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATg                                                                                                                                                           <  2:128188/250‑1 (MQ=255)
                                                                                               tCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCaa                                                                                                                                                       >  1:72043/1‑249 (MQ=255)
                                                                                               tCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGGGACCAGGCTGGGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCaacg                                                                                                                                                     >  2:315782/1‑249 (MQ=255)
                                                                                                  tGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGt                                                                                                                                                    <  1:103326/249‑1 (MQ=255)
                                                                                                      cgGATTCTGCTTCAAGGTTTTTCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCCCCTCGGCCAGCTCCTCGCGCTGTGGCTCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACt                                                                                                                                              <  2:438317/251‑1 (MQ=255)
                                                                                                      cgGATTCTGCTTCAAGGTTTTTACTTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGCGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCCATGTGGGCCTTCATGCCAAAGTGCCACt                                                                                                                                              >  2:410590/1‑251 (MQ=255)
                                                                                                      cgGATTCTGCTTCAAGGTTTTTACCTGGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCGTCATGCCAAAGTGCCACt                                                                                                                                              <  2:560275/251‑1 (MQ=255)
                                                                                                         aTTCTACTTCAAGGGTTTGACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGAGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTAGTTGGCCGCGGTGGTGACCAGGCTGAGGGTCAGGCCACTCTTGGCCTCGACACCAATGTGGGCCTTCATGCAAAAGTGCCACtaat                                                                                                                                           >  2:311774/1‑248 (MQ=11)
                                                                                                                    aGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTg                                                                                                                                <  1:276982/251‑1 (MQ=37)
                                                                                                                      gTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGACAGCTCCTCGCGCTGTGGCGCTACTTGGTAGCCGGCCTCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGt                                                                                                                              >  1:465061/1‑251 (MQ=255)
                                                                                                                        ttttACCTTGCTGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCTGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGAAGCCGGCATCGGCTGCGCCAAATTGCTCCTCTCCATGAAGCAGATTATCCGGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCGTGGCATCGACACCAATGGGGGCCTTCATGCCATATGGCCACTGATTGCCTTTCTTGGTCt                                                                                                                            >  2:1123/1‑251 (MQ=11)
                                                                                                                         tttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCt                                                                                                                            <  1:542287/250‑1 (MQ=25)
                                                                                                                          ttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTg                                                                                                                           <  1:169723/250‑1 (MQ=25)
                                                                                                                            aCCTTGCAGGGACGCTCGGCGATCAGCCAGTACACATCCACCTCGGCCAGATCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATg                                                                                                                        >  2:80803/1‑251 (MQ=255)
                                                                                                                             ccTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGc                                                                                                                       >  1:442596/1‑251 (MQ=21)
                                                                                                                              cTTGCCAGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGc                                                                                                                       >  1:459948/1‑250 (MQ=25)
                                                                                                                               ttGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAt                                                                                                                     >  1:122054/1‑251 (MQ=21)
                                                                                                                                     ggACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCgg                                                                                                               >  1:60932/1‑251 (MQ=21)
                                                                                                                                       aCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGAt                                                                                                             >  1:156142/1‑251 (MQ=18)
                                                                                                                                             ggCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGtt                                                                                                       >  1:158514/1‑251 (MQ=18)
                                                                                                                                             ggCGATCAACCCGTCCACATCAACCTCGGCCAGCTCCTAGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAACTTGCTCCTCTCCATGAAGCAGATTATACAGATGATTGAGGTCATGCTCGTTGGCCGCGGAGGCGACCAGGCTGTGGGTCAGGCAGCTCTTGGCATCGACACCACTGTGGGCCTTCATGACAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGtt                                                                                                       >  2:304854/1‑251 (MQ=11)
                                                                                                                                               cGATCAGCCAGTCCTCATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTtgc                                                                                                     <  2:439588/251‑1 (MQ=18)
                                                                                                                                                gATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTtgct                                                                                                    >  1:67831/1‑251 (MQ=18)
                                                                                                                                                    aGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGctc                                                                                                >  2:55869/1‑251 (MQ=18)
                                                                                                                                                       cAGTCCAAATCCACCTCGGACAGCTCCTCGCGCTGTGCCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCCGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGGGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCttt                                                                                             >  2:529024/1‑251 (MQ=255)
                                                                                                                                                           ccACATCCACCTCGTCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTc                                                                                         <  1:15660/251‑1 (MQ=18)
                                                                                                                                                                 cccccTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCGCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTGCTTGGTa                                                                                   <  2:69088/248‑1 (MQ=16)
                                                                                                                                                                  cACCTCAGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATGGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCGTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTag                                                                                  <  1:304854/251‑1 (MQ=255)
                                                                                                                                                                   aCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACACATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGCTGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGt                                                                                    >  2:469107/1‑248 (MQ=18)
                                                                                                                                                                     cTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGc                                                                               >  2:92909/1‑251 (MQ=21)
                                                                                                                                                                      tCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGAACAGGCTGGGGGTCAGGCCCCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCACAGTGCCACTGATTGCCTTTCTTGGTCTGCTGCATCTCCGGATCGCGTTGCCGCTCTTTGTTCTTGGTAGAGCt                                                                              >  2:418985/1‑251 (MQ=11)
                                                                                                                                                                       cGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                                             <  1:92909/251‑1 (MQ=21)
                                                                                                                                                                        ggCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTgg                                                                            >  2:531426/1‑251 (MQ=21)
                                                                                                                                                                        ggCCAGATCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATAGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTAATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTgg                                                                            >  2:524742/1‑251 (MQ=255)
                                                                                                                                                                         gCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTATAGCTggg                                                                           <  2:330856/251‑1 (MQ=21)
                                                                                                                                                                          ccAGCGCCTCGCGCTGTGGCGCTCCTGGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATGATCCAGCTGATTGAGGTCATGCTCGTTGGACGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGt                                                                          <  1:583299/251‑1 (MQ=255)
                                                                                                                                                                                   cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatg                                                                 >  1:368474/1‑251 (MQ=21)
                                                                                                                                                                                   cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCGGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAat                                                                  <  1:384856/250‑1 (MQ=25)
                                                                                                                                                                                   cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCGGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTTAGGTCATGCGCGTTGGCCGCGGTGGTTGCCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCGTTCTTGGTCTGATGCATCTCCGGATCGCGGTGCTGCTCTTTGTTCTTGGTAGATCTGGGTGCCTCaa                                                                   <  2:527839/249‑1 (MQ=11)
                                                                                                                                                                                    gcgcTGTGGCGCTCCTGGGTAGCCGTCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatga                                                                <  2:385715/251‑1 (MQ=21)
                                                                                                                                                                                        gggggCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggt                                                            <  2:165420/248‑1 (MQ=21)
                                                                                                                                                                                         gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGGTGCCTCAATGAtgg                                                             <  1:242922/249‑1 (MQ=14)
                                                                                                                                                                                         gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtg                                                           <  1:481223/251‑1 (MQ=21)
                                                                                                                                                                                              gcTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATc                                                      >  2:62263/1‑251 (MQ=21)
                                                                                                                                                                                              gcTCCTTGGTAGACGGCATCGGCTGAGACAACTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATc                                                      >  2:256865/1‑251 (MQ=21)
                                                                                                                                                                                                    tGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATccacc                                                  <  2:317573/249‑1 (MQ=25)
                                                                                                                                                                                                     ggTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaaa                                               <  1:486455/251‑1 (MQ=25)
                                                                                                                                                                                                     ggTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTGGTTCTTGGTAGAGCTGGGTGCCTCAATGATGTTGGCATCCACCaaa                                               <  2:107635/251‑1 (MQ=25)
                                                                                                                                                                                                              aatcgaatGAGACACATTGGTCCTCTCAATGAAGCAGATTCTCAAGCTGATAGATGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGCTCGGGCCACTCTTGACATCGACACCCATGTGGGCCTTCATGCCAACGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTgg                                      >  2:355672/8‑251 (MQ=11)
                                                                                                                                                                                                                 cGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTGTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGt                                    <  1:382771/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAACGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTTGCATCCACCACAGTGCCTTGGGtc                                   >  1:285162/1‑251 (MQ=255)
                                                                                                                                                                                                                        gaCAAATTGCTCCGCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGa                            <  1:196049/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCt                       <  2:442596/251‑1 (MQ=38)
                                                                                                                                                                                                                                 ctcctcTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCCAAGTGCCTTGGGTCATCATGACGCCTGCtt                   >  2:516861/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cctcTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGAGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCg                 >  1:210513/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cccATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTACCTTGGGTCATCATGACGCCTGCTTCGgcca             <  2:122054/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aTGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgccagcc          >  2:177533/1‑251 (MQ=38)
                                                                                                                                                                                                                                            gaatgaGATTATCCAGCTGAGTGAGGTCATTCTCGTTGGCCGCGGTGGTGGCCAGGCTGTGGGTCGGGCCACTCTTGGCATCGACACCAATGTGGGACTTCATGCCAAAGTGCCACTGATTGGCTTTCTTGGTCTGATGCATCGCCGGATCGCGTTGGTGGTCTTTGTGCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCATAGTGGCTTGGGTCAACGTGACGCCTGGGTCGtccagcca         <  2:56810/245‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCcagccagc       >  2:555200/1‑251 (MQ=37)
                                                                                                                                                                                                                                              aGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCg      >  1:441854/1‑251 (MQ=37)
                                                                                                                                                                                                                                                 aGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGAtt   <  1:62263/251‑1 (MQ=37)
                                                                                                                                                                                                                                                   aTTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTg  >  2:447696/1‑250 (MQ=37)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTG  >  NC_000913/1427860‑1428352

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