Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,524,596 C→T D139N (GAT→AAT)  yjhH ← putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,524,5960CT98.8% 268.8 / NA 81D139N (GAT→AAT) yjhHputative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (37/43);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524360‑4524842
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
cggggCAAAGTTAGCGCTGGCAGTTATCGCACCGTCGCCGCGCAGCAGCATCGTATTCAGCAAAGGATCATCGTAACCGCAGAATGCCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGGTGGTGTCTTTGATGCCAACGATAGTCTCGTTTTGCAGAGCCGGACGCGTCACGGTTACCGGGGTTAAGTCCTGACCCTTCAGATTCGGAAAGTTGTACa                                                                                                                                                                                                                                          <  2:518769/247‑1 (MQ=255)
cGGAGCAAAGTTAGCGCTGGCGGATATCGCTCCGTCGCCGCTCAGCATCATCGTAGTCAGTAAATGAGCATCGTAACCGCAGAATACCGAAAACGACGGGCGTGCCGACGTAGCTGTGTTGATCATCGTTCGCAAGTGACCAACTCCGTCGGTGGTGTCTTTGTTGGTAAGCATATTCTCTTTTGGCAGAGCCAGACAGGTCAGGGTTTCCGGGGTTAGGGCTTCACCCGGCAGAGCCGGAATTTTGTACa                                                                                                                                                                                                                                          <  2:65336/251‑1 (MQ=255)
  gAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCCTTGTAACCGCAGAATAACGAAAACAACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGTACCGTCAGATTCGGAAAGTTGTACAgg                                                                                                                                                                                                                                        >  2:306721/1‑251 (MQ=255)
  gAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCTCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAgg                                                                                                                                                                                                                                        <  1:264161/251‑1 (MQ=255)
        aggtAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTGGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCAcc                                                                                                                                                                                                                                  <  1:201524/248‑1 (MQ=255)
            aGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGt                                                                                                                                                                                                                               <  1:28056/250‑1 (MQ=255)
               gcTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAggg                                                                                                                                                                                                                           <  1:86101/251‑1 (MQ=255)
                  ggCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTg                                                                                                                                                                                                                         <  2:490021/250‑1 (MQ=255)
                        tATCGCACCTTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTGTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGAGCACCGGTAGGGTGACGCTa                                                                                                                                                                                                                   <  2:337664/250‑1 (MQ=255)
                            gCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGc                                                                                                                                                                                                              <  1:235083/251‑1 (MQ=255)
                             cACCGTCGACGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGATAAGTCCTGACCAGTCAGATTCGGAGAGATGTACAGGATCACCGGTAGGGTGACGCTActgg                                                                                                                                                                                                               >  1:312208/1‑246 (MQ=255)
                                      cgctcagcagCATCGTATTCAGCAAATGCTCATCGTAACCGCAGAATACCGAAAACGACGAGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGAGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAATTCCTGACCCGTCAGATTCGGAAAGTTGTCCAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTgg                                                                                                                                                                                                    >  2:6226/5‑251 (MQ=255)
                                           agcagcATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAAt                                                                                                                                                                                                >  1:310709/1‑250 (MQ=255)
                                           agcagcATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATa                                                                                                                                                                                               <  1:441897/251‑1 (MQ=255)
                                           agcagcATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTAGTACAGGATCACCGGAAGGGTGACGCTACGGGCGATCTGCTGGTAAt                                                                                                                                                                                                >  2:279706/1‑250 (MQ=255)
                                             cagcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgt                                                                                                                                                                                             <  2:310166/251‑1 (MQ=255)
                                               gcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcg                                                                                                                                                                                           <  1:517275/251‑1 (MQ=255)
                                                 aTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         <  2:142223/251‑1 (MQ=255)
                                                   cGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCa                                                                                                                                                                                        >  1:39505/1‑250 (MQ=255)
                                                    gtagtCAGCAAATGATCATCGTAACCGCAGAAGACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAg                                                                                                                                                                                      <  2:549509/247‑1 (MQ=255)
                                                          aGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCg                                                                                                                                                                                <  2:159226/251‑1 (MQ=255)
                                                               aTGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtggtg                                                                                                                                                                           <  1:446183/251‑1 (MQ=255)
                                                               aTGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtggtg                                                                                                                                                                           <  1:364698/251‑1 (MQ=255)
                                                                  atcatcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGCCCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCg                                                                                                                                                                         >  1:353833/1‑250 (MQ=255)
                                                                  atcatcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCgc                                                                                                                                                                        >  1:409900/1‑250 (MQ=255)
                                                                  atcatcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGTCTGTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGa                                                                                                                                                                        <  2:306532/251‑1 (MQ=255)
                                                                   tcatcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGAc                                                                                                                                                                       >  2:437713/1‑251 (MQ=255)
                                                                     ctcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGCCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACtt                                                                                                                                                                     >  2:431896/2‑251 (MQ=255)
                                                                           aCCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCACGATTTCGTGGTGCGACTTTCCagt                                                                                                                                                               >  1:25488/1‑251 (MQ=255)
                                                                           aCCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagt                                                                                                                                                               <  2:467297/251‑1 (MQ=255)
                                                                             cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtg                                                                                                                                                                              <  1:35092/234‑1 (MQ=255)
                                                                             cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtg                                                                                                                                                                              >  2:35092/1‑234 (MQ=255)
                                                                               cAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagtagt                                                                                                                                                            >  1:149609/1‑250 (MQ=255)
                                                                                  aaTACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAggg                                                                                                                                                        >  1:390800/1‑251 (MQ=255)
                                                                                        gAAAACGAAGGGCGTACCGAATTAACTGTGTTGATCATCGTACGCAAGTGACCCACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTAAAGGATAACCGGTAGGGTGACGCCACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGATGCGACTTTCCAGTAGTAGGGGTGGAt                                                                                                                                                  >  2:344018/1‑251 (MQ=255)
                                                                                          aaaCGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATgg                                                                                                                                                >  1:373272/1‑251 (MQ=255)
                                                                                                cgGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGcgacg                                                                                                                                           >  1:493821/1‑250 (MQ=255)
                                                                                                     gTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATAcc                                                                                                                                      >  1:167801/1‑250 (MQ=255)
                                                                                                       aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATc                                                                                                                                   >  1:463034/1‑251 (MQ=255)
                                                                                                       aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATc                                                                                                                                   >  2:573256/1‑251 (MQ=255)
                                                                                                           aCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAgcg                                                                                                                               >  1:76443/1‑251 (MQ=255)
                                                                                                                   gtgtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAg                                                                                                                         <  1:176026/249‑1 (MQ=255)
                                                                                                                    tgtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTCGGCt                                                                                                                      >  2:27524/1‑251 (MQ=255)
                                                                                                                     gtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCtt                                                                                                                     <  1:279706/251‑1 (MQ=255)
                                                                                                                     gtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCtt                                                                                                                     <  2:402915/251‑1 (MQ=255)
                                                                                                                     gtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCt                                                                                                                      <  1:347701/250‑1 (MQ=255)
                                                                                                                      ttGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTg                                                                                                                    >  1:384343/1‑251 (MQ=255)
                                                                                                                       tGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTg                                                                                                                    <  1:132941/250‑1 (MQ=255)
                                                                                                                                 aCGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCAtgctgc                                                                                                          <  1:335123/250‑1 (MQ=255)
                                                                                                                                   gCAAGTGACCAAAGCTGCCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGACAGACGCGTCACGGTTTACGGGGTTAAGTCCTGACCCGTAAGATTCGGAAAGTTGTACCGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTAGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACACTCAGCGCCGTCGGCTTGCGCATGCTGCGcc                                                                                                       >  2:340788/1‑251 (MQ=255)
                                                                                                                                          aCCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTg                                                                                                 >  2:433573/1‑250 (MQ=255)
                                                                                                                                          aCCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTg                                                                                                 >  1:429925/1‑250 (MQ=255)
                                                                                                                                             aaCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCg                                                                                             >  1:126897/1‑251 (MQ=255)
                                                                                                                                                 cTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTc                                                                                         >  2:284291/1‑251 (MQ=255)
                                                                                                                                                           tgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGGCAGATTCGGAACGTTGTACAGGATCACCGGAAGGGAGACGCTACGGGCGATCTGATGGTAATAGTCGTCAAGCTTTCGTGGTGCGACTTTCCAGTAGTAGGTGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGCCCGCTTCGTCAGTGGa                                                                                >  2:517018/1‑250 (MQ=255)
                                                                                                                                                                ttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAgggg                                                                           <  1:239113/250‑1 (MQ=255)
                                                                                                                                                                ttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGGTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAGAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAgggg                                                                           <  2:9060/250‑1 (MQ=255)
                                                                                                                                                                      ccAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCTGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGAc                                                                    <  1:51449/251‑1 (MQ=255)
                                                                                                                                                                        cacggtatTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTAAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAATCGTCAAGATTTCGTGTTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATAAGCGCCGTAGGCATGCGCATGCTGCGCCAGTTTGACCGCGTCGTCAGTGGCAGGGGAACCGACGc                                                                  >  2:145872/6‑251 (MQ=255)
                                                                                                                                                                           gATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCACGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCca                                                                >  1:108010/1‑250 (MQ=255)
                                                                                                                                                                                 ctcGTTTTGCAGAGCCAGACGCGTCACGGATTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGAGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGAAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCcaatcaat                                                          >  1:187969/1‑250 (MQ=255)
                                                                                                                                                                                          cAGAGCAAGACGCGTCACGGATTCCGGGTTTAAGACCTGACCCGTAAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGCCGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATAGCGACGATACCATCAGGGCCGTAGGATTGCGCATGCTGCGACAGGTTGACCCATTCGTCAGTGGAAGGGGAACCGATGCAAATCACTACCGGAACTc                                                >  2:440983/1‑251 (MQ=255)
                                                                                                                                                                                           agagCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCg                                               >  1:6164/1‑251 (MQ=255)
                                                                                                                                                                                              gCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGccc                                            <  2:482861/251‑1 (MQ=255)
                                                                                                                                                                                                  gACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCg                                        <  2:149609/251‑1 (MQ=255)
                                                                                                                                                                                                  gACGCGTCACGGTTGCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCg                                        <  1:487884/251‑1 (MQ=255)
                                                                                                                                                                                                    cgcgTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGAc                                      <  1:509963/251‑1 (MQ=255)
                                                                                                                                                                                                       gTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGCAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAAt                                   >  1:190189/1‑251 (MQ=255)
                                                                                                                                                                                                         cACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATg                                  <  1:573256/250‑1 (MQ=255)
                                                                                                                                                                                                                  cGGGGTTAAGTCCTGACCCGTCAGATTCGGAATGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATTCCGGCACTCGCCCGTCGACAATGGTTACAGctt                        <  1:122425/251‑1 (MQ=255)
                                                                                                                                                                                                                    gggTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttct                      <  1:317219/251‑1 (MQ=255)
                                                                                                                                                                                                                    gggTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttct                      >  2:61934/1‑251 (MQ=255)
                                                                                                                                                                                                                     ggTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCAGCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttctt                     <  1:524156/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggTTAAGTCCTGAACCGTCAGATTCGGAAAGTTGGACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCTCTCGCCCGTCGACAATGGTTACAGcttctt                     <  1:499376/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACACTGGTTACAGCttcttc                    >  1:359505/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATTGTTACAGCTTCTTCGGc                 <  2:390800/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCgg                  <  1:27524/250‑1 (MQ=255)
                                                                                                                                                                                                                                gACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTGCTTCGGCGAGTGcc          <  1:182566/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cccGGCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCAt        <  1:433573/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gTCAGATTCGGAAAGTTGTACAGGATCACCGGGAGGGTGACGCGACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTGCGTCAGTGGAAGGGGAACCGACGCCAATCAAGACCGGCACTCGCCCGTCGACAATGTTTACAGCGTCTTCGGCGAGTGCCATgcg     <  2:353833/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGAGCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTg  <  2:14445/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
CGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524360‑4524842

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