Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,751 A→G D59D (GAT→GAC)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,7510AG94.4% 417.9 / 9.2 143D59D (GAT→GAC) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (63/72);  total (66/77)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 143 total aligned reads displayed.

GCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAC  >  NC_000913/72504‑72997
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
     aaCGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGt                                                                                                                                                                                                                                                <  1:634120/251‑1 (MQ=255)
     aaCGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGt                                                                                                                                                                                                                                                >  1:378071/1‑251 (MQ=255)
      aCGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTc                                                                                                                                                                                                                                               <  1:953007/251‑1 (MQ=255)
            cgcgctgtcgcggACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGaa                                                                                                                                                                                                                                         <  2:842988/249‑1 (MQ=255)
             gccctggcgcttACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTGAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAAc                                                                                                                                                                                                                                        <  2:725541/244‑1 (MQ=255)
                  gtcgctgACCGCCGGAAAGTTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTTCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgct                                                                                                                                                                                                                                    <  1:588771/250‑1 (MQ=255)
                    cgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggc                                                                                                                                                                                                                                 <  1:543253/251‑1 (MQ=255)
                         aCCGCCGGCAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCg                                                                                                                                                                                                                            >  1:822365/1‑251 (MQ=255)
                               ggaacgccgcccagtgaaCCGCGCCCTTAAATTATCAATCCCCATCGGGCGGGCGATAGCGTGCAGTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGGCAGCGGACGGCGTGACGGCGGGGTCGTGGTGTGATCAACGCCGTCGATAGGCAGCGAACCGCTGGCTGGCGTCAg                                                                                                                                                                                                                        <  2:436913/233‑1 (MQ=255)
                                 aaaGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGCCGGCGTGACGGCGGCATCGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGaaa                                                                                                                                                                                                                     >  2:360433/1‑250 (MQ=255)
                                    gCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGTCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCCGAAAACCg                                                                                                                                                                                                                 >  1:660944/1‑251 (MQ=255)
                                       cGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                                                                              <  1:210217/251‑1 (MQ=255)
                                       cGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCGGGCGGGGGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAACCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGAGCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                                                                              <  2:240590/251‑1 (MQ=255)
                                        gccaggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                                                                              >  1:986626/5‑250 (MQ=255)
                                        gCCCGGTAACAGCGCCAATAAATTATCAATCCCCATCTGGCGGTCGATAGCGTGCATTTTCCCCTGCTGTGCCGCGTTCAGTTTCAAGCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGACCTACGCCATCGATAGTCAGCGAACCCCTGGCTGGCGTCAGAAAACCGGCGa                                                                                                                                                                                                             <  1:510127/251‑1 (MQ=255)
                                         ccatgtaccctcgccaatcaaTTATCAATCCCCATCTGGCGTGCGATAGCGTTCAATGTCGCCTGCTGTGCCGCGTTCAGTTTCAATCTCGGATTGATCCCCAGCCCGATGGTCTGCGCGACCGTCAGGTGACTGCACAGGTTGGTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGGCGAGAGGCAGCGACCCGCTGGCTGGCGTCAGAAATCCGGCGAt                                                                                                                                                                                                            <  2:778446/232‑1 (MQ=255)
                                         ccaggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATCGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGAt                                                                                                                                                                                                            >  2:697980/4‑251 (MQ=255)
                                         ccaggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGGCGTGGTGGGATCAACGCCGTCGATAGTCAGAGAACCGCTGGCTGGGGTAAGAAAACCGGCGac                                                                                                                                                                                                            >  2:724056/4‑250 (MQ=255)
                                           aggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCa                                                                                                                                                                                                          <  1:63614/250‑1 (MQ=255)
                                              tAACCGCGACATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAAACGGCGATCaagt                                                                                                                                                                                                       >  2:489696/1‑249 (MQ=255)
                                               aaCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGAGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGGGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCaagtt                                                                                                                                                                                                      >  1:751007/1‑248 (MQ=255)
                                               aaCCGAGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCaagtt                                                                                                                                                                                                      >  1:195512/1‑248 (MQ=255)
                                                  cgcgCCATTAAAGTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGAGTCAGAAAACCGGCGATCAAGTTcag                                                                                                                                                                                                   <  1:953926/251‑1 (MQ=255)
                                                   gcgcCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcagc                                                                                                                                                                                                  <  1:930543/251‑1 (MQ=255)
                                                       cATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAggg                                                                                                                                                                                              >  1:788689/1‑251 (MQ=255)
                                                          tAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTAAGCAggggac                                                                                                                                                                                           >  1:530863/1‑248 (MQ=255)
                                                           aaaTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTAc                                                                                                                                                                                           <  1:170253/250‑1 (MQ=255)
                                                             aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                                                                                        <  2:207389/251‑1 (MQ=255)
                                                             aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                                                                                        <  1:892727/251‑1 (MQ=255)
                                                             aTTACCAATCCCCTTCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAATCCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                                                                                        <  2:755925/251‑1 (MQ=255)
                                                              ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGCGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGAGGCAGCAAACCGGCGGTCACGTTCAGCAGGGTACtttt                                                                                                                                                                                       >  2:570215/1‑251 (MQ=255)
                                                                aTCAATCCCCATCTGGCGGGCGCTAGCGTGCAGTTTCTCCTGCTGTGCCGCGTTCAGGTTCAATCCCGGAGTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCGGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAc                                                                                                                                                                                     <  2:487185/251‑1 (MQ=255)
                                                                aTCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAc                                                                                                                                                                                     <  1:163138/251‑1 (MQ=255)
                                                                 tCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcc                                                                                                                                                                                    >  1:991509/1‑251 (MQ=255)
                                                                 ccAATCCCCATCGTGCGGGCGATAGCGTGCATTTTCTCCGGCTGTGCCGCGTGCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCGGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACgttttcc                                                                                                                                                                                    <  2:858061/250‑4 (MQ=255)
                                                                  cAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcc                                                                                                                                                                                    <  2:222964/250‑1 (MQ=255)
                                                                  cAATCCCCATCTCGCGGTCGATAGCGTGCATTTTCGCCTGCTGTGCCGCGTTCAGGTTCAATCCCGGATTTAGCCCCAGCCCGATTGTCTTCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAATTTCAGCAGGGTACTTTTAcct                                                                                                                                                                                   <  2:977775/251‑2 (MQ=255)
                                                                       cccATCTGGCTGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGcgcc                                                                                                                                                                              <  1:415159/251‑1 (MQ=255)
                                                                        ccATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAAAAGGTTGTTTTCCTTAACCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGCGTTATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCATACAACCGGCGATCAAATTCAGCAGGGTACTTTTAcctgagccg                                                                                                                                                                             >  1:549354/1‑244 (MQ=255)
                                                                            cTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTg                                                                                                                                                                         >  1:311053/1‑251 (MQ=255)
                                                                              ggcgggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCtt                                                                                                                                                                          <  1:695279/248‑1 (MQ=255)
                                                                              ggcgggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGc                                                                                                                                                                       <  1:841057/251‑1 (MQ=255)
                                                                                 gggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCg                                                                                                                                                                    >  1:892790/1‑251 (MQ=255)
                                                                                 gggcgATAGCGTGCATTTCCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCg                                                                                                                                                                                                                            <  1:35585/195‑1 (MQ=255)
                                                                                  ggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGCTGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGCTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGAGTCAGACAACCGGCAGTCAAGGTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCgg                                                                                                                                                                   >  2:93437/1‑250 (MQ=255)
                                                                                     gatatagTGCAATTTCTCCTGCTGTGCCGCGTTCAGTTCCAATCCCGTATTTAGCCCCAGCCCGATGGACTGCGCGACCTGCAGGTGACTGAACAGGTTGTTGGCCTGAAACAGCAGCGACACCGGACGGCGTGACGGCGGCGTCTTGGTGTGATCAACCCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCTGCGATCAAGTTCATCAGGTTACTTTAACCTGCGCCGCGTGGGCCGAgta                                                                                                                                                                <  2:134599/245‑3 (MQ=255)
                                                                                           gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCa                                                                                                                                                          <  1:215156/251‑1 (MQ=255)
                                                                                             gCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCAcc                                                                                                                                                        >  1:134547/1‑251 (MQ=255)
                                                                                             gCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCAcc                                                                                                                                                        <  1:453700/251‑1 (MQ=255)
                                                                                                   tctccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTcg                                                                                                                                                  >  1:702408/4‑251 (MQ=255)
                                                                                                      ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCTATGTTCTGCGCGACCGTCAGGTGACTGATCAGGTTGTTTTCCTTAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCCCCTGCTcgcct                                                                                                                                               <  1:663798/251‑2 (MQ=255)
                                                                                                      ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTcgccg                                                                                                                                               <  1:663796/251‑1 (MQ=255)
                                                                                                        ggctgTGCCGCGTTCAGGTTCAATCCCGGATGTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGGCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCcgcg                                                                                                                                             <  2:391861/250‑1 (MQ=255)
                                                                                                           ggtgccGCGTTCAGGTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTGTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTc                                                                                                                                          <  2:46414/247‑1 (MQ=255)
                                                                                                                  cgTTCCGTTTCAAGCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGtt                                                                                                                                   <  2:210375/251‑1 (MQ=255)
                                                                                                                  cgTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGtt                                                                                                                                   <  1:50923/251‑1 (MQ=255)
                                                                                                                  cgTTCAGTTTCAATACCGGATTTAGCCCCAGCCCGATGTTCTGCGCGCCCGTCAGGTGGCTGAACAGGTTGTTTTCCTGAAACCGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGCTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGGGGATCACCACCTGCTCGCCGCGTTCTACCGtt                                                                                                                                   >  2:133372/1‑251 (MQ=255)
                                                                                                                                cccGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCAGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCa                                                                                                                      <  1:360433/250‑1 (MQ=255)
                                                                                                                                   ggATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTCCTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCCACTGCTCGCCGCGTTCTCCCGTTAAGCTAAACCGCATCgg                                                                                                                  >  2:183990/1‑251 (MQ=255)
                                                                                                                                   ggATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCgg                                                                                                                  <  1:924442/251‑1 (MQ=255)
                                                                                                                                   ggATTTAGCCCCAGCCCGATGGTCTGCTCGACCGTCAGGGGACTGAACAGGTTGTTTTCCTGAAACAGCGTCGGCACCGGACGGCGTGACGGCGGCGTCGTGGTGGGATCAACTCCGTCGATAGTCAGCGAACCGCTGGCTGGGTTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTGTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCgg                                                                                                                  <  1:909621/251‑1 (MQ=255)
                                                                                                                                          gactcAGCCCGATGTCCTGCGCGCCCTTCAGGGGACTTATCAGGTTGTGTTCCTGAAACAGCCTCGTCACCGGACGGCGGTACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCTTTAAGCTAAAACGCATCGGCAAAtgt                                                                                                           <  2:533411/247‑2 (MQ=255)
                                                                                                                                          gCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgt                                                                                                           >  1:591522/1‑250 (MQ=255)
                                                                                                                                             acAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGCTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGGACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGGTAAGCTAAAACGCACAGGCAAAtgtggg                                                                                                        >  2:613640/2‑247 (MQ=255)
                                                                                                                                                 cccGATGTTCTGCGCGACCGTCAGGTGACTGCACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACTGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGAAAATGTTGGTaa                                                                                                     >  1:332812/1‑250 (MQ=255)
                                                                                                                                                    gATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGcc                                                                                                 >  1:289102/1‑251 (MQ=255)
                                                                                                                                                       gTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGGGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAgg                                                                                              <  2:368970/251‑1 (MQ=255)
                                                                                                                                                         tcttcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCTACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGt                                                                                             <  1:138613/245‑1 (MQ=255)
                                                                                                                                                           tgcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGat                                                                                          >  1:382230/4‑251 (MQ=255)
                                                                                                                                                                gACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcag                                                                                     <  1:735009/251‑1 (MQ=255)
                                                                                                                                                                 aCCGTCAGGTGACTGAACAGGTTGGTTTCCTGAAACAGCATCGACAGCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcagt                                                                                    <  1:8912/251‑1 (MQ=255)
                                                                                                                                                                   cGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTCCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTACAGCCAGGTGGTAtcagtc                                                                                   >  1:737438/1‑250 (MQ=255)
                                                                                                                                                                        ggggacTGAACAGGTAGTTTACCCGCAACAGCATCGTCACCAGACGGCGTGACGGCTGCGACGTGGTGCGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTGCTACCTTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCCGtt                                                                               <  2:281128/244‑1 (MQ=255)
                                                                                                                                                                         gtgacTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGtttt                                                                             >  1:76012/5‑250 (MQ=255)
                                                                                                                                                                            acTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGCCGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCTCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAAGTTGGTAAAGACAGGTGATATCAGACAGTTTTaaa                                                                          >  2:27655/2‑249 (MQ=255)
                                                                                                                                                                               gAACAGGTTGTTTGCCTGAAACAGCATCGGCACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAtct                                                                      <  2:378071/251‑3 (MQ=255)
                                                                                                                                                                                    ggTTGTTTTCCTGAAACGGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCCGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGCTCGCCACCTGCGCGCCGCGTTCTACCGTTAAGCTAAAACGCTTCGGCACATGTTGGTAAAGCCAGGTGATATCAGTCAGATTTAACATCTCGccc                                                                 >  2:18069/1‑251 (MQ=255)
                                                                                                                                                                                     gttgttTTCCTGAAACAGCATCGACACCGGACGGCGAGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTCCTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTCAGCTCAAACGCATCGGCAACAGTTGGTAACGCCCGGAGATATCTGTACGTTTTAACATCTCGcccc                                                                >  2:662714/1‑251 (MQ=255)
                                                                                                                                                                                      ttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCg                                                               >  1:824724/1‑251 (MQ=255)
                                                                                                                                                                                          tttcCTGAAACAGCATCGACACCGGAAGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAACAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTGAAGCCAGGACATATCAGTCAGTTTTAACATCTCCCCCCGGt                                                             >  2:629233/3‑249 (MQ=255)
                                                                                                                                                                                           ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggct                                                                                                                                                                                                                                <  1:560411/84‑1 (MQ=255)
                                                                                                                                                                                           ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGGACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCGTCGGCAACTGTTGGTAACGCCAGGTGATATCAGTCAGTTTTAAAATCTCGCCCCGGTAgg                                                          >  2:25699/2‑250 (MQ=255)
                                                                                                                                                                                            tcCTGAAACAGCATCGACGCCGGACGGCGTTACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGt                                                          <  2:331123/250‑1 (MQ=255)
                                                                                                                                                                                                gAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTc                                                     >  1:328612/1‑251 (MQ=255)
                                                                                                                                                                                                     aGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCa                                                <  1:19534/251‑1 (MQ=255)
                                                                                                                                                                                                           gACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCCAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCCCGGTa                                           >  1:97242/1‑250 (MQ=255)
                                                                                                                                                                                                           gACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTa                                           >  1:337507/1‑250 (MQ=255)
                                                                                                                                                                                                              acCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAc                                        >  1:257442/1‑250 (MQ=255)
                                                                                                                                                                                                                cGGACGGCGTGACTGCGGCGTCGTGGTGTGATCAACGCCGTCGTTAGTCAGCGAGCCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGTTCTCCACCTGCTCGCCTCGTTCTACCGTTAAGCTAAAACGCATCGGCATATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAGCTCGTCCCGGTAGTTTTTCAATCACGGTAAAcagc                                     <  2:596178/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagca                                    <  1:90919/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                     >  1:419630/1‑250 (MQ=255)
                                                                                                                                                                                                                    cGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGAAAACAGCAGa                                  >  1:845781/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaa                                 >  1:321959/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGAGATATCAGTCAGTTTTAAACTCTCGCCCCGGTAGTTTTTCAATCACggataaaagcagca                                 >  2:181801/1‑238 (MQ=255)
                                                                                                                                                                                                                         tGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATCTCAGTCAGTTTTAACATCTCGCCCCGGAAGTTTTTCAATCACGGTAAACAGCAGAAcaca                             >  1:556420/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAt              <  1:114224/249‑1 (MQ=255)
                                                                                                                                                                                                                                               tcaaCGCCGTCGATAGTCAGAGAACCGCTGGCTGGCGTCAGAAAACCGGCGATAAAGTTCAGCAGGGAACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGACAGGGGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCACTCACGGAATACAGCAGAAAACAGTGTAGCAGAAGAATTAAAGCgg      >  2:463193/4‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaCGCCGTCGATAGTCGGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg    >  1:214994/2‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGGTCAGCAGGGTACTTTTACCTGCGCCGCGTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg    <  2:741318/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   cGCCGTCGAGAGTCAGCGAACCGCTGGCTGGCGTCAGAAATCCGGCGAGCAAGTTCAGCAGGGTTCTTTTACCTGCGCCTCTTGGCTCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGAc  <  2:302217/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAC  >  NC_000913/72504‑72997

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