Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,691 C→T E79E (GAG→GAA)  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,6910CT93.9% 443.7 / 12.6 151E79E (GAG→GAA) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (62/77);  total (67/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 151 total aligned reads displayed.

TCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGT  >  NC_000913/72446‑72938
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
gcGAGCGCAGAGAACGGTTCACCGAGAAATAAAATCGGCTGTTCGCGCACCAGACAACGCTCTAACGCCACGCGCTGT‑CGCTGTCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTg                                                                                                                                                                                                                                                    <  2:798219/250‑1 (MQ=255)
          agaACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATc                                                                                                                                                                                                                                          >  2:19534/1‑251 (MQ=255)
           gaACGGTTCACCGACCAATAAAACCGGCGGTTCGCGCACCAGACAACGCGCTAACGCCCCGCGCTGT‑CGCTGCCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCg                                                                                                                                                                                                                                         <  2:42177/251‑1 (MQ=255)
            cgcacttccgatataaaTAAAAGCGGCTGGTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTAGCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATGTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTGGTTTTCCTGAAACAGCATCGa                                                                                                                                                                                                                                        <  2:394046/236‑1 (MQ=255)
              cGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGaca                                                                                                                                                                                                                                      <  1:900230/251‑1 (MQ=255)
                gtgcATCGAGAAAAAAACTCGGCTTCTCGCGCACCACACACCGCGCTTACGCCACTCGCGGT‑CGCAGACCTCCGGAAAGCTCGCCGTGGAACCGCGCCATTAAATTATCAATCCCCATCTGGCTGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTAAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCTACCGTCATGTTACTGAACAGGTTGTTTTCCTGAAACATCATCGACAcc                                                                                                                                                                                                                                    <  2:457535/248‑1 (MQ=255)
                gTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGCAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGCTGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACAcc                                                                                                                                                                                                                                    >  1:651348/1‑251 (MQ=255)
                 ttCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCAACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGAGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGATGGTTTCATGAAACAGCATCGACACCg                                                                                                                                                                                                                                   >  1:865158/1‑251 (MQ=255)
                  tCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCgg                                                                                                                                                                                                                                  >  1:242983/1‑251 (MQ=255)
                    aTCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAc                                                                                                                                                                                                                                >  1:554357/1‑251 (MQ=255)
                      agAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGA‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCAATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGAGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCAAGCCCGATGTTCTTCGCGACCGTCAGGTGCCTGGACAGGTTGTTTTCCTGAAACAGCATCGAACCCGGACgg                                                                                                                                                                                                                              >  2:663798/2‑251 (MQ=255)
                          cAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGACACAGCATCGa                                                                                                                                                                                                                                        >  1:394046/1‑237 (MQ=255)
                               aaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc                                                                                                                                                                                                                     <  1:719161/251‑1 (MQ=255)
                               aaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc                                                                                                                                                                                                                     >  1:160748/1‑251 (MQ=255)
                                aaTCGCCTGGTCGCGCCCCAGACAACGCGCTAACGCCCCGCGCTGT‑CGCTGACCGCCGGAAAGCGCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTTCATTTTCTCCTGCGGTGCCGCGTTCAGTGTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACTCCGGACGGCGTGAcggcg                                                                                                                                                                                                                    <  2:628311/251‑1 (MQ=255)
                                      cTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCACTCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcg                                                                                                                                                                                                              >  1:754285/1‑247 (MQ=255)
                                          gcgcgCACCAGACAACGCGCTAACGCCACGCGCGGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCAGCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt                                                                                                                                                                                                          <  2:53775/250‑9 (MQ=255)
                                           cgcgCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt                                                                                                                                                                                                          <  1:339767/250‑9 (MQ=255)
                                           cgcgCACCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt                                                                                                                                                                                                          <  1:881143/250‑9 (MQ=255)
                                                aCCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACTGCGTGACggcggc                                                                                                                                                                                                                  >  2:656031/1‑237 (MQ=255)
                                                aCCAGACAACGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggc                                                                                                                                                                                                                  <  1:656031/237‑1 (MQ=255)
                                                        acgctcgaaCGCCCCGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccgt                                                                                                                                                                                            <  2:580471/244‑3 (MQ=255)
                                                        aCGCGCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccgt                                                                                                                                                                                            <  1:841658/251‑3 (MQ=255)
                                                         cgcgCTAACGCCCCGCGC‑GTGCGCTGCCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccgtc                                                                                                                                                                                           <  2:440351/251‑4 (MQ=255)
                                                         cgcgCTAACGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTAGTGTGATCAACGccgtc                                                                                                                                                                                           <  2:617957/251‑4 (MQ=255)
                                                               aaCGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGt                                                                                                                                                                                     >  1:378071/1‑251 (MQ=255)
                                                               aaCGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGt                                                                                                                                                                                     <  1:634120/251‑1 (MQ=255)
                                                                aCGCCACGCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTc                                                                                                                                                                                    <  1:953007/251‑1 (MQ=255)
                                                                      cgcgctgt‑cgcggACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGaa                                                                                                                                                                              <  2:842988/249‑1 (MQ=255)
                                                                            gt‑cgctgACCGCCGGAAAGTTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTTCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgct                                                                                                                                                                         <  1:588771/250‑1 (MQ=255)
                                                                               cgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggc                                                                                                                                                                      <  1:543253/251‑1 (MQ=255)
                                                                                          ggaacgccgcccagtgaaCCGCGCCCTTAAATTATCAATCCCCATCGGGCGGGCGATAGCGTGCAGTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGGCAGCGGACGGCGTGACGGCGGGGTCGTGGTGTGATCAACGCCGTCGATAGGCAGCGAACCGCTGGCTGGCGTCAg                                                                                                                                                             <  2:436913/233‑1 (MQ=255)
                                                                                            aaaGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGCCGGCGTGACGGCGGCATCGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGaaa                                                                                                                                                          >  2:360433/1‑250 (MQ=255)
                                                                                                  cGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                   <  1:210217/251‑1 (MQ=255)
                                                                                                  cGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCGGGCGGGGGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAACCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGAGCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                   <  2:240590/251‑1 (MQ=255)
                                                                                                   gccaggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCg                                                                                                                                                   >  1:986626/5‑250 (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)
                                                                                                                        aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                             <  2:207389/251‑1 (MQ=255)
                                                                                                                        aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                                                                                                             <  1:242747/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)
                                                                                                                           aTCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCGGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACttt                                                                                                                             <  2:374772/248‑1 (MQ=255)
                                                                                                                            tCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcc                                                                                                                         >  1:991509/1‑251 (MQ=255)
                                                                                                                            ccAATCCCCATCGTGCGGGCGATAGCGTGCATTTTCTCCGGCTGTGCCGCGTGCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCGGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACgttttcc                                                                                                                         <  2:858061/250‑4 (MQ=255)
                                                                                                                               aTCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcctgc                                                                                                                      >  1:19153/1‑248 (MQ=255)
                                                                                                                                 ccccATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAAACGGCGATCAAGTTCAGCAGGGTACTTTTACCTgcgc                                                                                                                    >  1:801877/1‑251 (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)
                                                                                                                                            gggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCg                                                                                                         >  1:892790/1‑251 (MQ=255)
                                                                                                                                            gggcgATAGCGTGCATTTCCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCg                                                                                                                                                                 <  1:35585/195‑1 (MQ=255)
                                                                                                                                            gggcgATAGCGTGCATTTCCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGAAGGCGTGACGGCGGAATAGTTGTGTGATCTACGCCATCGATAGCCAGCGAACCGCTGGCTGGCg                                                                                                                                                                 >  2:35585/1‑195 (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)
                                                                                                                                                                 ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCTATGTTCTGCGCGACCGTCAGGTGACTGATCAGGTTGTTTTCCTTAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCCCCTGCTcgcct                                                                                    <  1:663798/251‑2 (MQ=255)
                                                                                                                                                                 ccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTcgccg                                                                                    <  1:663796/251‑1 (MQ=255)
                                                                                                                                                                    gctgTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGGTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTcgccgc                                                                                   <  1:686947/249‑1 (MQ=255)
                                                                                                                                                                      ggtgccGCGTTCAGGTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTGTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTc                                                                               <  2:46414/247‑1 (MQ=255)
                                                                                                                                                                          ccGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTAcc                                                                           <  1:703441/251‑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)
                                                                                                                                                                                               gATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGACAGCGAACCGCTGGCTGGCGTCCGAAAACCGGCGATCAAGTTCAGCAGGGAACTTTTACCTGCGCCGCTTGGACCGAGGATCGCCACCTGCTTGCCGCGTTCTACCGTTAAGATCAAACGCATCgg                                                       >  2:720988/1‑250 (MQ=255)
                                                                                                                                                                                                     gactcAGCCCGATGTCCTGCGCGCCCTTCAGGGGACTTATCAGGTTGTGTTCCTGAAACAGCCTCGTCACCGGACGGCGGTACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCTTTAAGCTAAAACGCATCGGCAAAtgt                                                <  2:533411/247‑2 (MQ=255)
                                                                                                                                                                                                     gCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAtgt                                                >  1:591522/1‑250 (MQ=255)
                                                                                                                                                                                                       cccAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGAGTGATCAACGCCGTCGCTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGAGCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGATAGGATAAAACGCATCGGCAAAtgttg                                              >  2:212865/1‑248 (MQ=255)
                                                                                                                                                                                                            cccGATGTTCTGCGCGACCGTCAGGTGACTGCACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACTGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGAAAATGTTGGTaa                                          >  1:332812/1‑250 (MQ=255)
                                                                                                                                                                                                               gATGTTCTGCGCGACCGTCAGGTGACTGAACAGGATGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCACGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAAAGCATCGGCAAATGTTGGTAAAGcc                                      >  2:819980/1‑251 (MQ=255)
                                                                                                                                                                                                                    tcttcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCTACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGt                                  <  1:138613/245‑1 (MQ=255)
                                                                                                                                                                                                                      tgcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGat                               >  1:382230/4‑251 (MQ=255)
                                                                                                                                                                                                                        cgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGACGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGatat                             <  2:746210/250‑1 (MQ=255)
                                                                                                                                                                                                                           gACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAATGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtc                            >  1:32640/1‑249 (MQ=255)
                                                                                                                                                                                                                           gACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcag                          <  1:735009/251‑1 (MQ=255)
                                                                                                                                                                                                                            aCCGTCAGGTGACTGAACAGGTTGGTTTCCTGAAACAGCATCGACAGCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcagt                         <  1:8912/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTCCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTACAGCCAGGTGGTAtcagtc                        >  1:737438/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ggTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGCCGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGtttt                  >  1:507172/1‑251 (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)
                                                                                                                                                                                                                                                   gttCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGt  <  1:833922/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGT‑CGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGT  >  NC_000913/72446‑72938

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: ‑