Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,391 A→G E160E (GAA→GAG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,3910AG96.5% 262.2 / ‑4.6 86E160E (GAA→GAG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (35/47);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTG  >  NC_000913/4626150‑4626634
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGc                                                                                                                                                                                                                                             >  1:32440/1‑251 (MQ=255)
     aaaGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGAt                                                                                                                                                                                                                                        <  1:262202/251‑1 (MQ=255)
      aaGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGAtt                                                                                                                                                                                                                                       <  1:133708/251‑1 (MQ=255)
       agagTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTg                                                                                                                                                                                                                                      <  1:599706/251‑1 (MQ=255)
       agagTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGGG‑CCATGCTGAGTGGCGGTAACCCTGGTGCGGGGAAAGCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCTCTCATCGCCGTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAATGCAGCATCGAGCAGATCTGCCTGATTg                                                                                                                                                                                                                                      <  2:95270/251‑1 (MQ=255)
         ggttatgccgcgtgtgaGGCGGCGGCGTTGTGTCATGAA‑TGTCCGTTCGGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCGGTTGCAGCAAACGCTGTGCAAGCTGGCCCAGCAGATGAGAACGCTGTATGTCACCGGCTAAGAGGCGCTGCAACAGGTGGCAATTCGCGCTCATCGCCTTGGTCTTCCGACTGTCAATCTCAAGAGGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGcc                                                                                                                                                                                                                                    <  2:64209/235‑1 (MQ=255)
         agTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGcc                                                                                                                                                                                                                                    >  1:265282/1‑251 (MQ=255)
          gTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCAGTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCTCGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCg                                                                                                                                                                                                                                   <  2:494635/251‑1 (MQ=255)
          gTTCGACCGCGGACTAGGCGGCGGCGTGGGGCCAGGAAGTG‑CCATTCTGATTGGCGTTATCCCGGGTGCGGAGAAATCCACGATGTTGCTGCAAACGCCGCGCAAACAGGCCCAGGAGAGGGAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCACCAGGGGGCTATGCGCGGTCACCGCCATGACCTGCCGACGGGCTAGCTCTTTATGATGTCGGAAACCAGGATCGATCAGATCTTCCTGATTGCCg                                                                                                                                                                                                                                   <  2:36402/251‑1 (MQ=255)
              gACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCgaaga                                                                                                                                                                                                                               <  1:469417/251‑1 (MQ=255)
               aCCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTAACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGCCTGACAATCTCAATATGTTGTCGGAAACCAGCCTCGAGCAGATCTGCCTGTTTGCCgcaga                                                                                                                                                                                                                               >  2:105339/1‑250 (MQ=255)
                 cgcgTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGCTGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACCATCTCACTTTGTTGTCGGAAACCAGCATCACGCAGATATGCCTGATTGCCgaagaggg                                                                                                                                                                                                                            >  2:498873/1‑248 (MQ=255)
                  gcgTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAagag                                                                                                                                                                                                                           <  1:613449/251‑1 (MQ=255)
                    gTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGACGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCGGCCGACTGACAATCTCAAAATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTTATTGCCGAAGAAGAGc                                                                                                                                                                                                                          >  2:196819/1‑250 (MQ=255)
                      aCTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAAc                                                                                                                                                                                                                       <  1:345861/251‑1 (MQ=255)
                         aGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGTCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTTCCAACGCTGTGCAAACTGGCCCAGCAGATGAAACCGCTGTCTGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCACGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGCTTGCCGAAGAAGAGCAAccc                                                                                                                                                                                                                     >  1:265528/1‑249 (MQ=255)
                                   gtggtgCCAGGAAGTG‑CCATTCTGATTGGCGGTAAACCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACAGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGAATGACAATCTAAATATGTTGTCGGAAACCAGCATAGCACAGATCTGCCTGATTGCCGCAGAAGAGCAACCGAAGCTGATGGt                                                                                                                                                                                                          >  1:392205/1‑251 (MQ=255)
                                        gCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGAATGCCGAAGAAGCGCAACCGAAGCTGATGGTAATTg                                                                                                                                                                                                     >  1:388471/1‑251 (MQ=255)
                                          cAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAc                                                                                                                                                                                                   <  1:350181/251‑1 (MQ=255)
                                            ggAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTc                                                                                                                                                                                                 <  1:48071/251‑1 (MQ=255)
                                               aGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGAt                                                                                                                                                                                              <  1:399401/251‑1 (MQ=255)
                                               aGTG‑CCATTCTGATTGGCGGGAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGAt                                                                                                                                                                                              <  2:602383/251‑1 (MQ=255)
                                                gTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATc                                                                                                                                                                                             >  1:371001/1‑251 (MQ=255)
                                                       ttCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTg                                                                                                                                                                                       >  1:342073/1‑251 (MQ=255)
                                                                 cGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATgg                                                                                                                                                                             >  1:310367/1‑251 (MQ=255)
                                                                        cctgggtcGGGGAAATCCACGTTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGATCATCGCCTGGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGATGATGGTAATTGACGCGATCCAGGTGTTGCATATGGCGGATGt                                                                                                                                                                      <  2:673215/244‑1 (MQ=255)
                                                                        ccTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGt                                                                                                                                                                      <  1:1596/251‑1 (MQ=255)
                                                                         cTGGTGCGGGGAAACCCACGCTGTTGCTGCGAACGCTGTGCAAACTGGCCCAGCAGATGAAAACACTGTATGTCACCGGCGAAGAGACGCAGCAACAGGTGGCAATGCGCGCTCAACGCCTTGGCCTGCCGGCTGACAATCTCAATATGTTGTCGGAAACCAGCATAGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCCGATGGTAATTGACTCGATACAGGTGATGCATATGGCGGgtgt                                                                                                                                                                      >  2:257533/1‑250 (MQ=255)
                                                                               cGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTc                                                                                                                                                                >  1:241219/1‑250 (MQ=255)
                                                                                ggggAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGTGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGCAGCTGATGGTAATTGACTCGATCCAGGTAATGCATATGGCGGATGTACCGTCAt                                                                                                                                                              >  2:492655/1‑251 (MQ=255)
                                                                                  ggAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCg                                                                                                                                                            <  2:265282/251‑1 (MQ=255)
                                                                                      aTCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTg                                                                                                                                                        >  1:275106/1‑251 (MQ=255)
                                                                                       tCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCTAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                                                                                                                                                       >  1:654811/1‑251 (MQ=255)
                                                                                         cACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGGTGTAccatcaacgactggcc                                                                                                                                                     >  2:505656/1‑236 (MQ=255)
                                                                                            gctgttgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCACATGGCGGAGGTAAAGTCATTGCCTGGCAGCg                                                                                                                                                  >  2:119060/8‑251 (MQ=255)
                                                                                              tgttgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTg                                                                                                                                                >  1:145562/6‑251 (MQ=255)
                                                                                               gttgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCACCTGGCAGCGTgg                                                                                                                                               >  1:495644/5‑251 (MQ=255)
                                                                                                ttgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAATAAGAGCAACCGAAGCTGATGGTAATTGACTCGAGCCAGGTGAt                                                                                                                                                                                     <  2:86926/209‑1 (MQ=255)
                                                                                                 tgcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACCGTCATCGCCTGGCAGCGTGgcg                                                                                                                                             >  1:376233/3‑251 (MQ=255)
                                                                                                  gcTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                             <  2:203936/249‑1 (MQ=255)
                                                                                                  gcTGCAAACGCGGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                             <  2:342073/249‑1 (MQ=255)
                                                                                                      cAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGGTGGGAATTGACTCGATCCAGGTGATGCATATGGCGGAGGTACAGTCATCGCCTGGCAGCgggggggcggg                                                                                                                                        >  2:557203/1‑241 (MQ=255)
                                                                                                      cAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGt                                                                                                                                        <  1:145257/251‑1 (MQ=255)
                                                                                                              gtgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGGGCGTGaa                                                                                                                                 >  1:407205/1‑250 (MQ=255)
                                                                                                                 cAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAAcgg                                                                                                                             <  1:637193/251‑1 (MQ=255)
                                                                                                                  aaaCGGGCCCAGCAGAGGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAAcggc                                                                                                                            <  2:145562/251‑1 (MQ=255)
                                                                                                                       ggCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCtt                                                                                                                       <  1:615940/251‑1 (MQ=255)
                                                                                                                          ccAGCCGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAtc                                                                                                                    <  2:349029/251‑2 (MQ=255)
                                                                                                                             gcaaaTGAAAACGCTGTATGTCCCCGGCGAAGAGTCGCTGCAACAGGTGGCAAGGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGa                                                                                                                 <  2:275106/247‑1 (MQ=255)
                                                                                                                                gATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGacgc                                                                                                              >  1:406871/1‑249 (MQ=255)
                                                                                                                                              attccACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAAcg                                                                                               <  2:253632/247‑1 (MQ=255)
                                                                                                                                              aTGTCACCGTCGACGAGTCGCTGCAACAGGTAGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGAATGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGCCTCGATCCAGGTGATGCATATGGCGGCTGTAACGTCATCGCCTGGACGCGTGGCGCAGGTGCGTGAAACGGGGGCTTATCTGAC‑GGCGTTCGCCAAAAcc                                                                                               >  2:431334/1‑250 (MQ=255)
                                                                                                                                              aTGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                             <  2:393171/206‑1 (MQ=255)
                                                                                                                                              aTGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGgcg                                                                                                                                             >  1:393171/1‑206 (MQ=255)
                                                                                                                                              aTGTCACCGGAGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGCCAATCTCAATATGTTGTCGGAAACCAGCCTCGAACAGCTATGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATAGCGGATGGACAATCTACACACGAGAGCGTGGCGCAGGTGCGTGAAACAGCGGGTTATTTGAACCGC‑TTTGCCaacaac                                                                                               >  2:607631/1‑247 (MQ=255)
                                                                                                                                                 tCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATGGCGGAAGTACAGGCATCGCCAGGCAGCATGGCGCCGGGGCGTGAAACGGCGGCTTATCTGAAGCGC‑TTAGCCCAAAagcgc                                                                                            >  2:560386/1‑251 (MQ=255)
                                                                                                                                                  cACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAAcgcgc                                                                                            <  2:104884/250‑1 (MQ=255)
                                                                                                                                                     cGGCGACGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGgtg                                                                                        <  2:149632/251‑1 (MQ=255)
                                                                                                                                                     cGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGACTGCCGACTGACAATCTCAATATGTTGTCTGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGgtg                                                                                        <  1:396613/251‑1 (MQ=255)
                                                                                                                                                        cGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTgg                                                                                     <  1:147855/251‑1 (MQ=255)
                                                                                                                                                          aaGAGTCGCTGCAACAGGTTGCTATTCGCGCTCATCGCTTTGGCCTGCCGACTGACAATCTCAAGATGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGGCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCg                                                                                   <  2:406088/251‑1 (MQ=255)
                                                                                                                                                           agagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGa                                                                                  <  1:565282/251‑1 (MQ=255)
                                                                                                                                                                  cTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCCGGGGCGTGAAACGGCGGCTTATCTGCCGCGC‑TTCGCCAAAACGCGCGGTGTGGCGGTTGTCAt                                                                           >  1:418580/1‑251 (MQ=255)
                                                                                                                                                                   tGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGGAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCGTCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGGCGCGC‑TTCGCCACAACTCGCGGTGTGGCGATTGGCAtg                                                                          >  1:117198/1‑251 (MQ=255)
                                                                                                                                                                      aaCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggtt                                                                       <  2:663535/251‑2 (MQ=255)
                                                                                                                                                                       aCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttg                                                                      <  1:196819/251‑3 (MQ=255)
                                                                                                                                                                          gggggCAATGCGCGCTAATCGCCGTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACGCGATCCAGGTGATGCATTTGGCGGATGTACAGTCATCGACTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAGCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGTCGATTGTCAtgtttggac                                                                   <  2:357696/248‑8 (MQ=255)
                                                                                                                                                                          ggTGGCAAAGCGCGCTCAACTCCTTGGCCTGCCGACGGCAAATCGCAATAGGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                   <  2:184053/251‑6 (MQ=255)
                                                                                                                                                                           gTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                <  1:565705/206‑1 (MQ=255)
                                                                                                                                                                              gCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGt                                                               <  1:557203/251‑1 (MQ=255)
                                                                                                                                                                                aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTaa                                                             <  1:363720/251‑1 (MQ=255)
                                                                                                                                                                                aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTa                                                              >  1:198512/1‑250 (MQ=255)
                                                                                                                                                                                     gcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGAGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGGGGCGGTTGTCAtggttggaaccgtaaccaa                                                         >  1:351375/1‑235 (MQ=255)
                                                                                                                                                                                           aTCGCCTTGGCCTGCAGACTGACAATCTTAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCTTGATTGCCGAAGAAGAGCAACTGAAGCTGATGGTAATTGACTCGATCCAGGTGATTCATATGGCTGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggc                                                  <  1:606464/251‑2 (MQ=255)
                                                                                                                                                                                            tCGCCATGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATAAGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGGTTATATGCCGCGC‑ATCGCCAAAACGCCCGGGGCGGGGATTGTCAtgggtgggaacgaaacaaaaaaagggc                                                 >  2:580576/1‑227 (MQ=255)
                                                                                                                                                                                                    gCCTTCCGACTGATAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCGGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATTTACAGTCATCGCCTGGCAGCGTGGCGGAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAATATGGCTCgctggca                                         <  2:308959/251‑2 (MQ=255)
                                                                                                                                                                                                     ccTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCgctggcgg                                        <  1:360040/251‑3 (MQ=255)
                                                                                                                                                                                                       tgcagACTAACAAGCTCAATATGTTGTCGGATACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGACGAGCAAACGAAGATGATGGAAATTGACTCGATCCAGGTGATGCATATGGGGGATGTACAGTCATCGCCTGGCAGCGTGGCGCTGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGATGGGGCGATTGTCCTAGTTGGACCCGTAACCAACGATGGGtggctggcgggg                                      >  2:355585/5‑247 (MQ=255)
                                                                                                                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGccc                                    >  1:558923/1‑251 (MQ=255)
                                                                                                                                                                                                               gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAAACAAAGATGGCTCGCTGGCGGGCCCTAAAGt                              >  1:476342/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                          <  1:492655/247‑1 (MQ=255)
                                                                                                                                                                                                                           atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACt                   <  1:431334/250‑1 (MQ=255)
                                                                                                                                                                                                                           atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTg                  >  1:578348/1‑251 (MQ=255)
                                                                                                                                                                                                                             gttgtCTGAAACCAGCATCGAGCAGATCTGCCGGATGGCCGAAGAAGAGCATCCGAAGCGGATGTGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGGTATCGCCTGGCATCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGATGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                <  2:328378/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggAAACCAGCATCGAGCAGATATGCATGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGGGGCGCAAGGGCGTGAAACAGCGGCATATCTGACGCGC‑TTCGCCAAAACGCGCGGGGGGGCGATTGTCAAGGGTGGGCACGTAAACAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGAc           >  1:365655/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg  <  1:105339/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTG‑CCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGC‑TTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTG  >  NC_000913/4626150‑4626634

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: ‑