Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 645,673 G→A R463R (CGG→CGA)  unknown → ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1645,6730GA97.5% 222.8 / NA 81R463R (CGG→CGA) unknownATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (31/48);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCAACGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGCTATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCTTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAA  >  1/645430‑645913
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
aGTGCAACGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATAtt                                                                                                                                                                                                                                              >  1:395317/1‑250 (MQ=255)
 gTGCAACGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCa                                                                                                                                                                                                                                            >  1:312152/1‑251 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:773503/251‑1 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:520052/251‑1 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:197247/251‑1 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:699668/251‑1 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:64481/251‑1 (MQ=255)
     aaCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:487345/251‑1 (MQ=255)
      aCGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:366011/250‑1 (MQ=255)
         tCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCg                                                                                                                                                                                                                                    >  1:459919/1‑251 (MQ=255)
          cTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGt                                                                                                                                                                                                                                   >  1:672792/1‑251 (MQ=255)
            tCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTc                                                                                                                                                                                                                                             <  1:201896/239‑1 (MQ=255)
               cATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:307795/241‑1 (MQ=255)
               cATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:165657/241‑1 (MQ=255)
               cATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAAt                                                                                                                                                                                                                              >  1:359915/1‑251 (MQ=255)
                          tttCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGgaagaa                                                                                                                                                                                                                   >  1:272382/1‑251 (MQ=255)
                            tCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTGGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAAtt                                                                                                                                                                                                                 >  1:859258/1‑251 (MQ=255)
                                   ccTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACata                                                                                                                                                                                                                                                <  1:4977/213‑1 (MQ=255)
                                      ccGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:646742/218‑1 (MQ=255)
                                           aaTCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTtggaatg                                                                                                                                                                                                  <  1:544411/251‑1 (MQ=255)
                                            aTCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTtggaatgg                                                                                                                                                                                                 <  1:846611/251‑1 (MQ=255)
                                                gCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        >  1:668928/1‑208 (MQ=255)
                                                gCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:74919/208‑1 (MQ=255)
                                                gCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAgg                                                                                                                                                                                             >  1:709408/1‑251 (MQ=255)
                                                        ggCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCg                                                                                                                                                                                     <  1:478454/250‑1 (MQ=255)
                                                                   cGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑cccc                                                                                                                                                                         <  1:572059/250‑1 (MQ=255)
                                                                     tgtgAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCaaa                                                                                                                                                                                 <  1:75476/241‑1 (MQ=255)
                                                                     tgtgAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑cccccc                                                                                                                                                                       <  1:389411/250‑2 (MQ=255)
                                                                     tgtgAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑cccccc                                                                                                                                                                       <  1:210536/250‑2 (MQ=255)
                                                                               ggCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        >  1:869039/1‑177 (MQ=255)
                                                                                gCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:801723/176‑1 (MQ=255)
                                                                                gCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGt                                                                                                                                                           >  1:702270/1‑251 (MQ=255)
                                                                                 cATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTc                                                                                                                                                                                                                                             <  1:526082/170‑1 (MQ=255)
                                                                                      gTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTc                                                                                                                                                      >  1:344591/1‑250 (MQ=255)
                                                                                             tCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTg                                                                                                                                                                                                                >  1:591470/1‑187 (MQ=255)
                                                                                             tCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTccccc                                                                                                                                                  <  1:339053/247‑1 (MQ=255)
                                                                                                  aTTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACat                                                                                                                                                                                                                                                 >  1:437462/1‑149 (MQ=255)
                                                                                                  aTTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                           <  1:563789/249‑1 (MQ=255)
                                                                                                   ttCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgg                                                                                                                                        <  1:769064/251‑1 (MQ=255)
                                                                                                     cGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggg                                                                                                                                      <  1:742171/251‑1 (MQ=255)
                                                                                                       gAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGCAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggggg                                                                                                                                    >  1:454942/1‑250 (MQ=255)
                                                                                                                  caacaaTCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTa                                                                                                                         <  1:231214/251‑1 (MQ=255)
                                                                                                                         cGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATc                                                                                                                  <  1:181570/251‑1 (MQ=255)
                                                                                                                          gCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCa                                                                                                                 >  1:541026/1‑251 (MQ=255)
                                                                                                                            ccGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCaga                                                                                                               <  1:637738/251‑1 (MQ=255)
                                                                                                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                           <  1:527926/220‑1 (MQ=255)
                                                                                                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggg                                                                                                                                      <  1:738921/225‑1 (MQ=255)
                                                                                                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATc                                                                                                                  <  1:222402/245‑1 (MQ=255)
                                                                                                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                        <  1:888539/251‑1 (MQ=255)
                                                                                                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                        <  1:525827/251‑1 (MQ=255)
                                                                                                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                        <  1:495214/251‑1 (MQ=255)
                                                                                                                                            cGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGGAGCCTTCCCATTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGGAAATCT‑CCCCCCTCCGCCTTCGGTGCTTCCACCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGctc                                                                                               <  1:411094/251‑1 (MQ=255)
                                                                                                                                                          aCCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAgg                                                                                 <  1:218417/251‑1 (MQ=255)
                                                                                                                                                             aaTTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGAt                                                                                                                                                                                                                                        <  1:680773/99‑1 (MQ=255)
                                                                                                                                                                ttGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATc                                                                           <  1:384957/251‑1 (MQ=255)
                                                                                                                                                                   gCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGTGGGATCTg                                                                         >  1:42004/1‑250 (MQ=255)
                                                                                                                                                                     gTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGcc                                                                       >  1:147602/1‑250 (MQ=255)
                                                                                                                                                                            aTTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTtggaat                                                                                                                                                                                                   >  1:89320/1‑121 (MQ=255)
                                                                                                                                                                            aTTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGc                                                               <  1:117252/251‑1 (MQ=255)
                                                                                                                                                                            aTTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGc                                                               <  1:461430/251‑1 (MQ=255)
                                                                                                                                                                                   aaaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgggg                                                        <  1:606023/251‑1 (MQ=255)
                                                                                                                                                                                    aaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggg                                                       >  1:458839/1‑251 (MQ=255)
                                                                                                                                                                                    aaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggg                                                       <  1:276967/251‑1 (MQ=255)
                                                                                                                                                                                      tGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggggg                                                     <  1:783634/251‑1 (MQ=255)
                                                                                                                                                                                         cGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAAt                                                  <  1:233624/251‑1 (MQ=255)
                                                                                                                                                                                           gtgtTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAggggg                                                                                                                                     >  1:848044/1‑166 (MQ=255)
                                                                                                                                                                                                  gTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                           >  1:234460/1‑153 (MQ=255)
                                                                                                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAg                                                                                  >  1:32710/2‑201 (MQ=255)
                                                                                                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGaa                                                   >  1:7202/2‑232 (MQ=255)
                                                                                                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGaa                                >  1:542910/2‑251 (MQ=255)
                                                                                                                                                                                                             aTTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGaaaa                              >  1:693197/1‑251 (MQ=255)
                                                                                                                                                                                                                   gCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATAc                                                >  1:519697/1‑227 (MQ=255)
                                                                                                                                                                                                                         ccGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTcac                  <  1:341354/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACa                                                                                   >  1:498687/1‑184 (MQ=255)
                                                                                                                                                                                                                            ttcggtgATGGACCGGATAAACGGCATATTCCGGGTGGGGTGAAATTGTTGGAAGGAGTGGGGGCGTTTGAATGGAAGGCCTA‑TCGCAAATCTGCCCCCTT‑CGCCTTCCGTTAGTCCCCCCCTTAACCGGGGGAGGGAAAGGGTATTTTATCAGGAGGAGGAACTTCGCTCCCCTCGTTTTCAAGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTcacac                <  1:444494/244‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTtc                      >  1:637673/1‑244 (MQ=255)
                                                                                                                                                                                                                                tGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGaa                                                                                                                               <  1:39958/135‑1 (MQ=255)
                                                                                                                                                                                                                                     gCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTAttt                                                                                                                      >  1:223578/1‑139 (MQ=255)
                                                                                                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAGTGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAAGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCaa  >  1:445562/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGGAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgg                                                          <  1:607284/194‑1 (MQ=255)
                                                                                                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgg                                                          <  1:340012/194‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
AGTGCAACGTCTTCCCATTGAACCCATTTCACAAGCCTCCGCCAATCAGCGTAAAGGGCGTTGCGGGCGTGTGAGCGAAGGCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGCTATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCTTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAA  >  1/645430‑645913

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