Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,692,311 T→C 8.4% E99G (GAG→GGG)  bcsB ← regulator of cellulose synthase, cyclic di‑GMP binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,692,3110TC 8.4% 423.3 / 2.4 158E99G (GAG→GGG) bcsBregulator of cellulose synthase, cyclic di‑GMP binding
Reads supporting (aligned to +/- strand):  ref base T (74/61);  new base C (4/12);  total (79/79)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCATTGATCAGTGGTTGCGTTGCTGGCGTCGCCTGCGTCATGAAA  >  NC_000913/3692077‑3692547
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
tGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGTACAGGCAGTAACGATGGCGATAGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCACCGCTGCGCATACCAAACTCAAt                                                                                                                                                                                                                                            >  1:709207/1‑239 (MQ=255)
  gCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCCCTTCCCGGCTGCGCATACCAAACACAAtgc                                                                                                                                                                                                                                          >  2:360590/1‑239 (MQ=255)
  gCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAAt                                                                                                                                                                                                                                            >  1:355027/1‑237 (MQ=255)
  gCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAAt                                                                                                                                                                                                                                            >  1:994063/1‑237 (MQ=255)
      aCAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGCTAAACCTTTAACTGCGACTGGACAGGCCGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATGGCCTTCCTCCCCCCTTCATCGCCTCGCATACCAAACCCAATGCTccc                                                                                                                                                                                                                                      >  2:448469/1‑239 (MQ=255)
         actccccgCCGTTCACGGTTGAAGTCGCTAAAAAACGGTGGGTTAATGGGCATTTGCGCCCGGGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGtt                                                                                                                                                                                                                                   <  2:872698/232‑1 (MQ=255)
         aaCTCCGGCCGTACCCGGTTGAAGTCCTTAATAAACCTGGGGTTAATGGGCATTTGGCCCGGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGtt                                                                                                                                                                                                                                   <  1:472969/239‑1 (MQ=255)
         aaCTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGtt                                                                                                                                                                                                                                   <  2:18996/239‑1 (MQ=255)
              cAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGGGTATTCGAGGTTGAGAATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGt                                                                                                                                                                                                                              >  1:534367/1‑239 (MQ=255)
                   cttcccggtgtaatgccgtaataaACGGTGGGTTAAGGGGCATTTGGCCCAGCGTTTTTTTACCCACCTGTTCTTTGGTCACTGGCAGGACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATg                                                                                                                                                                                                                         <  2:971599/221‑1 (MQ=255)
                      cacGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTTCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGca                                                                                                                                                                                                                       >  1:834091/1‑237 (MQ=255)
                        cggttgacgctcgaaagaaaCGGGGGGTTAATGGGCATTTGCACCAAGGTTTTTTTACCCAACTGGTCTTTGGTCACTGGCAGCATGCCCATCAGTTCATCATTGAGATAACACTTTACCTGCGACTGGAAAGGCAGTAACGATGGTGATGGGGTGTAGTCGAGGTTGAGCATCGCCTTCGTCACCTCTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACg                                                                                                                                                                                                                    <  1:627497/222‑1 (MQ=255)
                        cGGTTGAAGGCGCTAAAAAACTGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACg                                                                                                                                                                                                                    <  2:768021/239‑1 (MQ=255)
                             taaGTCGCTAATAAACTGGGGGTTAATGGGCATTTGCCCCAGCGTTTTTTTACCCAACTGTTCTTTTGTCCCTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCa                                                                                                                                                                                                               <  2:939804/238‑1 (MQ=255)
                             gAAGTCTTTAATAACCAGTTGTTTAATGGGCTTTTTCGCCCGCTTTTTTTTCCCCACCTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCTACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCa                                                                                                                                                                                                               <  2:105431/239‑1 (MQ=255)
                              aaGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTAATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGCTGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCTAACAAAACTCAATGCTGCCGTTCGGGTTAAGGACACGTAacac                                                                                                                                                                                                              >  1:588148/1‑239 (MQ=255)
                                   gCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGCGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGGGTATTCGAGGTTCAGCATCGCCTTCGTCACCACTTCATCGCTGGGCATACCAAGCCCAATGCTGCCGTTCGGGTTAATGCCACAAAGCAcccca                                                                                                                                                                                                          >  2:792778/1‑235 (MQ=255)
                                   gCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAtg                                                                                                                                                                                                          >  1:23316/1‑238 (MQ=255)
                                   gCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAtg                                                                                                                                                                                                          >  1:166537/1‑238 (MQ=255)
                                     taataaACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGtt                                                                                                                                                                                                                                   >  1:271773/1‑211 (MQ=255)
                                     taataaACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGAGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAtgct                                                                                                                                                                                                        >  1:850177/1‑238 (MQ=255)
                                          aaCATTGGGTTAATGGGCATTGGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcgg                                                                                                                                                                                                  <  1:222980/239‑1 (MQ=255)
                                              gTGGGTTAATGGGCTTTTGCGCCCGGGTTTTTTTTCCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCCGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcggcggc                                                                                                                                                                                              <  2:187975/239‑1 (MQ=255)
                                               tGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGATATAAACCTTTAACTGCGACTGGACCGGCAGTAACGATGGCGTTAGAGTGAATTCCAGGTTGGGCACCCCCGTCGTGACCACTTCATCGCTCACCGTACCAAACCCCATGCTGCCGGTCGGGGTAAACCCCCGTGACAACAAGCGCCCcggcggcg                                                                                                                                                                                             >  2:684080/1‑239 (MQ=255)
                                               tGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGtt                                                                                                                                                                                                                                   <  2:262131/201‑1 (MQ=255)
                                               tGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCAt                                                                                                                                                                                                           <  1:25235/225‑1 (MQ=255)
                                               tGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATCTGAGGTTGCGCCTCGCCTTCGTCACCACTTCTCCGCTGCCCATACCACACTCCATGCTGCCGTTCGGGTTAATGCCCCGCAACCCCAACCTGCCcgccggcg                                                                                                                                                                                             >  2:481320/1‑239 (MQ=255)
                                                gtgTTAATGGGCTTTTGCGCCAGCGTTTTTTTACCCAACTTTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCcggcg                                                                                                                                                                                                <  2:834091/233‑1 (MQ=255)
                                                 ggTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAAGGCCACGTAGCACCATGCCGCCcggggggggg                                                                                                                                                                                           >  2:670592/1‑238 (MQ=255)
                                                      aTGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGGGATGGGGTGTATTCGGGTTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCCGAAACCAAACTCAAGGCTGCCGTTGGGGGTAATCCCCCGCCCCCCCATCCCGCcc‑gcggcggggaaaaa                                                                                                                                                                                     >  2:732167/1‑232 (MQ=255)
                                                      aTGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAAt                                                                                                                                                                                      >  1:12926/1‑239 (MQ=255)
                                                      aTGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGGGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCGGCCGTTCGGGTGACTGCCACTGAGACCCAGGCTCCCcgggggcggggacat                                                                                                                                                                                      >  2:708237/1‑233 (MQ=255)
                                                         ggCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAAttt                                                                                                                                                                                    >  1:946840/1‑238 (MQ=255)
                                                            aTTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGGGCAATTTgcgc                                                                                                                                                                                >  1:755704/1‑239 (MQ=255)
                                                                  gcCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGt                                                                                                                                                                          <  1:670592/239‑1 (MQ=255)
                                                                              tttcccAATTGTTTTTTGGTCCCTGGGCGCACGCCCACTATGTCAGCATTGGGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATc                                                                                                                                                             <  2:754917/236‑1 (MQ=255)
                                                                              ttACCCAACTTTTTTTTTGTCCCTGGCAGCCCGCCCATCAGTTCATCAGTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATc                                                                                                                                                             <  2:148768/239‑1 (MQ=255)
                                                                              ttACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATc                                                                                                                                                             <  1:184106/239‑1 (MQ=255)
                                                                              tgtcccattttttCTTTGGGCACTGGCGCCGCGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTGACGGTGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATc                                                                                                                                                             <  2:283908/228‑1 (MQ=255)
                                                                                 cccAACTTGTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTGAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCAcg                                                                                                                                                          <  1:821093/239‑1 (MQ=255)
                                                                                   catcTGTTCTTTGGTCACTGGCCACCCGCCCATCAGTTCATCATTGAGATAACCCTTTAACGGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCAcgcg                                                                                                                                                        <  2:23316/236‑1 (MQ=255)
                                                                                   cAACTGTTCTTTGGTCACTGGCAGCCCGCCCATCAGTTCATCAGTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCAcgcg                                                                                                                                                        <  1:267713/239‑1 (MQ=255)
                                                                                   cAACTGTTCTTTGGGCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAGCCTTTAACTGCGACTGGACATGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTGCATCGCTGCGCATACCAAACGCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCAcgcg                                                                                                                                                        <  1:1050652/239‑1 (MQ=255)
                                                                                        ggtCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGa                                                                                                                                                   <  1:679041/237‑1 (MQ=255)
                                                                                              tGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGc                                                                                                                                                                                                                        >  1:609422/1‑165 (MQ=255)
                                                                                                   aCTGGCGGCCCCCCCAGCGTTTCATGATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGAGGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTgcgc                                                                                                                                        <  2:10913/239‑1 (MQ=255)
                                                                                                   aCTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGt                                                                                                                                                                          <  2:500269/206‑1 (MQ=255)
                                                                                                   aCTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGt                                                                                                                                                                          >  1:500269/1‑206 (MQ=255)
                                                                                                   aCTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACCTAGCACCATGCTGCCCGGCGGCGGGGCAATTTACGCAAAGGTCAGCTT‑CACATCACCCGAAGGACCGTTCTgcc                                                                                                                                         >  2:712270/1‑237 (MQ=255)
                                                                                                           caCGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCCCGACTg                                                                                                                                >  1:88146/1‑239 (MQ=255)
                                                                                                             cGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGCCGGGGCAATTTGCCCAAAGGTCAGCTTCCACTCCAC‑CAAAAGCCCGTTCTGCGccccacaaggg                                                                                                                              >  2:921274/1‑230 (MQ=255)
                                                                                                                        tcatcaTTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGt                                                                                                                                                                                                                                    <  1:534639/127‑1 (MQ=255)
                                                                                                                        tcatcaTTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGt                                                                                                                                                                                                                                    >  2:534639/1‑127 (MQ=255)
                                                                                                                          atcatTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATg                                                                                                                                                                                                                         >  2:1082616/1‑136 (MQ=255)
                                                                                                                          atcatTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATg                                                                                                                                                                                                                         <  1:1082616/136‑1 (MQ=255)
                                                                                                                          atcatTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTAACCCCTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCAGGCGGCGGTGCACTTTGCGCAAAGGTCAGCTC‑CACATCACGCGAAGGACCGTTCTCCGCCACGACTGTCGCATCACC‑Gccc                                                                                                                 >  1:640584/1‑238 (MQ=255)
                                                                                                                               ttgaaTAAACCTTTAACTGCGGCTGGTAAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACg                                                                                                           <  2:1011500/235‑1 (MQ=255)
                                                                                                                                   aTAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCCCGACTGGCCGATCAGC‑GCCCTGCACGCCCg                                                                                                       >  1:568099/1‑239 (MQ=255)
                                                                                                                                   aTAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTCCCCGCCGGCGGGGCAATTTGCGCAAAGGTACACTT‑CACATCACACCAAGGACGCTTCTGCGCCACGACAGgcccgccaaacccccgcccccccg                                                                                                        >  2:448041/1‑217 (MQ=255)
                                                                                                                                      aaCCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGt                                                                                                                                                                          >  2:241606/1‑171 (MQ=255)
                                                                                                                                            tAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CCCATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCCGCCCCGGCCTCCCt                                                                                               >  1:158473/1‑238 (MQ=255)
                                                                                                                                            tAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CCCATCACCCGAAAGACCGTTCCGCGCCCCCACTGGCGCATCAGC‑CCCCTGCCCGCCCGGcccccct                                                                                               >  2:1048552/1‑232 (MQ=255)
                                                                                                                                                  cGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGt                                                                                                                                                                          <  2:616094/159‑1 (MQ=255)
                                                                                                                                                  cGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCCCATCACC‑GCCCTGCACGCCCGGCATCCCTTGcccc                                                                                         >  1:305777/1‑238 (MQ=255)
                                                                                                                                                   gACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTTCACGCCCGGCATCACTTGccccc                                                                                        >  1:438292/1‑237 (MQ=255)
                                                                                                                                                   gACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCCGCGGCGCAATTTGCGCAAAGGTCAGCTT‑CACATCCCCCGCAGCACCGTTCTGCTCCACGACTGGCGCATCAGC‑GCCCTGCGCGCCCCGCACACCTTGCCCCa                                                                                        >  2:179692/1‑238 (MQ=255)
                                                                                                                                                        gACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACg                                                                                                                                                                                                                    >  2:302514/1‑111 (MQ=255)
                                                                                                                                                              gagtacggtgggcgagggGGTTTATTCGAGGTTGAGCATCGCCTTGGTCACCCGTTCACCGCGGTGCATACCAAACTCACTGGTGCCGTTCGGGGTAATGGCACTTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGAt                                                                               <  2:785940/226‑1 (MQ=255)
                                                                                                                                                              cAGTAACGGTGGCGATGGGGTGTATTCGAGGTTGAGCATGCCCTTCGTCCCCACTTCATCGCTGCGCATACCAAACGCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGAt                                                                               <  1:366321/236‑1 (MQ=255)
                                                                                                                                                              cAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGGCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGAt                                                                               <  1:444644/236‑1 (MQ=255)
                                                                                                                                                                  aaCGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTTCCCCACCTCCGGAATTTGtt                                                                         >  1:403153/1‑238 (MQ=255)
                                                                                                                                                                       tggcgatggGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGTCGCATCAAC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGCTCCg                                                                      >  1:435582/1‑236 (MQ=255)
                                                                                                                                                                       tggcgatggGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCAGCGGATTTTGTTCCg                                                                      >  1:916315/1‑236 (MQ=255)
                                                                                                                                                                           gatggGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGGGCCACGACTGG‑GCCGCCGCGGCCCTGCACGCCCGGCATCCCTTGCCCCCCCTGCGGCTTTTGTTCCGCCTgggg                                                               >  2:1040877/1‑238 (MQ=255)
                                                                                                                                                                                    ttttCGGGGTTGGCCCTCCCCTTCGTCACCACCTCATCGCTTCGCATCCCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCCCCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCCCGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGc                                                      <  2:849124/237‑1 (MQ=255)
                                                                                                                                                                                    ttttCGGGGTTGGACCTCCCCTTCGGCCCCACTTCATCGTTGCGCATACCAAACGCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGc                                                      <  2:355027/237‑1 (MQ=255)
                                                                                                                                                                                     aTTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCGGCGCATACCAAACTCAAGGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGc                                                      <  1:1040877/238‑1 (MQ=255)
                                                                                                                                                                                         ggggttgggccgccccTTCGTCCCCACTTCACCGCGGCGCAAACCAAACCCAATGCTGGCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGACGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCt                                                 <  2:663068/225‑1 (MQ=255)
                                                                                                                                                                                         ggggttgacactccccTTCGTCCCCCCTTCATCGGTGCGCATCCCAAACCCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCt                                                 <  2:734848/225‑1 (MQ=255)
                                                                                                                                                                                         gAGGTTGAGCCTCGCCTTCGTCACCACTTCATCGCTGCGCATTCCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCt                                                 <  1:850568/239‑1 (MQ=255)
                                                                                                                                                                                         gAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCCCCATGCTGCCCGGCGGCGGTGCAATTTTCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCCCGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTCCCCCCCCTGCGGATTTTGTTCCGTCCGGGCGGCTCCAGCCGGc                                                  >  1:1007888/1‑238 (MQ=255)
                                                                                                                                                                                         gAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGCTTTTGTTCCGTCTGGGCGGCTACAGCTGGc                                                  >  1:456907/1‑238 (MQ=255)
                                                                                                                                                                                          aGGTTGGGCGTCGGCGTGGGCACCACTTCAGCGCTGCGCAGCCCAAACCCAATGCTGCCGTTCGGGTTAATGCCTCGTAGCACCATGCTGCCCGGCGGCGGTGCAATTGGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCt                                                 <  1:448469/238‑1 (MQ=255)
                                                                                                                                                                                           ggTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCACTTTGCGCAAAGCTCAGCTT‑CACATCACGCGCAGGACCGTGCTGGGCCACGACGGGCGCATCAGC‑GCCCTGCCCGCCCGGCATCCATTGCCCCCCGCGCGGATTTTTTTCCCTCTGagcgggccccgccggccc                                                >  2:39433/1‑225 (MQ=255)
                                                                                                                                                                                              tGAGCCACGTCTTCGTACCCACTTGATCGCTGCGCATACCAAACTCAATGCGGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCa                                            <  1:1011264/239‑1 (MQ=255)
                                                                                                                                                                                              tGAGCATCGCCTTCGTCACCCCTTCATCGCTTCGCATACCAAACCCAATGCTGCCGGTCGGGTTAATGCCACTCAGCCCCATTGTGGCCGGCGGCGGTGCAATTGGCGCAAAGGTCAGCTT‑CACATCACGCGGAGGCCCGTTCGTCGCCACGGCTGGCGCATCATC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTACGGATTTTGTTCCGTCTTGGCGGCTACAGCTGGTTCAGCa                                            <  2:82864/239‑1 (MQ=255)
                                                                                                                                                                                              tGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGGATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGGGGCAACACCGGGCTCAGCa                                            >  2:850568/1‑239 (MQ=255)
                                                                                                                                                                                               gAGCATCGCCTTCGTCACCACTTCATCGCTGCCGATCCCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCa                                            <  1:1043492/238‑1 (MQ=255)
                                                                                                                                                                                               gAGCAGCGCCTTCGTCCACAGTGCAACGGTGTGCATTCCAAACTCAATGCCGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCa                                            <  1:732167/238‑1 (MQ=255)
                                                                                                                                                                                                aGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCAtt                                          >  1:319439/1‑239 (MQ=255)
                                                                                                                                                                                                       ccTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCCTCACTTGCCCCACCCGCGGATTTTGTTCCGTCCGGGCGGCTACAGCCGGCCCAACATTGAACAgg                                   >  2:437649/1‑238 (MQ=255)
                                                                                                                                                                                                          tctccaccccTTCCACGCTGCCCAGACCAAACCAAATGCTGCCGTTCGGGTTAATGCCACGTAGCCCCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCCCGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCATTGATCAGTGgt                                <  2:783853/235‑1 (MQ=255)
                                                                                                                                                                                                            gTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCaa                                                                                                                                                                              >  1:18536/1‑97 (MQ=255)
                                                                                                                                                                                                             tCACCACTTCATCGGTGCGCAACCCAAACGTAATGCTGGCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGGGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCATTGATCAGTGgttgc                             <  1:921274/239‑1 (MQ=255)
                                                                                                                                                                                                                              cgcATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCCCGACTGGCGCATCAGC‑GCCCTGCACACCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGCCTGGGCGGCTCCACCCGGCCCCCCCTTGACCAGGGGTTGCGTTGCTggggcggcccg            >  2:559889/1‑230 (MQ=255)
                                                                                                                                                                                                                              cgcATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCCCCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGCCCCCGCATTGTTCCGTGGTTGCGTTGCTGGCGTGGccca            >  2:182293/1‑237 (MQ=255)
                                                                                                                                                                                                                              cgcATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCCGCGGATTTTGTTCCGTCTGGGCGGCTACCGCTGGCTCAGCATTGACACGGGGTTGCGGTTCTGGCGGCGccca            >  2:713058/1‑237 (MQ=255)
                                                                                                                                                                                                                                        ccTCAATGCTGCCGTTCGGGTTAATGCCCCGTACCACCCTGCTGCCCGGGGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCCCGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCATTGATCAGTGGTTGCGTTGCTGGCGTCGCCTGCGTCATGaaa  <  2:1033502/238‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTTGGCGGCTACCACTGGCTCCCCCTTGGTCAGTGGTTGCCTTGCTGGCGGCCCCtcgcgactgaaa  >  2:374859/1‑228 (MQ=255)
                                                                                                                                                                                                                                        aCTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCt                                                 <  2:599041/192‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
TGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGTAGCACCATGCTGCCCGGCGGCGGTGCAATTTGCGCAAAGGTCAGCTT‑CACATCACGCGAAGGACCGTTCTGCGCCACGACTGGCGCATCAGC‑GCCCTGCACGCCCGGCATCACTTGCCCCACCTGCGGATTTTGTTCCGTCTGGGCGGCTACAGCTGGCTCAGCATTGATCAGTGGTTGCGTTGCTGGCGTCGCCTGCGTCATGAAA  >  NC_000913/3692077‑3692547

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