Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,359 C→T I377I (ATC→ATT)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,3590CT97.9% 294.1 / NA 93I377I (ATC→ATT) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (44/47);  minor base G (0/2);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAACTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619115‑4619605
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
cgcgcgcgAAATTGCAAGCAGTGCTGGACAACGTTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCGGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTAc                                                                                                                                                                                                                                                        <  2:3506/251‑1 (MQ=255)
  cgcgcgGAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGAAGACGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGATGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTAc                                                                                                                                                                                                                                                        >  2:581872/1‑249 (MQ=255)
   gcgcgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGc                                                                                                                                                                                                                                                     <  2:133835/251‑1 (MQ=255)
       gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACAGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                                                                                                                                                                                                 >  2:410557/1‑251 (MQ=255)
       gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGAATTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                                                                                                                                                                                                 >  1:583067/1‑251 (MQ=255)
        aaaTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCgg                                                                                                                                                                                                                                                >  1:428612/1‑251 (MQ=255)
            tgcaagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCtt                                                                                                                                                                                                                                             <  1:67854/242‑1 (MQ=255)
             gcaagcagTGCGGGACAACGGGAAAGCGGCAGAAGTCGTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCttt                                                                                                                                                                                                                                            <  2:190454/243‑1 (MQ=255)
                agcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTg                                                                                                                                                                                                                                        <  2:544252/247‑1 (MQ=255)
                 gcagTGCTGGACAACGGTAAAGCGGCAGAATTCTTTGGTCGTATGGTAGCGGAACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGCTGCTGACGAAAGCAGTCTAAGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGAAGGGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGCTTACAGCGTCGGCTTTACTg                                                                                                                                                                                                                                        >  2:465965/4‑250 (MQ=255)
                 gcagTGCTGGACAACGGTAAAGCGGCAGAAGACTTTGGTCGTATGCTAGCGGCACAAAAAGGCCCAACCGACTACGTTGCGAACTACGCGACGAATCTGCCGCCAGCGATGCTGACGAAAGCAGTCTATGACGATACCGAAGGTTTTGTCAGTGAAATGGATACCAGCGCGCTGGGGATGGCAGTGGTTACAATGGGCGGCGGTCCTCGTCAGGCATATGACACCATTGATAACAGCGTCGGCTTTACTg                                                                                                                                                                                                                                        >  1:253412/4‑250 (MQ=255)
                   agTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGata                                                                                                                                                                                                                                     <  1:504908/250‑1 (MQ=255)
                      gCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATgg                                                                                                                                                                                                                                  <  2:428612/251‑1 (MQ=255)
                               aggTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGCCTTAGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGACAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGTCAGTGGTTGCAATGGGCGGCGGTCGGCGTCAGGCATCTTACAACATTGAGTACAGCGCCGGCTTTACTGATTTGGcgcggtgtg                                                                                                                                                                                                                         >  2:135173/2‑246 (MQ=255)
                                     aGCGGCAGAAGGCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGttttgccgacc                                                                                                                                                                                                                   <  2:46069/251‑11 (MQ=255)
                                              aGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACg                                                                                                                                                                                                          >  1:262465/1‑251 (MQ=255)
                                               gtcgttGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGGTGCGATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTGGGGCGACCAGGGAGACgg                                                                                                                                                                                                         <  2:257609/247‑1 (MQ=255)
                                               gTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACgg                                                                                                                                                                                                         >  2:391037/1‑251 (MQ=255)
                                                 cTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTc                                                                                                                                                                                                       >  2:265315/1‑251 (MQ=255)
                                                          tATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGc                                                                                                                                                                                              >  1:88941/1‑251 (MQ=255)
                                                           aTGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGACGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCt                                                                                                                                                                                             >  1:593931/1‑251 (MQ=255)
                                                            ggggagcggccaaaaaaGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGGTTTGGCAGTGAAATGGATACCCGCGCGCGGGGGATGGCAGGGGTTGCAATGGGCGGCGGTCGTCGTAAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCtt                                                                                                                                                                                            <  2:80769/239‑2 (MQ=255)
                                                               tAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg                                                                                                                                                                                         >  1:212345/1‑251 (MQ=255)
                                                                    gCACAAAAAGGCCCAACCGACTTCGTTGATATCTACGCGAAGTATCTGCCGACAGCGATGCGGACGAAAGCAGTCTATGCTGATGCCGAAGGTTTTGTCAGTGAAATTGATACCCGCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGGCGACAGGCATCTGACACCATTGATGACAGCGTCGGCTTTCCTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTag                                                                                                                                                                                    <  2:354407/251‑2 (MQ=255)
                                                                      acaAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATcc                                                                                                                                                                                  >  1:26591/1‑251 (MQ=255)
                                                                      acaAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTAGCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCTTCCGCTGGCGGTTATcc                                                                                                                                                                                  <  2:532582/251‑1 (MQ=255)
                                                                         aaaaGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAcg                                                                                                                                                                               <  2:477577/251‑1 (MQ=255)
                                                                                 caaccgacTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAc                                                                                                                                                                       <  1:413603/249‑1 (MQ=255)
                                                                                           cGTTGAGATCTACGCGAAGTATCTGCCGACAGCGAGGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAaacagct                                                                                                                                                            <  2:139266/251‑4 (MQ=255)
                                                                                              tGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCCGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTGTTCCCCGCGaaaaacagaaacagccggg                                                                                                                                                          >  2:197471/1‑235 (MQ=255)
                                                                                               gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑c                                                                                                                                                        <  2:450275/250‑1 (MQ=255)
                                                                                               gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCACGCGCAAGACA‑AAAACAGCTGG‑cg                                                                                                                                                       >  2:522965/1‑250 (MQ=255)
                                                                                                  aaCTACGCGAAGTATCTGCCCACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCCGG‑CAgg                                                                                                                                                     >  1:195096/1‑250 (MQ=255)
                                                                                                    cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGAAAGACG‑AAAACAGCGGG‑CAGGa                                                                                                                                                    >  1:49819/1‑249 (MQ=255)
                                                                                                    atACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTAAGTGCAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTCGACGGTCAGCGTCCGCTGGCGGGTCTCCCAGGGAAAGGCG‑AAAACAGCTGG‑CAggg                                                                                                                                                    >  2:386883/2‑248 (MQ=255)
                                                                                                      aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAg                                                                                                                                                  <  1:222720/249‑1 (MQ=255)
                                                                                                      aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGAGTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAg‑c                                                                                                                                                <  2:363772/250‑1 (MQ=255)
                                                                                                       cgcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATTCACCCGGAAGACG‑AAAACCGCGGG‑GAGGAAg‑cgg                                                                                                                                              >  2:544205/1‑251 (MQ=255)
                                                                                                           aaGTATCTGCCGACAGCGATGCCGACGAAAGCAGTCTGTGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACGCCATGGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAG‑AAAACAGCTGG‑CAGGAAG‑CGGAGaa                                                                                                                                          <  2:173445/251‑1 (MQ=255)
                                                                                                              gatcggCCGACAGCGGTGCTGGGGAAAGCAGGCTATGGTGATACCGAAGGTGTGGTCAGTGAAGGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGGCGTCAGGCATCTGACGCCGTTGATTACAGCGGCGGCTGTACTGATATGGGGCGTTTGGGCGACCAGGGAGGCGGTCAGCGTCCGCTGGCGTTTATGCATGCGAAAGACG‑AAAACAGCTGG‑C‑GGTAGTCGGCGAAAGc                                                                                                                                       <  2:90751/246‑1 (MQ=255)
                                                                                                                   gCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTg                                                                                                                                   <  1:293573/250‑1 (MQ=255)
                                                                                                                   gCCGAAAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGt                                                                                                                                    <  2:26591/249‑1 (MQ=255)
                                                                                                                     cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGCAAGAAG‑GAAACAGCTGG‑CAGGAAA‑CGGCGAAAGCGGGGaa                                                                                                                                 >  1:183051/1‑250 (MQ=255)
                                                                                                                     cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCCCGCGGAAGACGGAAAACG‑CTGG‑CCGGAAg‑cggcgcaaacggggaa                                                                                                                                 >  2:220253/1‑239 (MQ=255)
                                                                                                                         aGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGAACAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCACGCGAAAGACG‑GAAACAGCTGG‑CAGGGAG‑CGGCGAAAGAGGGGAAAACgg                                                                                                                            >  1:382379/1‑251 (MQ=255)
                                                                                                                               gCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAAtt                                                                                                                       <  2:212345/250‑1 (MQ=255)
                                                                                                                                cTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTa                                                                                                                      >  1:503334/1‑250 (MQ=255)
                                                                                                                                   aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑Ct                                                                                                                 <  1:48800/251‑1 (MQ=255)
                                                                                                                                     gAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGTTGAAAGCGGCAATTAAA‑CTTg                                                                                                               <  1:487044/251‑1 (MQ=255)
                                                                                                                                      aaaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGGAAGCGGTGAAAACGGAAATTAAA‑CTTGc                                                                                                              >  1:613023/1‑251 (MQ=255)
                                                                                                                                         gCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCCATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGCTTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTCTCCCCGCGGCAGACG‑AAAAACGCTGG‑CAGGAAG‑CGGGGGAAGCGGTGAAAACGGCacatta                                                                                                                     >  2:209318/1‑237 (MQ=255)
                                                                                                                                         gCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTaa                                                                                                                     <  1:209318/242‑1 (MQ=255)
                                                                                                                                          cAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCCTCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAGG‑GAAAAAGCTGG‑CAGGAAG‑CGGCGGAAGCGGTGAAAACGGCAATTAAA‑ATTGCCGa                                                                                                           >  1:112672/1‑250 (MQ=255)
                                                                                                                                          cAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAGACG‑AAAACAGCCGG‑CCGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGAt                                                                                                          >  1:3793/1‑251 (MQ=255)
                                                                                                                                           aGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATa                                                                                                         <  2:88941/251‑1 (MQ=255)
                                                                                                                                            gTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATaa                                                                                                        <  1:265315/251‑1 (MQ=255)
                                                                                                                                                aTGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGGCCGCTGGCGGGTATCCACCCGGAAGACG‑AAAACGGCTGG‑GCGGAAA‑ACGCGAAAGCGGAGGAAGGGGCCATTAAA‑Ctgtccgcacaata                                                                                                     >  2:147101/1‑238 (MQ=255)
                                                                                                                                                  gCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAAACG‑AAAACAGCGGG‑CAGGAAG‑CGGCGGAGACGGGGGAAACGGCCATTAAA‑CTTGCCGATAAAGAAc‑c                                                                                                 >  2:538598/1‑251 (MQ=255)
                                                                                                                                                      aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGCTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGGGAAAGACG‑ACAACAGCTGG‑CCGGCAG‑CGGCGTAATCGGGGCAAGCTGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGa                                                                                              >  2:140167/1‑250 (MQ=255)
                                                                                                                                                      aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGCAAGACG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGGGAAAGGGGCCATTAAC‑CTTTCCGATAAAGCAC‑Cggg                                                                                              >  2:117941/1‑249 (MQ=255)
                                                                                                                                                       tACCGTAGGATTGGTCAGTGAAATGGATACCCGCGCGCTGAGGCTGGCAGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCAACCCGGTAGACGGTCAGCGCCCGCTGGCGGTTTTCCCCACGCAAGACG‑AAAACAGCTGG‑CCGGGAG‑CCGCGAAAGAGGGGAAAAGGGCAATTAAA‑CTTGCCGGTTAAGCAA‑CGGaa                                                                                             >  2:368570/1‑250 (MQ=255)
                                                                                                                                                         ccGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGc                                                                                          >  1:106884/1‑251 (MQ=255)
                                                                                                                                                            aaGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CCGGAAG‑CGGCGAAAGCGGTGAAAACGGCAATTAAA‑CTTGCCGGTAAAGCAC‑CGGAAAGcac                                                                                        >  1:246822/1‑250 (MQ=255)
                                                                                                                                                            aaGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTAAA‑ATTGCCGGTAAAGCAC‑CGGGAAGcac                                                                                        >  1:60900/1‑250 (MQ=255)
                                                                                                                                                                 tgtGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCTTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑aa                                                                                  <  1:391794/248‑1 (MQ=255)
                                                                                                                                                                   tGTCAGTGAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCGTCGTCAGGCATCGGACACCATTGATTACAGCGTCGGTGTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTAAGCCACGCGAAAGACA‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACCCC‑AACt                                                                                <  2:106884/250‑1 (MQ=255)
                                                                                                                                                                     tCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAGGCG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGGGAAAGCGGCCATTAAA‑CTTGCCGATAAAGCAA‑CGGAAAGAACACC‑AACTGTc                                                                             >  1:343345/1‑251 (MQ=255)
                                                                                                                                                                           aaaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGc                                                                       <  1:410557/251‑1 (MQ=255)
                                                                                                                                                                               ggATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGGGGTCGTCGTCAGGCATCTGACCCCATGGATTACAGCGTCGTCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCggt                                                                   <  2:108740/251‑3 (MQ=255)
                                                                                                                                                                                 aTACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCACGCGAAAGACG‑AAAACAGCCGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTCAA‑CTTTCCGCAAAAACAC‑CGGCAAGCAAACC‑AACTGTCTATCCCCGCATc                                                                 >  1:630224/1‑251 (MQ=255)
                                                                                                                                                                                    cccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGc                                                              <  1:42293/251‑1 (MQ=255)
                                                                                                                                                                                          cgcTGGGGAAGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGCGCAAGAAG‑GCAACAGCGGG‑CAGGAAG‑CGGGGGAAGGGGTGGAAGCGGGAATTAAAATTTGCGGACAAA‑CAC‑CCGAAAGCACCCCCAACTGTC‑ATCGCCGGAGCAGCGaaaa                                                         >  2:183349/1‑248 (MQ=255)
                                                                                                                                                                                           gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAAc                                                       <  2:593931/251‑1 (MQ=255)
                                                                                                                                                                                             tGGGGATGGCCGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGTGAAAACGGCAATTAAA‑CTTGCCGGTAACACAC‑CGGAAAGCAAACC‑AACTGTCTATCGCCGTATCAGCGAATAACgg                                                     >  1:193068/1‑251 (MQ=255)
                                                                                                                                                                                              ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGt                                                    <  1:558002/251‑1 (MQ=255)
                                                                                                                                                                                               ggaaTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGATATCCACGCGGAAGACG‑AAAACAGCTGGACAG‑AAG‑CGGCGAAAGCGGTGAAAGCGGCAAATAAA‑CCTACCGATTAAG‑ACACAGAAAGCCCACC‑AACTGTCTATCGCCGTAtacgcggaaaacgga                                                    >  2:120252/4‑236 (MQ=255)
                                                                                                                                                                                                       aGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatc                                           <  2:49819/251‑1 (MQ=255)
                                                                                                                                                                                                           gTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTAAA‑CTTGCCGGTAAAGCAC‑CGGAAAGCACACC‑AACCGTCTATCGCCGGCACAGCGAATAAAGGTATActtatctgat                                       >  1:359375/1‑251 (MQ=255)
                                                                                                                                                                                                               cAATGGGGGGCGGGCGGCGTCAGGCATCTGACACCAGTGATTACAGCGGCGGCTTTAGTGATATGGAGCGTTTGGGCGACCAGGTAGACGGTCAGCGTGCGCTGGCGGTTATGCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACCAC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                   <  1:107936/251‑1 (MQ=255)
                                                                                                                                                                                                               cAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                   <  1:107911/251‑1 (MQ=255)
                                                                                                                                                                                                                   gggCGGCGGTTGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTaaa                               <  2:351969/251‑1 (MQ=255)
                                                                                                                                                                                                                   gggCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTaa                                <  1:380991/250‑1 (MQ=255)
                                                                                                                                                                                                                    ggcggcGGTCGTCGTCAGGCAACTGACACCATTGATTACAGCGTGGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCACGCGAAAGACG‑GAAAAAACTGG‑CGGGGAG‑CGGGGAAAGCGGTGAACACGGCAATTTAA‑CCTGGCGCTAAAACAC‑CGGAAAAGACCCC‑AAATTTCTATCGCCGTATCAACGAGTAACGGGATActgatctgcacattttaa                               >  2:171003/1‑240 (MQ=255)
                                                                                                                                                                                                                      cggcggGCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTAGCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAttt                            <  2:600478/251‑1 (MQ=255)
                                                                                                                                                                                                                         gggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGaa                         <  1:581872/244‑1 (MQ=255)
                                                                                                                                                                                                                           gtcgtcGTCAGGCATCTGACACCAGTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                       <  1:332586/246‑1 (MQ=255)
                                                                                                                                                                                                                             cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGGAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑CACTGTCt                                                                            >  2:224561/4‑196 (MQ=255)
                                                                                                                                                                                                                             cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCt                                                                            <  1:224561/193‑1 (MQ=255)
                                                                                                                                                                                                                                       cATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGga           <  1:220253/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGga           <  1:117941/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACata    <  1:616974/251‑1 (MQ=255)
                                                                                                                                                                                                                                               acCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACatat   <  2:193068/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCGGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATg  <  2:215275/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
CGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAACTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATG  >  NC_000913/4619115‑4619605

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