Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,351,590 C→T V25I (GTC→ATC)  dcuS ← sensor histidine kinase DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,5900CT98.1% 346.2 / ‑4.7 103V25I (GTC→ATC) dcuSsensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (46/55);  minor base G (0/2);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 103 total aligned reads displayed.

CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGA  >  NC_000913/4351343‑4351833
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
cGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAt                                                                                                                                                                                                                                                   >  1:600266/1‑251 (MQ=255)
 gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGCGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                  >  2:320711/1‑251 (MQ=255)
    gACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCGGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAg                                                                                                                                                                                                                                               <  2:600266/251‑1 (MQ=255)
      cTTTGCATATCGTTAACGACCATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                                                                                                 >  2:146727/1‑247 (MQ=255)
              atCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTgg                                                                                                                                                                                                                                     <  1:180078/251‑1 (MQ=255)
               ccGGTAACGCCAATAACCAGCAGAGCGTTGGGTTTGCGTACGGCTGCCGCGAGGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                                                                                                    <  2:133565/250‑1 (MQ=255)
                cggcaaCGACAATTAACAGCAGATCGTTGCGTCTGCGTACGGCTTCCGCGACGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCATGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGGTGCACCACCAATAGCACAGAGAACAGGACCGCACTGGTCATTAAGATCACTGTGGt                                                                                                                                                                                                                                    >  2:543262/5‑250 (MQ=255)
                cGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                                                                                                    >  1:352509/1‑250 (MQ=255)
                 ggTAACGACAATAATCAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTGTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGCAAATCAGATGCACCATCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTAc                                                                                                                                                                                                                                  <  1:653573/251‑1 (MQ=255)
                  gTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCATTAGCCCCGAGAACAGTACCGCACTGATCATTAAGACCACTGTGGTACt                                                                                                                                                                                                                                 >  2:295030/1‑251 (MQ=255)
                  gTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACt                                                                                                                                                                                                                                 <  1:320711/251‑1 (MQ=255)
                    aaCGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTc                                                                                                                                                                                                                                >  1:725482/1‑250 (MQ=255)
                      cGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGACACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAt                                                                                                                                                                                                                             >  2:635522/1‑251 (MQ=255)
                       gACAATAAACCGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCCCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAACCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGCTGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                                                                                                            >  2:697471/1‑251 (MQ=255)
                       gACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                                                                                                            >  1:616265/1‑251 (MQ=255)
                       gACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                                                                                                            <  2:558474/251‑1 (MQ=255)
                         cAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTc                                                                                                                                                                                                                          >  1:81588/1‑251 (MQ=255)
                          aaTAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCa                                                                                                                                                                                                                         <  1:279478/251‑1 (MQ=255)
                             aaccaacagATCGTTGCGTTGGCGGGCGGCTACCGGGATGGCCTGGATGCCACTATCCTGCAGGTTTTACTGCAAGCCCTGGCGTATGTCCGGCGAGTCGGCAAGGTTACGCGGCACTGCCAGTTCATTGGTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCTCCATCAATAGCTCCGAGAACAGTGCCGCACTCAGCATTAAGATCTCTTTTGTACTATATTTCATCg                                                                                                                                                                                                                      <  2:704404/245‑1 (MQ=255)
                                cagcagATCGTTGCGTTGGCGTACGGCTCCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGGCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCGCTGATCATTAAGATCACTGTGGTACTCAATGTCATCGGAc                                                                                                                                                                                                                   <  2:504365/251‑1 (MQ=255)
                                        cGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGt                                                                                                                                                                                                           <  2:53887/251‑1 (MQ=255)
                                         gTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGt                                                                                                                                                                                                           >  1:156079/1‑250 (MQ=255)
                                           tgcgtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTaa                                                                                                                                                                                                         <  1:551362/250‑1 (MQ=255)
                                              gtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAg                                                                                                                                                                                                                                               <  1:236541/209‑1 (MQ=255)
                                              gtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAg                                                                                                                                                                                                                                               >  2:236541/1‑209 (MQ=255)
                                                ttgcgtACGGCTTCCGCGATGGCCTGGATGCCTCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGc                                                                                                                                                                                                    <  1:543137/250‑1 (MQ=255)
                                                ttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGc                                                                                                                                                                                                    <  2:486704/250‑1 (MQ=255)
                                                 tgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTTGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAATATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCg                                                                                                                                                                                                   <  1:61040/250‑1 (MQ=255)
                                                  gcgtCCGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTGTGCTTAACATGCGGt                                                                                                                                                                                                 <  1:610629/251‑1 (MQ=255)
                                                    gtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAATAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGt                                                                                                                                                                                                 >  2:125850/1‑249 (MQ=255)
                                                      aCGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGGTTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACCTGCGGTAggg                                                                                                                                                                                             >  2:218961/1‑251 (MQ=255)
                                                       cGGCTTCCGCGATGACCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGACATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGCGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGc                                                                                                                                                                                            >  1:514949/1‑251 (MQ=255)
                                                          cTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaat                                                                                                                                                                                         >  1:254093/1‑251 (MQ=255)
                                                             ccGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatga                                                                                                                                                                                       >  1:87537/1‑250 (MQ=255)
                                                                    ggCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtctc                                                                                                                                                                                <  1:146727/250‑1 (MQ=255)
                                                                     gCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCgg                                                                                                                                                                                                                     <  2:574082/212‑1 (MQ=255)
                                                                     gCCTGGATGCCACTCTCATGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCgg                                                                                                                                                                                                                     >  1:574082/1‑212 (MQ=255)
                                                                         ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatca                                                                                                                                                                            <  1:9036/249‑1 (MQ=255)
                                                                         ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatca                                                                                                                                                                            <  1:250309/249‑1 (MQ=255)
                                                                          gATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGCTTTCCGGCGAGTCGGCGAGGGCACGCGCCACTGCCAGTGCCTTGTGGTCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGCAGTAAATCAGTTGCACCACCAATAGCACCGAGAACAATACCGCCCTGATCATTAAGATCACTGGGGTACTCAATTTCATCGGCCGTTTGCGGAACATGCGGTAGGGAGATGCATGTCtcataacct                                                                                                                                                                         >  2:273559/1‑246 (MQ=255)
                                                                               cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCtt                                                                                                                                                                    <  2:531591/251‑1 (MQ=255)
                                                                                     ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtg                                                                                                                                                                 <  2:10561/248‑1 (MQ=255)
                                                                                         cGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                                                                                                          <  2:210801/251‑1 (MQ=255)
                                                                                         cGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                                                                                                          <  1:313909/251‑1 (MQ=255)
                                                                                                 tcggCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTa                                                                                                                                                   <  2:456606/247‑1 (MQ=255)
                                                                                                   tGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGACTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGACCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGc                                                                                                                                                >  2:439896/1‑251 (MQ=255)
                                                                                                    tcAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCGCCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTaa                                                                                                                                                  <  1:56310/247‑1 (MQ=255)
                                                                                                    gCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAg                                                                                                                                                 <  2:184108/249‑1 (MQ=255)
                                                                                                     cAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAt                                                                                                                                              >  2:448060/1‑251 (MQ=255)
                                                                                                      aaGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                                                                                                             <  2:465059/251‑1 (MQ=255)
                                                                                                      aaGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                                                                                                             >  1:101820/1‑251 (MQ=255)
                                                                                                        gCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAt                                                                                                                                           <  1:484806/251‑1 (MQ=255)
                                                                                                        gCCCGGACGGATTTCCGGCGAGTCTGAGAGGGTACGCGCCTCTGCCAGTGCCGTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACTGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGGCTCATCAGCTTCCTGGTGTGACAAATTTCTTAAGCATTAt                                                                                                                                           <  2:279461/251‑1 (MQ=255)
                                                                                                          ccTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtct                                                                                                                                         >  2:392730/1‑251 (MQ=255)
                                                                                                             gACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtg                                                                                                                                      <  1:392339/251‑1 (MQ=255)
                                                                                                              aCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtga                                                                                                                                     <  2:595146/251‑1 (MQ=255)
                                                                                                                       cGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATTAGGCggg                                                                                                                            >  1:179995/1‑251 (MQ=255)
                                                                                                                                ggCGAGGGTCCGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTcga                                                                                                                    >  2:226177/1‑248 (MQ=255)
                                                                                                                                ggCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAgg                                                                                                                                >  1:178281/1‑238 (MQ=255)
                                                                                                                                ggCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAgg                                                                                                                                <  2:178281/238‑1 (MQ=255)
                                                                                                                                ggCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAtttca                                                                                                                    >  2:547114/1‑247 (MQ=255)
                                                                                                                                  ggAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATGTGCGAGAAGTAAATCAGATGCACCACCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAgg                                                                                                                 <  2:664978/250‑1 (MQ=255)
                                                                                                                                    aGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAgggg                                                                                                               >  1:453431/1‑250 (MQ=255)
                                                                                                                                            cgcCACGGCCAGTGCTTTGTTGGCTAGCCCATCTCGTGTCATATCACTGACTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACAGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATGTATTAAGCATTATCTCTGCTGAGGCGGGTAATTCAAAGGGAGTAAGAAt                                                                                                       <  2:563081/251‑1 (MQ=255)
                                                                                                                                                ccTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTAGTTCTCAATGTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCAGCAGCTTCCTTGTGTTACAAATGTCTGAAGCATTATCTCTGATGAGGCTGGTAATGCAAAGGGAGTAAGAATGAtt                                                                                                   <  2:36259/250‑1 (MQ=255)
                                                                                                                                                       gTGCCTTGTTGGCGAGCCGATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCTCCAATAGCACCGAGAACAGGTCCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACATTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCtata                                                                                            <  1:543262/251‑1 (MQ=255)
                                                                                                                                                              gttgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAgggg                                                                                      <  2:616265/250‑1 (MQ=255)
                                                                                                                                                                tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGaa                                                                                   <  1:302664/250‑1 (MQ=255)
                                                                                                                                                                                tgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                                                                                                             <  1:642939/176‑1 (MQ=255)
                                                                                                                                                                                tgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                                                                                                             >  2:642939/2‑177 (MQ=255)
                                                                                                                                                                                   cATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGTAGTAAGAATGTTTTGGTTTATTTGGGAAGCGAATCTGGCAACGGAAAc                                                               >  2:31662/1‑251 (MQ=255)
                                                                                                                                                                                    atatCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGAATCTGGCAACGGAAACt                                                              >  1:345409/1‑251 (MQ=255)
                                                                                                                                                                                     tatCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCACAGGGAGGAAGAATTATTTGGTATATAGGGGGAGAGAATCTGGCAACGGaaaatg                                                             >  1:255285/1‑248 (MQ=255)
                                                                                                                                                                                         aCTGATTTGCGAGAAGAAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAGTTCACAGGGGGGCAAAATGATTGGCTATATAGGGGGGGAGACGCTTGCAACGGaaaacgcccg                                                         >  2:385631/1‑244 (MQ=255)
                                                                                                                                                                                             aTTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCt                                                     <  2:81588/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTAt                                                 <  1:392835/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTAt                                                 <  1:542211/251‑1 (MQ=255)
                                                                                                                                                                                                   gagaAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAGTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAATAATGATTGGCTATATAgg                                                                                        <  2:662232/211‑1 (MQ=255)
                                                                                                                                                                                                    agaAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGAATCTGGCAACGGAAAATGCCAGGGCTGTATGa                                               >  2:245604/1‑250 (MQ=255)
                                                                                                                                                                                                     gaAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAg                                             <  1:218961/251‑1 (MQ=255)
                                                                                                                                                                                                      aaGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCCGTGCTGTATGa                                               >  1:258211/1‑248 (MQ=255)
                                                                                                                                                                                                        gTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGAGCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACGCTGGCAACGGAAACTGCCAGTGCTGTATGAAGAtt                                          <  2:725482/251‑1 (MQ=255)
                                                                                                                                                                                                              cAGAGGCACCACCAATAGCACCGAGAACAGGACCGCACTGATCATTAAGATCCCTGTGGGACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATCAATGTCTCATCAGCTTCCTGGTGTGACAATTTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAGAGGGGATGAGACTCTGGCAACAGAAACTGCCAGTTCTGTATGAAGATTTCggg                                     <  2:476281/250‑1 (MQ=255)
                                                                                                                                                                                                                   gCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTAt                                <  2:152089/250‑1 (MQ=255)
                                                                                                                                                                                                                    caccacCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGGAAAAGATATGGCAACGGAAAATGACCGTTCTGGATTAAGATTCCTGGGcctagc                              >  1:442890/1‑246 (MQ=255)
                                                                                                                                                                                                                     accaccAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCCCTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGAAATTTCAAGGGGGTAAGAATGATTGGCTATATAGGGGAAGAGAATCTGGCAAAGGAAACAGCCCGGGGTGTATGAAGATTACCGGGCTAttct                             >  2:730300/1‑248 (MQ=255)
                                                                                                                                                                                                                      ccaccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCtt                            <  2:262587/251‑1 (MQ=255)
                                                                                                                                                                                                                       caccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATAACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCCGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTAAAAGGGAGGAAGAATGCTTTGCCATATAGGGGAAGAGAATCTGGAAACGGGAACCGGCCGGGCTGGATGAAGGTTTCGGGGCCATggct                            >  1:680953/1‑247 (MQ=255)
                                                                                                                                                                                                                           aaTAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGc                       <  1:40585/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCg                      <  1:385631/251‑1 (MQ=255)
                                                                                                                                                                                                                               gccccGAGAACAGTACCGCATTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGAGTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCGGTATGAAGATTCCGGGGCTATGCTTATAGCGATa                   <  2:518478/248‑1 (MQ=255)
                                                                                                                                                                                                                                   cGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGACTGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGCTGAGGCGGGTCATTCAAAGGGAGTAAGAATGATTGGCTATATAAGGGGAGAGACTCTGGCAAACGAAAAAGCCAGTGCTGTATGAAGATTCCGGGGGTATGGTTATAGCGATAAtt                >  1:722088/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gagaAAAGAACCGCACTGATCATTAAGATCACTGTGGCACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGGATTCCCAAGGAGTAAGAATGATTTGCTTTATAGGGGAAGAGACTCTGGCAAAGGAAACTGCCAGTGCGGTATGAAAATTCCGGGGCTATGCTTTTTGCGATAATCac              >  2:434500/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gaACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAGGGAGGAAGAAAGATTGGCTTTATCGGGGGAGAGACTCTGGCAACGGGAAATGCCAGTGCTGTATGGAGATTCCGGGGCTATGCTTATAGCGATAATCATAc            >  1:706534/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaCAGTACCGCACTGATCATTAAGATCCCGGTTGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAGCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACGGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATAc            <  2:680953/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTATGGCAACGGGAACTGCCCGTGCTGTTTGAAGATTCCCGGGCTATGCTTATAGCGATAATCATACtg          >  1:119994/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aCCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAAGGAGTAAGAATGATTGGCTATATTGGGGAAGAGACTCTGGCAACGGAAAATGCCAGTGCTGTATGAAGATTTCGGGGCTATGCTTATAGCGATaa                  >  1:581367/1‑238 (MQ=255)
                                                                                                                                                                                                                                             aCCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATaa                  <  2:581367/238‑1 (MQ=255)
                                                                                                                                                                                                                                             aCCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACtgatga      <  1:42028/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATgagaga  <  1:448060/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGA  >  NC_000913/4351343‑4351833

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