Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,513,841 A→G R156R (CGT→CGC fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,513,8410AG96.0% 289.6 / 1.0 100R156R (CGT→CGCfecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (58/38);  total (60/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

CGCGATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCA  >  NC_000913/4513594‑4514084
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
cgcgATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGt                                                                                                                                                                                                                                                  >  1:670623/1‑251 (MQ=255)
 gcgATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTg                                                                                                                                                                                                                                                 >  1:12906/1‑251 (MQ=255)
      aGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGc                                                                                                                                                                                                                                            >  1:680353/1‑251 (MQ=255)
        tGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCAt                                                                                                                                                                                                                                          >  1:62340/1‑251 (MQ=255)
         gggCACCCAGCATCCGGGCAGGATGGACCGCCAGCAGTTGCTCCAGGCCGAGGGACGGCATGGACGCCCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCAtc                                                                                                                                                                                                                                         <  1:527360/251‑1 (MQ=255)
         gggCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCAtc                                                                                                                                                                                                                                         >  1:45339/1‑251 (MQ=255)
             aaCCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCAtctctc                                                                                                                                                                                                                                     >  1:611141/1‑251 (MQ=255)
                          gcaggATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGGCGCACCCGCCATCGCAGCGGGAACGTTCAGACCCAGAGAGGCCAGCACGCTGCCGGTACCGGTCTCCTGAGTATGCAGGTTCAATTGCTGTTCGCGTGATGTGCCAAAAGCCCCGCGGGTCCCATTGGGAAGCTGGCTGGCCCACGGCGCCATCCCCTCTTTATGTTGTTCCAGGGGTCCCTCCATCTCTCGCTTTTTAccccc                                                                                                                                                                                                                        >  2:369588/1‑249 (MQ=255)
                            aggATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACccaccc                                                                                                                                                                                                                      >  2:369228/1‑250 (MQ=255)
                              ggttGACCGCCAGCATTTGCGTTAGGCCTCTGGACGTAATGGACGCCCCCGCCAGCGCAGCGGGAACGTTCAGCCCCAGAGAGGTCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGGCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTGTATGTTGTTCCAGGCGTGTCTGCATCTCTCGCTTTTTACCCACCAtt                                                                                                                                                                                                                    <  1:210810/249‑1 (MQ=255)
                                  gACCGCCAGCATTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCg                                                                                                                                                                                                                 <  1:327901/250‑1 (MQ=255)
                                   aCCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc                                                                                                                                                                                                               >  1:115032/1‑251 (MQ=255)
                                      gCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgct                                                                                                                                                                                                            >  1:651097/1‑249 (MQ=255)
                                        cagcagTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatga                                                                                                                                                                                                          <  1:544789/251‑1 (MQ=255)
                                           cagTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGCGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGGTATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCCTCCTCCCTTTATGTTGTTCCAGGCGTGCCTGCCTCTCTCGCTTTTTACCCACCCTTTCTCCgctgct                                                                                                                                                                                                         >  2:522435/1‑244 (MQ=255)
                                           cagTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAg                                                                                                                                                                                                       <  1:67471/251‑1 (MQ=255)
                                             gTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGcc                                                                                                                                                                                                     >  1:36839/1‑251 (MQ=255)
                                                gCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGc                                                                                                                                                                                                   >  1:61034/1‑250 (MQ=255)
                                                 cgccAGGCCAGTGGACGACATGGCCGACCCCGCCATCTCATCGGGAACGTTCAGCCCCCGAGATGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCATGTTGAATTGCTGTTTGCTTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAg                                                                                                                                                                                                 <  2:473663/249‑1 (MQ=255)
                                                  cccggcccgaggggtggCATGGACACACCCGCCATCGCAGCGGGAACGTCCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGGGTATGCAGGTTTAATTGCTGTTCGCGTGATGGGCCAAAGGCCACGCGTTTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTTCCCACCATTTCGCCGATGATAGACGCa                                                                                                                                                                                                  <  2:574950/234‑1 (MQ=255)
                                                    cAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCATTTTACCCACCATTTCGCCGATGATAGCCGCAGAtt                                                                                                                                                                                              >  1:584936/1‑251 (MQ=255)
                                                     ggcgcgagggagggacgtgcggaccccGCCAGCGCGGCGGGACCGTTAAGCCCCAGAGAGGCCAGCCCGCTGCCGGTCCAGGTCTCTGGATTATGCAGGTTGAATGGCTGTTCGCGTGATGTTCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCTATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTg                                                                                                                                                                                             <  2:674581/227‑1 (MQ=255)
                                                                cGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGCGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGATGGCCCACTGCGCCATCCTCTATTTATGTTTTTCCAGGAGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGCTCGCCGCAGATTGCAAATTTTc                                                                                                                                                                                    >  2:41052/1‑249 (MQ=255)
                                                                   cATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGACAGCACGATGCCGGTCCAGGTCTCCTGCGTCTGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGt                                                                                                                                                                               >  1:607066/1‑251 (MQ=255)
                                                                    aTGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGt                                                                                                                                                                               >  1:69926/1‑250 (MQ=255)
                                                                      ggcccccccctccATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAgg                                                                                                                                                                            <  2:638460/240‑1 (MQ=255)
                                                                       gACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGt                                                                                                                                                                           >  1:548429/1‑251 (MQ=255)
                                                                            aCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCAAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCATCGGAGGTTTCGt                                                                                                                                                                      >  2:614584/1‑251 (MQ=255)
                                                                             cccGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGtt                                                                                                                                                                     >  1:492075/1‑251 (MQ=255)
                                                                               cGCCATCGCAGCGGGAACGTTCAGCCCCATAGAGGCCAGCACGCTGTCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCTAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGc                                                                                                                                                                   <  1:124310/251‑1 (MQ=255)
                                                                                tccataacAGCGGGGACGTTCAGCCCCAGAACGACCAGCACGCTGCCGGTACAGGTCTCCTGAGTCTGCAGGTTGAATTGCTGATCGCGTGATGTGCCAAAGGACACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGGGCCATCCTCTCTTTATGTTGTTCAAGGCGGGGATGCTTCTCTCGCTTTTGACCAACCATTTAGCCGATGACAGCCGCAGCTTGCCAGATTTCAGCGTATGTTTCGTTg                                                                                                                                                                    >  2:126072/8‑249 (MQ=255)
                                                                                gCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCg                                                                                                                                                                  >  1:371182/1‑251 (MQ=255)
                                                                                ccctgcgctgtgggAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCGGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTTGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCGTCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCg                                                                                                                                                                  <  2:680353/240‑1 (MQ=255)
                                                                                      gCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACtt                                                                                                                                                            >  1:629279/1‑251 (MQ=255)
                                                                                          cGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTa                                                                                                                                                           >  1:409995/1‑248 (MQ=255)
                                                                                           gggAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAagc                                                                                                                                                        >  1:507925/1‑250 (MQ=255)
                                                                                                cggtCAGCCCCAGAGAGGCCAGCACGCTGCAGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGt                                                                                                                                                  <  2:36839/248‑1 (MQ=255)
                                                                                                    aaGCCCCAGAGAGGTCCGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTGGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCg                                                                                                                                              <  2:670907/250‑1 (MQ=255)
                                                                                                      gCCCCAGAGAGGCCAGCTCGCTTCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGGGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgc                                                                                                                                            <  2:574080/251‑1 (MQ=255)
                                                                                                      gCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgc                                                                                                                                            >  1:484194/1‑251 (MQ=255)
                                                                                                       ccccAGAGAGTCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACGGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgcg                                                                                                                                           <  2:629279/251‑1 (MQ=255)
                                                                                                       ccccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgcg                                                                                                                                           >  1:601652/1‑251 (MQ=255)
                                                                                                        cccAGAGAGGCCAGCACGCTGCCGGTCCAGTTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGtt                                                                                                                                                                          <  2:465669/219‑1 (MQ=255)
                                                                                                        cccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGtt                                                                                                                                                                          >  1:465669/1‑219 (MQ=255)
                                                                                                         ccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                                                                                                                       >  1:605016/1‑173 (MQ=255)
                                                                                                          cAGAGAGACCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGACCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGcgcgg                                                                                                                                        >  2:626754/1‑250 (MQ=255)
                                                                                                           agagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGAt                                                                                                                                       >  1:193939/1‑251 (MQ=255)
                                                                                                            gagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGCTGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGcgcgctt                                                                                                                                      >  2:563708/1‑248 (MQ=255)
                                                                                                              gagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtg                                                                                                                                    <  1:237902/251‑1 (MQ=255)
                                                                                                                   cAGCACGCCGCAGGTCCAGGTCACCTGAGAATGAAGGTTGAATTGCTGTTCGCGTGATGAGCCAAAGGCCACGAGTGTCACTTAGGGAAGCTGGCTAGACCACCGAGCCATCCTCTCTTTAAGTTGTTCCTGGCGTGCCTGCGGCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAACTTTTCATCGTAGGATTCGTTGCGGGACTTCAGCAGACGTACCGgcgagtttatctgc                                                                                                                                >  2:221865/1‑239 (MQ=255)
                                                                                                                    aGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGACATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCCTCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGAAGGTTTCGTTGCGGGACTTAAGAAGCATTACAGACGAGATTTACTGCa                                                                                                                               >  2:492452/1‑250 (MQ=255)
                                                                                                                    aGCACGCTGCCGATCCAGGCCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCa                                                                                                                               >  2:183811/1‑250 (MQ=255)
                                                                                                                     gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAg                                                                                                                             <  1:369228/251‑1 (MQ=255)
                                                                                                                                  ccAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAcc                                                                                                                 <  1:626754/250‑1 (MQ=255)
                                                                                                                                     ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCg                                                                                                             <  1:659377/251‑1 (MQ=255)
                                                                                                                                          ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATgg                                                                                                        >  2:542221/1‑251 (MQ=255)
                                                                                                                                              aGTATGCAGGTTAAATTGCTGTTCGCGTGATGTGCCAAAGGCAACGCGTGTCACTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGGTGTAACCCCCCGCATGGCGAc                                                                                                    >  2:68949/1‑251 (MQ=255)
                                                                                                                                                    cAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGGCCCTTTGGGAAGCGGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGt                                                                                              <  1:66500/251‑1 (MQ=255)
                                                                                                                                                     cgggttaagtgctggccgcgTGAGGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGTGGGCCCACTGCGCCATCCGCTCTTTATGTTGTTCCAGGCGTGCCTGCATCGCTCGCTTGTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCTCATGGCGACTGCTGTc                                                                                             <  2:665944/235‑1 (MQ=255)
                                                                                                                                                     aGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTc                                                                                             >  1:18990/1‑251 (MQ=255)
                                                                                                                                                      ggTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGGTTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGCCTGCTGTCg                                                                                            >  2:498572/1‑251 (MQ=255)
                                                                                                                                                      ggTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCg                                                                                            >  1:198660/1‑251 (MQ=255)
                                                                                                                                                        ttGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCCTCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGCTTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGAGTGGCGACTGCTGTCGGc                                                                                          >  1:649540/1‑251 (MQ=255)
                                                                                                                                                           aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCCGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAAt                                                                                       >  1:253915/1‑251 (MQ=255)
                                                                                                                                                              tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgat                                                                                                                                                                                                            >  1:2030/1‑131 (MQ=255)
                                                                                                                                                              tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgat                                                                                                                                                                                                            <  2:2030/131‑1 (MQ=255)
                                                                                                                                                               gCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGAt                                                                                    <  1:439372/250‑1 (MQ=255)
                                                                                                                                                                cTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTCAGCAGCAGTACCGGCGCGCTTTGCTGCAAGGCGATGTAACCCCCCGCATGGCGACTGCTGTCGGCAATGATc                                                                                   >  1:382053/1‑250 (MQ=255)
                                                                                                                                                                  gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAg                                                                                 <  1:542221/250‑1 (MQ=255)
                                                                                                                                                                     cgtgtgttgtgCCAACGGCCACTCGTGTCCTTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTGATGTTGGCCAAGGCGTGCCTGCATCTCTCGCTGTTTACCCTCCATTTCGCCGATGATAGCCGCAGATTGCAAATGTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTTCCGGCGCGATTTGCTGCAAGGCGATGTAACCCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCt                                                                             <  2:45339/244‑1 (MQ=255)
                                                                                                                                                                      gcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTg                                                                            <  1:446166/251‑1 (MQ=255)
                                                                                                                                                                         ggATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGtt                                                                         <  1:131543/250‑1 (MQ=255)
                                                                                                                                                                           aTGTGCCAACTGCCACGAGTGTCCCTTTGGGAAGATGGCTGGCCCACTGCACCATCCTCTCTTTATGTTGATCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGCTAGCCGCCGATTGCAAATTTTCAGCGTAGGTTTCGTTGCCGGACTTAAGCATCAATACCGGCGCGCTTTGCTGCAAGGCGATGTAAACCCCCGCCTGGCGACTGCTGTCGGCAATGCTCAGATCTGGtttt                                                                       >  2:625926/1‑250 (MQ=255)
                                                                                                                                                                            tgtgCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTc                                                                       >  1:103028/1‑250 (MQ=255)
                                                                                                                                                                                ccAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGc                                                                  >  1:427989/1‑251 (MQ=255)
                                                                                                                                                                                  aaaGGCCACGCGTGTCCCTTGGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCGTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGTGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCgg                                                                <  2:103028/251‑1 (MQ=255)
                                                                                                                                                                                    tggCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCCCTGCCCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCTGCGTAGGTTTCTTTGCGGGACt                                                                                                                                                             <  2:441508/155‑1 (MQ=255)
                                                                                                                                                                                    aGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACt                                                                                                                                                             >  1:441508/1‑156 (MQ=255)
                                                                                                                                                                                      gCCACGCGGGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGTTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAATCCCCCGCATGGCGACTGCGGTCGGTAATGATCAGGTCTGGTTGCAGAGCGGCAAt                                                            <  1:225338/251‑1 (MQ=255)
                                                                                                                                                                                           gcgTGTCCCTGTGGGACGCGTTCTGGCCCTCTGCGCCATCCTCTCTTTATGTTGTTCCAGGAGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGTCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGTGCGGCAATGGCtt                                                       <  2:670623/251‑1 (MQ=255)
                                                                                                                                                                                             gtgtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGGAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGGCTGGTTTCAGAGCGGCAATGGCTTcc                                                     >  2:362525/1‑251 (MQ=255)
                                                                                                                                                                                              ggtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCa                                                    <  1:183811/250‑1 (MQ=255)
                                                                                                                                                                                                 cccTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTc                                                                                                                                                                        <  1:329359/132‑1 (MQ=255)
                                                                                                                                                                                                 cccTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTc                                                                                                                                                                        >  2:329359/1‑132 (MQ=255)
                                                                                                                                                                                                 cccTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGc                                                 >  2:332783/1‑251 (MQ=255)
                                                                                                                                                                                                 cccTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGc                                                 >  1:270015/1‑251 (MQ=255)
                                                                                                                                                                                                     gtGGGAATCGGGCTGGCCCACTGCGCCATCCTCCCTTTATTTTTGTCCCGGCGTGCCTGCATCTCTCTCTTTTTACCCCCCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCTGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCTACTGCTGTCGGCAATGGTCAGGTCGGGTTTCAGAGCGGCAATGGCTTCCAGGCTCg                                              <  2:384130/249‑1 (MQ=255)
                                                                                                                                                                                                         gAAGCTGGCTGGCCCACTGCGCAATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGAAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGCCGACTGCTGTCGGCAATGATAAGGTCTGGTTTCAGAGCGGCGATGGCTTCCAGGCTCGGccgc                                         >  2:485059/1‑248 (MQ=255)
                                                                                                                                                                                                                gctggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCCTCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAATGCTGTCGGCAATGATCAGGTATGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGcgcgcgcg                                  >  2:567197/1‑251 (MQ=255)
                                                                                                                                                                                                                  tggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGtt                                >  1:378017/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTc                               <  1:226354/251‑1 (MQ=255)
                                                                                                                                                                                                                    gcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCCATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTcc                              >  2:40197/1‑251 (MQ=255)
                                                                                                                                                                                                                        aCTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACg                          >  1:434112/1‑251 (MQ=255)
                                                                                                                                                                                                                            cgcCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGATTTGTACCCACCATTTCGCCGATGATAGCCGGAGATTGCAAATTTTCAGCGTAGGGTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTGGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGAGCAGGTCGTGTATCAGAGCGGCAATGGGTTCCAGGCGCGGCTGCGCGCGGGTTCCGACGGACt                      <  2:96043/251‑1 (MQ=255)
                                                                                                                                                                                                                            cgcCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCGCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACt                      <  1:431458/251‑1 (MQ=255)
                                                                                                                                                                                                                               cATCCTCTCTATATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAATGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGCCggtacacc                   >  2:467403/1‑245 (MQ=255)
                                                                                                                                                                                                                                aTCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGCCGGAATGccc                  >  1:56306/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cctcgttagtgtgttCCAGGCGTGCATGCATCTTTCGCTTTGTCCCCTCCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGATGGACTGCCACGGt              <  2:611141/240‑1 (MQ=255)
                                                                                                                                                                                                                                         ttATGTTGTTCCAGGCGTGCCTGCATCTCTAGCTTTTTACCCACCATTTCGCCGATGCTAGCTGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGTCTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCCGGTTCCGACGGACTGCCACGGTTTc           >  2:256827/1‑249 (MQ=255)
                                                                                                                                                                                                                                                gttCCAGGCGTGCCTGCATCTCTCGCTTTGTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCTGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCGGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGcgca  <  1:41264/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
CGCGATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCA  >  NC_000913/4513594‑4514084

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: