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→AGG) fecEferric 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.

AAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGC  >  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)
                    tGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCCCTTCCACCCTGAAAACTGTTTTCAGC‑aa                                                                                                                                                                                                                              >  1:127492/1‑239 (MQ=255)
                      cGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCGCGGGGTGGATTTCCCCTCTCACCCTCAATACCGTTTTCatccaa                                                                                                                                                                                                                              >  2:631871/1‑233 (MQ=255)
                      cGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTAAGGCCTGGGCAATCTACCCCATTATGCCACTCGGCCTGCCCGATCCCCCCTCGGGGGGTTTTTCCCCCTCCCCCCTgcaaacatctttcaacaccc                                                                                                                                                                                                                             >  2:1017118/1‑220 (MQ=255)
                              aTGAACCTGGCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGtc                                                                                                                                                                                                                    <  1:747094/239‑1 (MQ=255)
                               tGAACCTGTACGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGtca                                                                                                                                                                                                                   <  1:585555/239‑1 (MQ=255)
                                 aaCCTGTCAGACTAAAGCAAGCGGATGGGCGATGAGTATTGGTAAGCTTTCAGAGGCCGGAAAAGTTCAGCCCCAGTCGGAGCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGGGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCAATAGGCACATCGGCCTGCCAGATGCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGtcatc                                                                                                                                                                                                                 <  1:655023/239‑1 (MQ=255)
                                  aCCTGTCAGCCTAAGGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGGCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGtcatc                                                                                                                                                                                                                 <  1:584882/238‑1 (MQ=255)
                                   ccTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCTTTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCAc                                                                                                                                                                                                               <  1:11741/239‑1 (MQ=255)
                                     tGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGCCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGTCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCCGGGGGGATTTTCGCCTTCCCCCTCAATACGGTCCTCAGC‑AACCCTGGGGtcgtaccct                                                                                                                                                                                                             >  1:660710/1‑232 (MQ=255)
                                        cAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACCCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCtctt                                                                                                                                                                                                          <  1:631871/239‑1 (MQ=255)
                                             tAAAGCAAGCGGATGGACGAGGAGTATTGGTAACCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑Ggt                                                                                                                                                                                                    <  1:390191/239‑1 (MQ=255)
                                                   aagcGGATGGACGTTGGGTATTGGTAATCTTTCAGGGTCCGGAAAAGTTCAGCCCCAGGTTGAACAGGCTTGGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTTCCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCt                                                                                                                                                                                             <  2:252975/239‑1 (MQ=255)
                                                    agcctgtGGTCGGTGAGTTTTGGTAACTTTTCGGGGGCCCGAAAAGTTAAGCCCGCGCTTTGACAGGGTTGCTGGCGCCAGTCCAGTTTCATTCAGCTGTGGGTTGGGACTTACGGCCTGTGCAATCTACTTTATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACTCTGAATACGGTTCTCATC‑ACTCCTGGGTTGATCACCTCTTCT‑GGT‑GTGCCt                                                                                                                                                                                             <  2:776752/232‑1 (MQ=255)
                                                    agcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCtt                                                                                                                                                                                            <  1:1003420/239‑1 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACCTCGGCCTGCCCGATACCGGCTCGGGGGGTATTTCCCCTCCCCCCCGGCATACTGTTCTCACC‑ACACCCGGGGTCCTCCCCCCTCTT‑CGG‑GTGCCtt                                                                                                                                                                                            >  2:217440/1‑238 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTCCGCTTCCACGCTGCATACTGTTCTCCGC‑AATCCTGGGGGCATCACCCCTTCT‑TGT‑GTGCCtt                                                                                                                                                                                            >  1:931741/1‑238 (MQ=255)
                                                     gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGGATTTCCGGTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGGCCTCACCCCTTTTGGGT‑GTGctt                                                                                                                                                                                             >  1:616107/1‑236 (MQ=255)
                                                         aTGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTgcgc                                                                                                                                                                                        <  1:786697/238‑1 (MQ=255)
                                                           ggACGATTAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTTAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCa                                                                                                                                                                                      <  1:269452/238‑1 (MQ=255)
                                                                 tgggTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACAt                                                                                                                                                                                <  1:379448/235‑1 (MQ=255)
                                                                         ggtattcttgtcgggTCCGGAAAATTTCAGCCCCAGTCGGAACAGGTTTGCGGGCGCCGGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATTCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTg                                                                                                                                                                       <  2:417573/225‑1 (MQ=255)
                                                                            aatttttgaggggacGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTGTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTGTGCCa                                                                                                                                                                    <  2:83097/225‑1 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCCGATACCGGCTCGGGGTGTTTTTTCCCTTCCACCCTGCATACTGTTTTCCGC‑CATCCCGGGGTCCTCACCTCTTCT‑TGC‑GTGCCTTGCGcccacaaaatcccgtttgcca                                                                                                                                                                    >  2:197439/1‑219 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTGCGCTTTCACCCTGGATACCGTTCTCCGC‑AATCCTGGGGGCCTTCCCCCTTCC‑CGTGGTGCCTT‑CGCCATAAAAATTCCCTTTggcca                                                                                                                                                                   >  2:136262/1‑235 (MQ=255)
                                                                             aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTTCCCTTCCACGCTGAAAACCGTTCTCAGC‑AATCCCTGGGTCCTCACCCCTTTT‑TGT‑GTGCCTTGCGCCCTAACACTTCCGTTTTccct                                                                                                                                                                   >  1:899113/1‑237 (MQ=255)
                                                                                tttCCGGGTTCGGAAAATTTCCGCCCCTGTTTGAACAGGCTTGCGGGCGGCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTTCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTac                                                                                                                                                                <  2:130859/239‑1 (MQ=255)
                                                                                   cAGGGTTCGTAAAAGTTCACGCCCTGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCGGTCGTGGGTTGGTTCTTACGGCCTGTGCAATCTGCCTCATTACGCACATCGGCCTGCCAGATCCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTaccac                                                                                                                                                             <  2:337819/239‑1 (MQ=255)
                                                                                        tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCTGCCTGCCAGATACCGGCTCGGGGGTTATTTCCGCTTCCCCCCTGCATACTGTTCTCAGC‑ACACCTGGGGGCACAACCCCTTTT‑GGG‑GTGCCTTGTGCCCAAAAAAGTCCGTTTTCCCTACCCACCCGtt                                                                                                                                                        >  2:223683/1‑239 (MQ=255)
                                                                                        tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACctc                                                                                                                                                                                                            <  1:638375/189‑1 (MQ=255)
                                                                                         ccGTAAAAGGTCAGGCGCAGTCTGAACCGGCTTGCTGGCACCTGTCCAGTTTCATTCAGTCGTGGTTTGGTTCCTACGGGCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:77626/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACCGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:217440/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                                                                                                                                                                                            <  2:261602/157‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGtt                                                                                                                                                        <  1:638815/238‑1 (MQ=255)
                                                                                         ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTAATTCCGCTTCCACGCTGAATa                                                                                                                                                                                                                                            >  1:261602/1‑157 (MQ=255)
                                                                                          cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCCGGGGGGATTTTCCCTTCCACCCTGCATACCGTTTTCCGCAAATCCT‑GGGTCATCCACCCCTCT‑GGT‑GGTCCTTTTGCCCTAAAATGTTCGTTTTCCCTTTcccccaactat                                                                                                                                                      >  2:240405/1‑234 (MQ=255)
                                                                                          cGGAAAAGGTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGTCCTTCCAGATACCAGACCAGGGTGCAATTTCGCTTCCACCCTGGATCCCGTTCTCAGC‑AACCCTGGGGTCAACACCCGCTCT‑TGT‑GTGCCTTTAGCCATAAAATCTACCGTTGCCATTAACACCCGGTGa                                                                                                                                                      >  2:538930/1‑239 (MQ=255)
                                                                                            gaaattttctcccccCGTCTGAACCGGCGTTCTGGCGCCCGTCCAGTTTAACTCATGCGTGGTTTGGTTCTTCCGGCCTGTGCAATCTACCTCATTAGTCCCATCCGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATc                                                                                                                                                    <  2:942269/228‑1 (MQ=255)
                                                                                              aaaGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCg                                                                                                                                                   <  1:359677/238‑1 (MQ=255)
                                                                                                           cgtttgaaagggttttcTGGCGCCAGTCTATTTTTTTTCCGTCGTGGTGTCGTTCTTCCGGCCTGGGAAATCTACCTCATTAGGCCCAACGGCCTGCCAGATACGGGCTCGGGGGGTATTTCCGCTTCCGCGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  2:926425/223‑1 (MQ=255)
                                                                                                            ggcTGAACAGGCTTGGTGGCGCCAGTCTAGGTTCATTGAGGCGTGGTTTGGTTCTTACGGGCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTGCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  1:1038488/236‑1 (MQ=255)
                                                                                                            gTCTGAACCGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCtt                                                                                                                                                                                            <  1:1096688/183‑1 (MQ=255)
                                                                                                            gTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGTATTTCCGCTTCCACGCTGCATACTGTTCTCAGC‑CATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCCTTACCACCCGGTGATCGCAGTACCGGCTAg                                                                                                                                      >  1:907149/1‑237 (MQ=255)
                                                                                                              cTGAACAGGCTTTGTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                                                     <  1:1011995/236‑1 (MQ=255)
                                                                                                               tGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                                                                                                                                                                                                  <  1:418378/145‑1 (MQ=255)
                                                                                                                     ggtttGCTGGGTCCAGTCCAGTTTCCTTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCCCCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAgg                                                                                                                          <  1:486506/236‑1 (MQ=255)
                                                                                                                     gggtgtgcgtcccccgttccgtTTTATTCAGTTGTTGTTTGGGTTTTACGGGCTGTGCAATCTTCCTGAATAGGCAAATCGGCCTGCCAGATTCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAgg                                                                                                                          <  2:907149/219‑1 (MQ=255)
                                                                                                                     agggacagggCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAgg                                                                                                                          <  1:969246/231‑1 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCGGC‑AATCCTTGGGTCATAACCTCTTCT‑GGG‑GGGCCTTGCGCCATAAc                                                                                                                                                                                  >  2:1096688/1‑183 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGGATTTTCCCTTCCACGCCGCAAACCGTTTCCAGG‑AATCCCGGGGTCACAACCTCTTCT‑GGG‑GTGCCTTTCGCCCTAACATGTCCGTTTTCCATTTCCCCCCGTTTATCGCCGTAcccacccacccccataaagg                                                                                                                          >  2:539187/1‑221 (MQ=255)
                                                                                                                      gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGGGGATTTCCGCTTCCACGCTGGATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTCCGCCCTAACATGTCCGTTTGCCATTaccac                                                                                                                                                             >  1:936610/1‑204 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTc                                                                                                                        <  1:897187/238‑1 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTc                                                                                                                        <  1:539187/238‑1 (MQ=255)
                                                                                                                           cTGGCTCCATTCCAGTTTCATTCAGGCGTGTTTTGGTTCTTACGGCCTGTCCAATCTACCCCATTAGGCCCATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTgc                                                                                                                    <  2:616107/239‑1 (MQ=255)
                                                                                                                           cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCg                                                                                                                                                   <  1:170447/209‑1 (MQ=255)
                                                                                                                                        attttctttcAGTCTTTGTTTTGTTCTTTCGGCCTGGGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATCCCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGacca                                                                                                       <  2:931741/237‑1 (MQ=255)
                                                                                                                                               ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGG‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTTCCACCAGTTGACCGCACGACCGGCTAGCC‑TGACTAAGGGCGTGCAGCAAAGCCGCCACAGttt                                                                                                 >  1:413564/1‑238 (MQ=255)
                                                                                                                                                    tctttttttggttCTTTCGGCCTTTGTAATCTCCCTCATTACGCACCACGGCCGGCCAGATTCCGGTCCTGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCATC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTGCCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCt                                                                                           <  2:899113/233‑1 (MQ=255)
                                                                                                                                                                ttttcaggccggttgaactcgccTTATTAGGCACAGCGGCCTGGCAGGTCCCGGCTCGGGGTGGAGTTCCGCTTCCAGGGTTAATACTGGTCTCAGC‑AATCCTTGGGTCATCACCTCTTCT‑GGT‑GGGCCTTGCGCCATAACATGTCCGTTTGCCATTTACACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGtt                                                                               <  2:1022169/218‑1 (MQ=255)
                                                                                                                                                                tCTTACGGCCTGTGCAATCTACCCCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGtt                                                                               <  1:240405/239‑1 (MQ=255)
                                                                                                                                                                  ttgcGGTCTGTGCAATCTTCCGCATGAGGCCCATCGGCCTGCCAGGTACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCAGCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                             <  2:950454/236‑1 (MQ=255)
                                                                                                                                                                  ctACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCg                                                                             <  1:605661/238‑1 (MQ=255)
                                                                                                                                                                   tACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCCCTTCT‑GGT‑GTGCCTTGCGCCCTAACATGTCCGGTTTCCCTTACCACCAGTTTGTCGCCGTACCGGCCAACC‑CGATTTAAGGCCTGCAGCACAGCGACCCCCGTTTTCCCCtcaatctcgggtttg                                                                             >  2:907422/1‑224 (MQ=255)
                                                                                                                                                                            tgCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTCTTTCCGCTTCCCCCCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCCGCCAGCCAGGATTAA‑GTCGCGCAGACCAGCGACCACCGTTTTTCCCCTAGGCCgggggtgcccccacaaa                                                                    >  2:907225/1‑223 (MQ=255)
                                                                                                                                                                            tgCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACACAGACCACCGTTTTCCCCCTAGTCCGGGGTTCTCCCCTCAAc                                                                    >  1:956046/1‑238 (MQ=255)
                                                                                                                                                                               aaTCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATCCTGTTCGCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCa                                                                <  1:270355/239‑1 (MQ=255)
                                                                                                                                                                                      cTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACCTGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAACC‑TGATTAAGGTCGTGCAACCCAGCGACCACCCTTTTCCCCTTAGTCCGGAGTTTGCCCAaccacccccaccagt                                                          >  1:1061476/1‑229 (MQ=255)
                                                                                                                                                                                         ctttggCCCACCGGCCTGCCTGATTCCGGCTCGGGGTGTTTTTCCGTTTCCACGCTGAATACTGTTCTCAGC‑AAGCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCAc                                                      <  2:127492/235‑1 (MQ=255)
                                                                                                                                                                                          ttAGGCACGTCGGCCTGCCAGATACCCGCTCGGTGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTgcgc                                                                                                                                                                                        <  2:251435/109‑1 (MQ=255)
                                                                                                                                                                                            agggcaatgcgCCTGCCAGATTACGGCTCGGGGTGTATTTCCGGTTCCACGCTGAATACTGTTCCCAGG‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑TTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCAc                                                      <  1:125316/226‑1 (MQ=255)
                                                                                                                                                                                             ggCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTgg                                                  <  1:907422/239‑1 (MQ=255)
                                                                                                                                                                                               cacaTCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCCTCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCCCCCGTTGATCGCAGTACCGGCCAGCC‑TGATTTAAGGCGTTCAGCACAGCAACCACCGTTTTCCCCCGAGGCCCGGGGTCGCCCACCaaaccccacaaggtcccccgggg                                                >  2:328440/1‑218 (MQ=255)
                                                                                                                                                                                                  aTCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCCCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATAAGGTCCACCCgggggtt                                             >  1:752536/1‑234 (MQ=255)
                                                                                                                                                                                                   tCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑CATCCTTGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGCACTGCAGCACAGCGACCACCGCTTTCCCCTCAGTCCGGGGGTGCCCCATCAACAGCAACAGGTCCCCCTGGggattg                                            >  1:45236/1‑239 (MQ=255)
                                                                                                                                                                                                    cGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACTCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGG‑CTGCCTTGCGCCCTAACTTGt                                                                                                                                                                              >  1:251435/1‑109 (MQ=255)
                                                                                                                                                                                                     ggttgtccagttccacgCTCGGTGTGTTTTTCCGCTTCCACGTTGAATTCTGCTCTCAGA‑AATCGTGGGGTCATCACCTCGTCT‑GGT‑GTCCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCACATCAACCGCATCAGGTCCACCTGGtgatt                                             <  2:660710/220‑1 (MQ=255)
                                                                                                                                                                                                            cAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAGCC‑TGCTTAAGGTCGTTGAGCACAGCGACCCCCCTTTTCCCCTGCGTCCGGGGTTCGCCCCCCAAACCCACCAGGCCCACCCGGTGAATGATATCaaaa                                   >  2:417032/1‑237 (MQ=255)
                                                                                                                                                                                                                     gCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGGGTTCGCCCATCAACCGCATCAGGTCCACCTGGGGATTGATATCAAGCTAGGGGGtt                          >  1:23114/1‑239 (MQ=255)
                                                                                                                                                                                                                            gttgaTTTCCGCTTCCACGCTGAATTCTGTTCCCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                    <  1:851532/234‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtATTTCCGTTTCCAAGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCCCCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtca                    <  1:1017118/238‑1 (MQ=255)
                                                                                                                                                                                                                              gtATTTCCGCTTCCACGCTTAATACTGTTCCCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcatc                  <  1:350299/238‑1 (MQ=255)
                                                                                                                                                                                                                              gtATTTCCGCTTCCACGCTGAATTCTGTTCTCAGC‑AATCCTGGGGTCATCCCCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtcatc                  <  1:277376/238‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGgga                                               >  1:224593/1‑200 (MQ=255)
                                                                                                                                                                                                                                         tCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGTGGATTGATATCAAGATAGGTGGTTGGCTCATCAGGTAATACAAc      >  1:823093/1‑239 (MQ=255)
                                                                                                                                                                                                                                         tCCACGCTGAATACTGGTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCa                                                                         <  2:281757/172‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTGAATACTGTTCCCAGC‑AATCCTGGGGTCATCGCCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGc  <  1:567046/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 
AAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGC‑AATCCTGGGGTCATCACCTCTTCT‑GGT‑GTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCC‑TGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGC  >  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: ‑