Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,026 A→G E266E (GAA→GAG)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,0260AG94.1% 245.9 / 1.0 86E266E (GAA→GAG) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (39/41);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAC  >  NC_000913/4618781‑4619268
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
gCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAATCGGTGCAGTTCCTGACGTGTGAGTATCg                                                                                                                                                                                                                                                  <  1:141819/251‑1 (MQ=255)
 ccTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGt                                                                                                                                                                                                                                                 >  1:373146/1‑251 (MQ=255)
          tGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGGATCGTAACCCGCGt                                                                                                                                                                                                                                        >  1:133928/1‑251 (MQ=255)
            gCGAAGAAACTGGTGGAAGGTCGGGACGCGCTGGTGATGGACGGTAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCGCTGAAGCCCGTGCCGAAGCGATTTTTGGCGAGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTGCGGTAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt                                                                                                                                                                                                                                      <  2:346346/251‑1 (MQ=255)
             cGAAGAAACTGGCGGAAGGACTGGCCGAGCTGGTTATGGAAGTGAAAGTGGGTAGCGGCGCGTTTATGCCGCCCTACGCACTCTCTGCAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACAGGAATCCGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGCAGTTCGTGCAGCGGTGCAGTTCCTGACGGGTGGTTATCGTAACCCGCGtatg                                                                                                                                                                                                                                     >  2:415265/1‑248 (MQ=255)
                   aaCTGGCGTAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGGCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGAt                                                                                                                                                                                                                               <  1:426169/251‑1 (MQ=255)
                   aaCTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGAt                                                                                                                                                                                                                               >  1:416079/1‑251 (MQ=255)
                        gcggaaggtcatgcctcgcTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTCTGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCTTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTCCTCACCGACATGAAGCAGGTACTTGCCGCCAGTGCAGTTAACGCGGTTGAAGTTCGTGAAGCGGTTCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCa                                                                                                                                                                                                                           <  2:35302/234‑1 (MQ=255)
                        gCGGAAGGTCTGGACGCCCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCGCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCAc                                                                                                                                                                                                                          <  1:250672/251‑1 (MQ=255)
                          ggAAGGTCTGGACGCGCTTGGGATGGCCGTGAACGTGGGGAGCGGCGCGGTTATGCCGTCCTACGAACTCTCTGATGCTCTTGCCGAAGCGATTGTGGGAGTGGCTAACGGCTCTGGCGTGCGCCCCACCGCGCTGCTCACAGACATGAATCAGGGACGGGCCTCCAGGGTAGGTAACGCGGTTTAAGGGCGCGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGa                                                                                                                                                                                                                        <  2:361638/251‑1 (MQ=255)
                           gAAGGTCTGGCCGCGCTGGTGATGAACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTCCGAACTCTCTGAAGCACTTGCCGAAGCGAATGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATGAATCAGGAACTGGCCTCCAGGGCCGGTAACGCGGTTGAAGTTCATGAAGCGGTGCAGTTACTGTCGGGGGAGTATCGTAACCCGCGTCTGTTTGGTGACACGAt                                                                                                                                                                                                                       >  2:608535/1‑251 (MQ=255)
                            aaGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTCGCGGCGCGTTGATGCCGACCTACGAACTCTCTGAAGCCCTGGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTTCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGGCACGATg                                                                                                                                                                                                                      <  2:48980/251‑1 (MQ=255)
                              ggTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGCCGGGTGAGGATCGTAACCCGCGTCTGTTTGATGTCCAGATGgc                                                                                                                                                                                                                    >  2:542502/1‑251 (MQ=255)
                                    gACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCtgtg                                                                                                                                                                                                              <  1:155558/251‑1 (MQ=255)
                                              tGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGGGGCTAAAGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGAAATGCATCAGGTACTGGACTCCAGTGCAGGTGACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACACGCGTCTGTTTGATTTCCCGCTGGCGCTGTGCGTgggggtg                                                                                                                                                                                                    >  2:74842/1‑246 (MQ=255)
                                               gATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATg                                                                                                                                                                                                    <  1:150452/250‑1 (MQ=255)
                                                          aaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAGTAACGTAACCCGCGTCTGTTTGATTTCACGATTGCGCTGTGCGTGGAGAAGCTTATCTCCGGc                                                                                                                                                                                        >  2:366763/1‑251 (MQ=255)
                                                            gTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATAAATCAGGGACTGGCCTCCAGTGCCGGGAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGGGGATCTTAACCCGCGTCTGTTTGATGTCACGGTGGCCCTGTGCGTGGAGATTCTTATCTCCGGCa                                                                                                                                                                                       >  2:411152/1‑250 (MQ=255)
                                                                      gcgcGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGaaa                                                                                                                                                                            <  1:542502/251‑1 (MQ=255)
                                                                          gTTTATGCCGACCTGCGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATg                                                                                                                                                                        >  1:504690/1‑251 (MQ=255)
                                                                          gTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATg                                                                                                                                                                        >  1:652778/1‑251 (MQ=255)
                                                                                gCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGcc                                                                                                                                                                   <  1:26380/250‑1 (MQ=255)
                                                                                gCCGACCTACGAACTCGCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCTACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                                                                                                                                  <  1:611497/251‑1 (MQ=255)
                                                                                    aCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt                                                                                                                                                                                                                                      <  1:41193/179‑1 (MQ=255)
                                                                                    aCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt                                                                                                                                                                                                                                      >  2:41193/1‑179 (MQ=255)
                                                                                     ccTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGTCAAACTGGCGAAAGATGACGCCGAAgcg                                                                                                                                                             >  1:498393/1‑251 (MQ=255)
                                                                                     ccTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcg                                                                                                                                                             <  1:469526/251‑1 (MQ=255)
                                                                                           aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGAAATGGCCTCCAGTGCAGGTAACGCGGCTGAAGTTCGTGAAGCGGTGCAGTTCATGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTAACTCCGTCAAACTGGCGAAAGCTGACGCAGAATcgcgcgcg                                                                                                                                                       >  1:57357/1‑251 (MQ=255)
                                                                                                     aGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑Caag                                                                                                                                            <  1:45394/251‑3 (MQ=255)
                                                                                                      gCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑Caagc                                                                                                                                           <  1:37382/251‑4 (MQ=255)
                                                                                                        ccTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑Caagcag                                                                                                                                         <  1:376692/251‑6 (MQ=255)
                                                                                                           ggCCGAAGCGATTGTGGGCGTTGCTAACGGCGCTGGCGTGCGCTCCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑Caagcagtg                                                                                                                                       <  2:676389/249‑8 (MQ=255)
                                                                                                             ccGAAGCGATTGTTGCGGTGGCTACCGGCGCTGGCGTGCGCACCACCGCGCTGCCCACCGACATGAATCAGGTACTGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCTTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTG‑CAAGCAGTGCTg                                                                                                                                    <  2:362594/251‑1 (MQ=255)
                                                                                                               gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAAACCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGGGGAGATGCTTATCGCCGGCAAACTGGCGCAAGATGACGCCGAAGCGCGCGCGACATTG‑Caagcaaggctgg                                                                                                                                   >  2:109937/1‑239 (MQ=255)
                                                                                                                aaGCGATTGTGGGCGTGGCTAACGGCGCTGGCGTGCGCTCCACCGCGCTGCTCACCGAAATGAATCAGGTACTGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTgg                                                                                                                                   <  2:652778/249‑1 (MQ=255)
                                                                                                                  gCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACaa                                                                                                                               <  1:366763/251‑1 (MQ=255)
                                                                                                                   cGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAAc                                                                                                                              >  1:559928/1‑251 (MQ=255)
                                                                                                                                     aCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGATCCTGACGGGTGAGGATCGAAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGAAAACGGTaagacggccgaagac                                                                                                            >  2:196144/1‑238 (MQ=255)
                                                                                                                                        gcgcTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCttt                                                                                                         >  1:330476/1‑251 (MQ=255)
                                                                                                                                           cTGGCGTGCGCACCACCGCGCTGCTCACCGACATGACTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGTCCAACTGGCGAAAGATGACGCCGCAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCCGAAGGCTTTGGt                                                                                                      >  2:112782/1‑251 (MQ=255)
                                                                                                                                                 tGCGCACCACCGCGCTGCTCACCGACATGAATCCGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATg                                                                                                >  1:473201/1‑251 (MQ=255)
                                                                                                                                                 tGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATg                                                                                                >  1:622657/1‑251 (MQ=255)
                                                                                                                                                     cccccccgagctgtgcACCGACATTGATCGGTTACTGGCCTCCAGTGCAGGTATCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTGCTGACGGGTGAGTGTCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTTTGCGTGGTGATGCTTATCTCCGGCAAACTGGCGAAAGAGGACGCCGAAGCGCGCTCGTAATTG‑CAAGCAGTGCTGGACAACGGTATACCGGCAGAAGTCTTGGGTCGTATGGTAg                                                                                            <  2:336807/236‑1 (MQ=255)
                                                                                                                                                     caccacCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAg                                                                                            <  1:549798/251‑1 (MQ=255)
                                                                                                                                                     caccacCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAg                                                                                            <  1:251017/251‑1 (MQ=255)
                                                                                                                                                      cccaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑c                                                                                          <  1:566719/250‑1 (MQ=255)
                                                                                                                                                              gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGACGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑aa                                                                                  >  2:546532/1‑250 (MQ=255)
                                                                                                                                                              gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑aaa                                                                                 <  1:625342/251‑1 (MQ=255)
                                                                                                                                                              gctgctCACCGAAATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGCCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGAAC‑aa                                                                                  >  1:382488/1‑250 (MQ=255)
                                                                                                                                                                  ctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAgg                                                                             >  1:171774/1‑251 (MQ=255)
                                                                                                                                                                    cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGcc                                                                           <  1:588932/251‑1 (MQ=255)
                                                                                                                                                                      ccGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGccc                                                                          >  1:274604/1‑250 (MQ=255)
                                                                                                                                                                        gACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGGCGTATGGTAG‑CGGCACCAAAAAGGccg                                                                          >  1:275355/1‑248 (MQ=255)
                                                                                                                                                                        gACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGGAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCccggg                                                                       >  1:47426/1‑249 (MQ=255)
                                                                                                                                                                        gACATGAATCAGGTACTGGCCGCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCccg                                                                         <  2:275355/249‑1 (MQ=255)
                                                                                                                                                                        acaaTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCccg                                                                         <  2:47426/246‑1 (MQ=255)
                                                                                                                                                                             gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACt                                                                  >  1:480116/1‑251 (MQ=255)
                                                                                                                                                                             gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACt                                                                  <  1:109937/251‑1 (MQ=255)
                                                                                                                                                                                  aGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCATCCGACTTCGtt                                                             <  1:577704/251‑1 (MQ=255)
                                                                                                                                                                                      aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTgaga                                                         <  1:196144/251‑1 (MQ=255)
                                                                                                                                                                                      aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGAAGTTCCTGCCGGGTGAGTATCGTAACCCGCGTCTGTTTGCTGTCACGATGGCGCTGTGCGGGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGAAGCCGAAGCGCGAGAGAAATTGACAG‑CAATGCTGGACACCGGGAAAGCGGCAGAAGTCTTTGGGCGTATGG‑AGACGGCAC‑AAAAAGACCCAACCGAATTCGTTgtg                                                          >  2:406726/1‑248 (MQ=255)
                                                                                                                                                                                       cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTg                                                                                                                                    <  2:303906/177‑1 (MQ=255)
                                                                                                                                                                                       cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTg                                                                                                                                    >  1:303906/1‑177 (MQ=255)
                                                                                                                                                                                               cAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTAcgcg                                                 >  1:48413/1‑250 (MQ=255)
                                                                                                                                                                                                aGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCAGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTACGGCAAACTGGCGAAAGATGACGCCGAAGCGAGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGACAAAAGACTTTGATCGTATGGTAA‑CGGCAC‑AAAAAGGCCCAACAGACTTCGtgtaggactagtcga                                                >  2:327934/1‑236 (MQ=255)
                                                                                                                                                                                                aGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGa                                                >  1:135213/1‑250 (MQ=255)
                                                                                                                                                                                                  tGCAGGTACCGCGGTTGACGTTAGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTATTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGAATTAGTTGAGAACTACGCGAAgg                                             >  2:320613/1‑250 (MQ=255)
                                                                                                                                                                                                  tGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                             >  1:176690/1‑251 (MQ=255)
                                                                                                                                                                                                     ccgattgtcgGTTGAAGTTCTGGAAGCGGTGCAGTTCCGGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                          <  2:601410/243‑1 (MQ=255)
                                                                                                                                                                                                         aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                                                                                                                                  >  1:676633/1‑130 (MQ=255)
                                                                                                                                                                                                         aaCGCGGTTGAAGTTCGAGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                      <  1:493174/251‑1 (MQ=255)
                                                                                                                                                                                                          aCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCg                                     <  1:608535/251‑1 (MQ=255)
                                                                                                                                                                                                            gcgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                          >  1:601410/1‑244 (MQ=255)
                                                                                                                                                                                                             cgGTTGAAGTTCGTGAAGCTGGGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCGGTTTGATGTCACGATGGCTCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGTGAAATTG‑CAAGCAGTGCTGGACAACGGTAAATCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACa                                  <  2:675936/251‑1 (MQ=255)
                                                                                                                                                                                                             cgGTTGAAGTTCGTGAAGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGACTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGAAGAAGTATTTGGTCGTAAGATAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAATGATCTGCCGAcc                                  >  2:630675/1‑250 (MQ=255)
                                                                                                                                                                                                               gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGaa                                                                                                                                                                >  2:458664/1‑126 (MQ=255)
                                                                                                                                                                                                               gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGaa                                                                                                                                                                <  1:458664/126‑1 (MQ=255)
                                                                                                                                                                                                                        cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                       <  1:204290/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                       >  1:665587/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGt              <  1:660192/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATTATCGCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCa                <  2:204937/249‑1 (MQ=255)
                                                                                                                                                                                                                                     agttccgtaCGGTTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGGTGGGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTAt          <  2:504690/243‑1 (MQ=255)
                                                                                                                                                                                                                                     aGTTCCTGACGGGAGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGCAGTATCTGCCGACCGCGATGCTGACGAAAGCAGTCTAt          >  1:354986/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGAt    <  1:623405/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAc  <  1:546532/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCACCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAc  >  2:596344/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTG‑CAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAG‑CGGCAC‑AAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAC  >  NC_000913/4618781‑4619268

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