Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,351,519 C→T T48T (ACG→ACA dcuS ← sensor histidine kinase DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,5190CT99.0% 346.0 / NA 102T48T (ACG→ACAdcuSsensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (43/58);  minor base G (0/1);  total (43/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 102 total aligned reads displayed.

CATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAA  >  NC_000913/4351272‑4351752
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
cATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGTTTACCGTTTCCTCTATGGCCTTGATGCCACTCTCCTGTGGTTTTTTCTGCATGCCCTGACGGATTTACGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTc                                                                                                                                                                                                                                        >  2:542211/1‑251 (MQ=255)
cATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTc                                                                                                                                                                                                                                        >  1:10561/1‑251 (MQ=255)
    cAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATAGCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTGCCGCGATGGCCGGGATGCCACTCGCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTTGCTAGCCCATCTCTTGTCata                                                                                                                                                                                                                                     <  2:96505/250‑1 (MQ=255)
     aGCGCTTTAAGGCTGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATc                                                                                                                                                                                                                                   >  1:288436/1‑251 (MQ=255)
      gcgcTTTAAGGATGTCATCACCTGTAAATGGCTGACCAATACGCTGGGCTTCAGTATGCGAGTAGCGAAGACTTTTCATATCGGTAACGACAATAAACAGCATATCGTTGCGTTTGCTTACGGCTTCCGCGATGGCCTGGATGCCTCTCTCCTGCGGTTTTTTCTGCAATCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGCCATATCa                                                                                                                                                                                                                                  <  1:416339/251‑1 (MQ=255)
       cgcTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCAc                                                                                                                                                                                                                                 >  1:456606/1‑251 (MQ=255)
        gcTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATCTCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTATTTTCTGCAAGCCCTGATGGTTTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTCTTGGCTAGCCCATCTCGTGTCATATCAc                                                                                                                                                                                                                                 >  1:476281/1‑250 (MQ=255)
           ttAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTg                                                                                                                                                                                                                               <  1:151229/249‑1 (MQ=255)
           ttAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATAGCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGa                                                                                                                                                                                                                              <  1:456988/250‑1 (MQ=255)
             aaGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCata                                                                                                                                                                                                                                     >  1:128191/1‑241 (MQ=255)
             aaGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGAGGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCGGTCGGATGTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGGGCCTTGTTGGCTAGCCCATCTCGTGTCata                                                                                                                                                                                                                                     <  2:128191/241‑1 (MQ=255)
              aGGATGTCATCACCTTTAAATGGCTGACCAAGACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGGTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACTCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTg                                                                                                                                                                                                                          <  2:563758/251‑1 (MQ=255)
               ggATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGc                                                                                                                                                                                                                         >  2:40585/1‑251 (MQ=255)
               ggATGTCATCACCTTTAAATGGCTGACCAATACGCGGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTGGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGc                                                                                                                                                                                                                         <  2:288436/251‑1 (MQ=255)
                gATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCg                                                                                                                                                                                                                        <  1:198345/251‑1 (MQ=255)
                  tGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCg                                                                                                                                                                                                                        >  1:210801/1‑249 (MQ=255)
                    tcatcaCCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCgaga                                                                                                                                                                                                                     >  1:486704/1‑250 (MQ=255)
                         cccTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATAGCACTGATTTGCGAGAAGTaaa                                                                                                                                                                                                               <  2:697318/250‑1 (MQ=255)
                             ttAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACGCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                             <  1:288073/249‑1 (MQ=255)
                                   ggCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGcac                                                                                                                                                                                                     <  1:333002/251‑1 (MQ=255)
                                       gACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACACTAAACAGCAGATCGTTCCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCATGCCCTGACGGATTTCCGGCGAGTCGGCGCGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGCTTTGCGAGAAGTAAATCAGATGCaccacc                                                                                                                                                                                                 >  1:36259/1‑251 (MQ=255)
                                       gACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGGAGTAAATCAGATGcacca                                                                                                                                                                                                   <  2:463205/249‑1 (MQ=255)
                                       gACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGcaccac                                                                                                                                                                                                  <  1:353554/250‑1 (MQ=255)
                                          cAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCGCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGGGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAAt                                                                                                                                                                                              <  1:57999/251‑1 (MQ=255)
                                           aaTACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATa                                                                                                                                                                                             <  2:6593/251‑1 (MQ=255)
                                                    ggCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCgaga                                                                                                                                                                                    <  1:158827/251‑1 (MQ=255)
                                                       ttCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAAc                                                                                                                                                                                  >  2:550945/1‑250 (MQ=255)
                                                        ccAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCGCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATGTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTGGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCCCCAATAGCACCGAGAAc                                                                                                                                                                                  <  2:4954/248‑1 (MQ=255)
                                                                     aGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTGTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGAGGCACCACCAATAGCACCGAGAACAGTACCGCACTGat                                                                                                                                                                    <  2:142684/250‑2 (MQ=255)
                                                                       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                                                                                                                                          <  2:558474/251‑1 (MQ=255)
                                                                                              gACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                          >  1:616265/1‑251 (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:250309/249‑1 (MQ=255)
                                                                                                                                                 gATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGCTTTCCGGCGAGTCGGCGAGGGCACGCGCCACTGCCAGTGCCTTGTGGTCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGCAGTAAATCAGTTGCACCACCAATAGCACCGAGAACAATACCGCCCTGATCATTAAGATCACTGGGGTACTCAATTTCATCGGCCGTTTGCGGAACATGCGGTAGGGAGATGCATGTCtcataacct                                                                                       >  2:273559/1‑246 (MQ=255)
                                                                                                                                                      cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                             >  1:406531/1‑128 (MQ=255)
                                                                                                                                                      cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                             <  2:406531/128‑1 (MQ=255)
                                                                                                                                                      cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCtt                                                                                  <  2:531591/251‑1 (MQ=255)
                                                                                                                                                            ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtg                                                                               <  2:10561/248‑1 (MQ=255)
                                                                                                                                                                cGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                        <  1:313909/251‑1 (MQ=255)
                                                                                                                                                                cGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                        <  2:210801/251‑1 (MQ=255)
                                                                                                                                                                        tcggCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTa                                                                 <  2:456606/247‑1 (MQ=255)
                                                                                                                                                                         cTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTa                                                                 >  1:564373/1‑249 (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                                                           >  1:101820/1‑251 (MQ=255)
                                                                                                                                                                             aaGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                           <  2:465059/251‑1 (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                                              <  2:178281/238‑1 (MQ=255)
                                                                                                                                                                                                       ggCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAgg                                              >  1:178281/1‑238 (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)
                                                                                                                                                                                                                       ccTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATAGCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATTAATGTCTCATCAGCTTCCTTGTGTGACAAAGTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGAtt                 <  2:322319/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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
CATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAA  >  NC_000913/4351272‑4351752

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