Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,510,751 G→C 11.1% S236R (AGC→AGG fecE ← ferric citrate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,510,7510GC11.1% 283.1 / 2.7 105S236R (AGC→AGGfecEferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (26/62);  new base C (13/2);  total (40/65)
Fisher's exact test for biased strand distribution p-value = 4.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

AAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGC  >  NC_000913/4510513‑4510988
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 
caGCCTGTACAGCTGTGTGTTGCGTTGATTATGAACCTGTCAGCCGAAAGCAAGCGGATGGACGATGAGTATTGGTAATCGTTCAGATTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                                                                                                                                                                                                   <  1:223480/238‑1 (MQ=255)
aagcctgtaacggctgggTTTGGGTTGATTATGAACCTGGCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                                                                                                                                                                                                   <  2:787134/223‑1 (MQ=255)
aaGCCTGTAAAGCTTTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGGCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                                                                                                                                                                                                   <  1:563175/239‑1 (MQ=255)
    cTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCt                                                                                                                                                                                                                                                 <  1:968134/237‑1 (MQ=255)
    cTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCt                                                                                                                                                                                                                                                 <  1:56417/237‑1 (MQ=255)
            cttgggTTTGCGTTTATTAGGAACCTGTCCGCCTAAAGCAAGCGGATGGACGGTGAAGATTGGTAATTCTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGtt                                                                                                                                                                                                                                       <  2:231080/235‑1 (MQ=255)
            ctctGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCCCGCTGAATACTGt                                                                                                                                                                                                                                        >  1:252975/1‑238 (MQ=255)
            ctctGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGGATTTCCGCTTCCACGCTGAATACTGt                                                                                                                                                                                                                                        >  1:173806/1‑238 (MQ=255)
            ctctGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGGTCTTACGGCCGGTGAAATCTCCCTCCTTAGGCACATAGGCCTCCCAGAAACCGGCCCGGGGGGGATTTTCCCCTCCACGCGGaaaccagtt                                                                                                                                                                                                                                       >  2:1037897/1‑232 (MQ=255)
                 gTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCCTTAAGCCCCTCGGCCTTCCAGATACCGGCCCCGGGGGGTTTTCCCCGTCCACCCttgaaaccgttttctc                                                                                                                                                                                                                                  >  2:11741/1‑224 (MQ=255)
                 gTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGTGGTGTATTTCCGATTTCCCCCTGAAAACTCTTTTCa                                                                                                                                                                                                                                   >  1:656475/1‑238 (MQ=255)
                 gTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                                                                                                                                                                                                  <  1:421351/239‑1 (MQ=255)
                 gTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGAAACCGGCTCGGGGGGGATTTCCGCTTCCCCGCTGAATatggtcttcgc                                                                                                                                                                                                                                  >  1:231080/1‑229 (MQ=255)
                 gTTTGCGTTGATTAAGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                                                                                                                                                                                                  <  1:1037897/239‑1 (MQ=255)
                   ttGCGTTGATTATGAACCGGTCAGGCAAAAGCAAGCGGATGGACGATGAGGATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTACCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGc                                                                                                                                                                                                                                 <  1:1030908/238‑1 (MQ=255)
                    tGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCCCTTCCACCCTGAAAACTGTTTTCAGCaa                                                                                                                                                                                                                              >  1:127492/1‑239 (MQ=255)
                      cGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCGCGGGGTGGATTTCCCCTCTCACCCTCAATACCGTTTTCatccaa                                                                                                                                                                                                                              >  2:631871/1‑233 (MQ=255)
                      cGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTAAGGCCTGGGCAATCTACCCCATTATGCCACTCGGCCTGCCCGATCCCCCCTCGGGGGGTTTTTCCCCCTCCCCCCTgcaaacatctttcaacaccc                                                                                                                                                                                                                             >  2:1017118/1‑220 (MQ=255)
                              aTGAACCTGGCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtc                                                                                                                                                                                                                    <  1:747094/239‑1 (MQ=255)
                               tGAACCTGTACGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtca                                                                                                                                                                                                                   <  1:585555/239‑1 (MQ=255)
                                 aaCCTGTCAGACTAAAGCAAGCGGATGGGCGATGAGTATTGGTAAGCTTTCAGAGGCCGGAAAAGTTCAGCCCCAGTCGGAGCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGGGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCAATAGGCACATCGGCCTGCCAGATGCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                                                                                                                                                                                                 <  1:655023/239‑1 (MQ=255)
                                  aCCTGTCAGCCTAAGGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGGCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                                                                                                                                                                                                 <  1:584882/238‑1 (MQ=255)
                                   ccTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCTTTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                                                                                                                                                                                               <  1:11741/239‑1 (MQ=255)
                                     tGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGCCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGTCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCCGGGGGGATTTTCGCCTTCCCCCTCAATACGGTCCTCAGCAACCCTGGGGtcgtaccct                                                                                                                                                                                                             >  1:660710/1‑232 (MQ=255)
                                        cAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACCCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCtctt                                                                                                                                                                                                          <  1:631871/239‑1 (MQ=255)
                                             tAAAGCAAGCGGATGGACGAGGAGTATTGGTAACCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGgt                                                                                                                                                                                                    <  1:390191/239‑1 (MQ=255)
                                                   aagcGGATGGACGTTGGGTATTGGTAATCTTTCAGGGTCCGGAAAAGTTCAGCCCCAGGTTGAACAGGCTTGGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTTCCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCt                                                                                                                                                                                             <  2:252975/239‑1 (MQ=255)
                                                    agcctgtGGTCGGTGAGTTTTGGTAACTTTTCGGGGGCCCGAAAAGTTAAGCCCGCGCTTTGACAGGGTTGCTGGCGCCAGTCCAGTTTCATTCAGCTGTGGGTTGGGACTTACGGCCTGTGCAATCTACTTTATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACTCTGAATACGGTTCTCATCACTCCTGGGTTGATCACCTCTTCTGGTGTGCCt                                                                                                                                                                                             <  2:776752/232‑1 (MQ=255)
                                                    agcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCtt                                                                                                                                                                                            <  1:1003420/239‑1 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACCTCGGCCTGCCCGATACCGGCTCGGGGGGTATTTCCCCTCCCCCCCGGCATACTGTTCTCACCACACCCGGGGTCCTCCCCCCTCTTCGGGTGCCtt                                                                                                                                                                                            >  2:217440/1‑238 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTCCGCTTCCACGCTGCATACTGTTCTCCGCAATCCTGGGGGCATCACCCCTTCTTGTGTGCCtt                                                                                                                                                                                            >  1:931741/1‑238 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGGATTTCCGGTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGGCCTCACCCCTTTTGGGTGTGctt                                                                                                                                                                                             >  1:616107/1‑236 (MQ=255)
                                                         aTGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTgcgc                                                                                                                                                                                        <  1:786697/238‑1 (MQ=255)
                                                           ggACGATTAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTTAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCa                                                                                                                                                                                      <  1:269452/238‑1 (MQ=255)
                                                                 tgggTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt                                                                                                                                                                                <  1:379448/235‑1 (MQ=255)
                                                                         ggtattcttgtcgggTCCGGAAAATTTCAGCCCCAGTCGGAACAGGTTTGCGGGCGCCGGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATTCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTg                                                                                                                                                                       <  2:417573/225‑1 (MQ=255)
                                                                            aatttttgaggggacGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTGTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTGTGCCa                                                                                                                                                                    <  2:83097/225‑1 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCCGATACCGGCTCGGGGTGTTTTTTCCCTTCCACCCTGCATACTGTTTTCCGCCATCCCGGGGTCCTCACCTCTTCTTGCGTGCCTTGCGcccacaaaatcccgtttgcca                                                                                                                                                                    >  2:197439/1‑219 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTGCGCTTTCACCCTGGATACCGTTCTCCGCAATCCTGGGGGCCTTCCCCCTTCCCGTGGTGCCTTCGCCATAAAAATTCCCTTTggcca                                                                                                                                                                   >  2:136262/1‑235 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTTCCCTTCCACGCTGAAAACCGTTCTCAGCAATCCCTGGGTCCTCACCCCTTTTTGTGTGCCTTGCGCCCTAACACTTCCGTTTTccct                                                                                                                                                                   >  1:899113/1‑237 (MQ=255)
                                                                                tttCCGGGTTCGGAAAATTTCCGCCCCTGTTTGAACAGGCTTGCGGGCGGCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTTCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTac                                                                                                                                                                <  2:130859/239‑1 (MQ=255)
                                                                                   cAGGGTTCGTAAAAGTTCACGCCCTGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCGGTCGTGGGTTGGTTCTTACGGCCTGTGCAATCTGCCTCATTACGCACATCGGCCTGCCAGATCCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTaccac                                                                                                                                                             <  2:337819/239‑1 (MQ=255)
                                                                                        tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCTGCCTGCCAGATACCGGCTCGGGGGTTATTTCCGCTTCCCCCCTGCATACTGTTCTCAGCACACCTGGGGGCACAACCCCTTTTGGGGTGCCTTGTGCCCAAAAAAGTCCGTTTTCCCTACCCACCCGtt                                                                                                                                                        >  2:223683/1‑239 (MQ=255)
                                                                                        tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACctc                                                                                                                                                                                                            <  1:638375/189‑1 (MQ=255)
                                                                                         ccGTAAAAGGTCAGGCGCAGTCTGAACCGGCTTGCTGGCACCTGTCCAGTTTCATTCAGTCGTGGTTTGGTTCCTACGGGCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:77626/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACCGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:217440/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                                                                                                                                                                                            <  2:261602/157‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:638815/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTAATTCCGCTTCCACGCTGAATa                                                                                                                                                                                                                                            >  1:261602/1‑157 (MQ=255)
                                                                                          cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCCGGGGGGATTTTCCCTTCCACCCTGCATACCGTTTTCCGCAAATCCTGGGTCATCCACCCCTCTGGTGGTCCTTTTGCCCTAAAATGTTCGTTTTCCCTTTcccccaactat                                                                                                                                                      >  2:240405/1‑234 (MQ=255)
                                                                                          cGGAAAAGGTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGTCCTTCCAGATACCAGACCAGGGTGCAATTTCGCTTCCACCCTGGATCCCGTTCTCAGCAACCCTGGGGTCAACACCCGCTCTTGTGTGCCTTTAGCCATAAAATCTACCGTTGCCATTAACACCCGGTGa                                                                                                                                                      >  2:538930/1‑239 (MQ=255)
                                                                                            gaaattttctcccccCGTCTGAACCGGCGTTCTGGCGCCCGTCCAGTTTAACTCATGCGTGGTTTGGTTCTTCCGGCCTGTGCAATCTACCTCATTAGTCCCATCCGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATc                                                                                                                                                    <  2:942269/228‑1 (MQ=255)
                                                                                              aaaGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCg                                                                                                                                                   <  1:359677/238‑1 (MQ=255)
                                                                                                           cgtttgaaagggttttcTGGCGCCAGTCTATTTTTTTTCCGTCGTGGTGTCGTTCTTCCGGCCTGGGAAATCTACCTCATTAGGCCCAACGGCCTGCCAGATACGGGCTCGGGGGGTATTTCCGCTTCCGCGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  2:926425/223‑1 (MQ=255)
                                                                                                            ggcTGAACAGGCTTGGTGGCGCCAGTCTAGGTTCATTGAGGCGTGGTTTGGTTCTTACGGGCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTGCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  1:1038488/236‑1 (MQ=255)
                                                                                                            gTCTGAACCGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCtt                                                                                                                                                                                            <  1:1096688/183‑1 (MQ=255)
                                                                                                            gTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTCCGCTTCCACGCTGCATACTGTTCTCAGCCATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCCTTACCACCCGGTGATCGCAGTACCGGCTAg                                                                                                                                      >  1:907149/1‑237 (MQ=255)
                                                                                                              cTGAACAGGCTTTGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  1:1011995/236‑1 (MQ=255)
                                                                                                               tGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                                                                                                                                                                                                  <  1:418378/145‑1 (MQ=255)
                                                                                                                     ggtttGCTGGGTCCAGTCCAGTTTCCTTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCCCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAgg                                                                                                                          <  1:486506/236‑1 (MQ=255)
                                                                                                                     gggtgtgcgtcccccgttccgtTTTATTCAGTTGTTGTTTGGGTTTTACGGGCTGTGCAATCTTCCTGAATAGGCAAATCGGCCTGCCAGATTCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAgg                                                                                                                          <  2:907149/219‑1 (MQ=255)
                                                                                                                     agggacagggCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAgg                                                                                                                          <  1:969246/231‑1 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCGGCAATCCTTGGGTCATAACCTCTTCTGGGGGGCCTTGCGCCATAAc                                                                                                                                                                                  >  2:1096688/1‑183 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGGATTTTCCCTTCCACGCCGCAAACCGTTTCCAGGAATCCCGGGGTCACAACCTCTTCTGGGGTGCCTTTCGCCCTAACATGTCCGTTTTCCATTTCCCCCCGTTTATCGCCGTAcccacccacccccataaagg                                                                                                                          >  2:539187/1‑221 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGGATTTCCGCTTCCACGCTGGATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTCCGCCCTAACATGTCCGTTTGCCATTaccac                                                                                                                                                             >  1:936610/1‑204 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTc                                                                                                                        <  1:897187/238‑1 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTc                                                                                                                        <  1:539187/238‑1 (MQ=255)
                                                                                                                           cTGGCTCCATTCCAGTTTCATTCAGGCGTGTTTTGGTTCTTACGGCCTGTCCAATCTACCCCATTAGGCCCATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTgc                                                                                                                    <  2:616107/239‑1 (MQ=255)
                                                                                                                           cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCg                                                                                                                                                   <  1:170447/209‑1 (MQ=255)
                                                                                                                                        attttctttcAGTCTTTGTTTTGTTCTTTCGGCCTGGGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATCCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGacca                                                                                                       <  2:931741/237‑1 (MQ=255)
                                                                                                                                               ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGGGTGCCTTGCGCCATAACATGTCCGTTTGCCATTTCCACCAGTTGACCGCACGACCGGCTAGCCTGACTAAGGGCGTGCAGCAAAGCCGCCACAGttt                                                                                                 >  1:413564/1‑238 (MQ=255)
                                                                                                                                                    tctttttttggttCTTTCGGCCTTTGTAATCTCCCTCATTACGCACCACGGCCGGCCAGATTCCGGTCCTGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCATCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTGCCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCt                                                                                           <  2:899113/233‑1 (MQ=255)
                                                                                                                                                                ttttcaggccggttgaactcgccTTATTAGGCACAGCGGCCTGGCAGGTCCCGGCTCGGGGTGGAGTTCCGCTTCCAGGGTTAATACTGGTCTCAGCAATCCTTGGGTCATCACCTCTTCTGGTGGGCCTTGCGCCATAACATGTCCGTTTGCCATTTACACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGtt                                                                               <  2:1022169/218‑1 (MQ=255)
                                                                                                                                                                tCTTACGGCCTGTGCAATCTACCCCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGtt                                                                               <  1:240405/239‑1 (MQ=255)
                                                                                                                                                                  ttgcGGTCTGTGCAATCTTCCGCATGAGGCCCATCGGCCTGCCAGGTACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCAGCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                             <  2:950454/236‑1 (MQ=255)
                                                                                                                                                                  ctACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                             <  1:605661/238‑1 (MQ=255)
                                                                                                                                                                   tACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCCCTTCTGGTGTGCCTTGCGCCCTAACATGTCCGGTTTCCCTTACCACCAGTTTGTCGCCGTACCGGCCAACCCGATTTAAGGCCTGCAGCACAGCGACCCCCGTTTTCCCCtcaatctcgggtttg                                                                             >  2:907422/1‑224 (MQ=255)
                                                                                                                                                                            tgCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTCTTTCCGCTTCCCCCCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCCGCCAGCCAGGATTAAGTCGCGCAGACCAGCGACCACCGTTTTTCCCCTAGGCCgggggtgcccccacaaa                                                                    >  2:907225/1‑223 (MQ=255)
                                                                                                                                                                            tgCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACACAGACCACCGTTTTCCCCCTAGTCCGGGGTTCTCCCCTCAAc                                                                    >  1:956046/1‑238 (MQ=255)
                                                                                                                                                                               aaTCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATCCTGTTCGCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCa                                                                <  1:270355/239‑1 (MQ=255)
                                                                                                                                                                                      cTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACCTGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAACCTGATTAAGGTCGTGCAACCCAGCGACCACCCTTTTCCCCTTAGTCCGGAGTTTGCCCAaccacccccaccagt                                                          >  1:1061476/1‑229 (MQ=255)
                                                                                                                                                                                         ctttggCCCACCGGCCTGCCTGATTCCGGCTCGGGGTGTTTTTCCGTTTCCACGCTGAATACTGTTCTCAGCAAGCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCAc                                                      <  2:127492/235‑1 (MQ=255)
                                                                                                                                                                                          ttAGGCACGTCGGCCTGCCAGATACCCGCTCGGTGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTgcgc                                                                                                                                                                                        <  2:251435/109‑1 (MQ=255)
                                                                                                                                                                                            agggcaatgcgCCTGCCAGATTACGGCTCGGGGTGTATTTCCGGTTCCACGCTGAATACTGTTCCCAGGAATCCTGGGGTCATCACCTCTTCTGGTTTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCAc                                                      <  1:125316/226‑1 (MQ=255)
                                                                                                                                                                                             ggCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTgg                                                  <  1:907422/239‑1 (MQ=255)
                                                                                                                                                                                               cacaTCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCCTCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCCCCCGTTGATCGCAGTACCGGCCAGCCTGATTTAAGGCGTTCAGCACAGCAACCACCGTTTTCCCCCGAGGCCCGGGGTCGCCCACCaaaccccacaaggtcccccgggg                                                >  2:328440/1‑218 (MQ=255)
                                                                                                                                                                                                  aTCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCCCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATAAGGTCCACCCgggggtt                                             >  1:752536/1‑234 (MQ=255)
                                                                                                                                                                                                   tCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCCATCCTTGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGCACTGCAGCACAGCGACCACCGCTTTCCCCTCAGTCCGGGGGTGCCCCATCAACAGCAACAGGTCCCCCTGGggattg                                            >  1:45236/1‑239 (MQ=255)
                                                                                                                                                                                                    cGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACTCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGGCTGCCTTGCGCCCTAACTTGt                                                                                                                                                                              >  1:251435/1‑109 (MQ=255)
                                                                                                                                                                                                     ggttgtccagttccacgCTCGGTGTGTTTTTCCGCTTCCACGTTGAATTCTGCTCTCAGAAATCGTGGGGTCATCACCTCGTCTGGTGTCCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCACATCAACCGCATCAGGTCCACCTGGtgatt                                             <  2:660710/220‑1 (MQ=255)
                                                                                                                                                                                                            cAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAGCCTGCTTAAGGTCGTTGAGCACAGCGACCCCCCTTTTCCCCTGCGTCCGGGGTTCGCCCCCCAAACCCACCAGGCCCACCCGGTGAATGATATCaaaa                                   >  2:417032/1‑237 (MQ=255)
                                                                                                                                                                                                                     gCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGGGTTCGCCCATCAACCGCATCAGGTCCACCTGGGGATTGATATCAAGCTAGGGGGtt                          >  1:23114/1‑239 (MQ=255)
                                                                                                                                                                                                                            gttgaTTTCCGCTTCCACGCTGAATTCTGTTCCCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                    <  1:851532/234‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtATTTCCGTTTCCAAGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCCCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                    <  1:1017118/238‑1 (MQ=255)
                                                                                                                                                                                                                              gtATTTCCGCTTCCACGCTTAATACTGTTCCCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcatc                  <  1:350299/238‑1 (MQ=255)
                                                                                                                                                                                                                              gtATTTCCGCTTCCACGCTGAATTCTGTTCTCAGCAATCCTGGGGTCATCCCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcatc                  <  1:277376/238‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGgga                                               >  1:224593/1‑200 (MQ=255)
                                                                                                                                                                                                                                         tCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGTGGATTGATATCAAGATAGGTGGTTGGCTCATCAGGTAATACAAc      >  1:823093/1‑239 (MQ=255)
                                                                                                                                                                                                                                         tCCACGCTGAATACTGGTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCa                                                                         <  2:281757/172‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTGAATACTGTTCCCAGCAATCCTGGGGTCATCGCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGc  <  1:567046/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 
AAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGC  >  NC_000913/4510513‑4510988

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