Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,124,034 T→C 12.1% E141G (GAG→GGG)  cytR ← Anti‑activator for CytR‑CRP nucleoside utilization regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,124,0340TC12.1% 182.6 / 2.6 80E141G (GAG→GGG) cytRAnti‑activator for CytR‑CRP nucleoside utilization regulon
Reads supporting (aligned to +/- strand):  ref base T (39/26);  new base C (4/6);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTT  >  NC_000913/4123812‑4124269
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
aTGCCGCAGCGCCGCAGCGCCGGAACAGAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCaaa                                                                                                                                                                                                                                 <  1:553739/239‑1 (MQ=255)
    cGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCACCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACgg                                                                                                                                                                                                                              <  2:646750/238‑1 (MQ=255)
             gCAGGGGCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑Gaa                                                                                                                                                                                                                   <  2:818970/239‑1 (MQ=255)
                      gAACATAGCCTTGCAGGCCGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGcaac                                                                                                                                                                                                          <  2:525887/239‑1 (MQ=255)
                              ccTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGc                                                                                                                                                                                                   <  2:50017/238‑1 (MQ=255)
                                         gtattggctacgCGGCATCTTTTTGGGACCGGCTTTTCCGCCAATTCGTTTATGCCCTTGCTCATATAAATAATTTACTGGCTCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCt                                                                                                                                                                                       <  1:330375/228‑1 (MQ=255)
                                          ttgTGCCACAGCGGCTTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCAACAAAGGCGGCGGTCAGATTGGCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCt                                                                                                                                                                                       <  1:872892/236‑1 (MQ=255)
                                          tAGGGACACAGCGGCATCCCCTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCt                                                                                                                                                                                       <  2:522654/238‑1 (MQ=255)
                                                    gCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATTAACTGTAGGCAACTCCAGCTCCGGGGCAAATTCGTTCCCCATCACCATCGGCGGCGGATTAACCTGTTCCTCAATGCTGGC‑ACCAAACGCCAGCCTT‑GAACCCACAAACAAAAATCCACCAAATTCCTT‑GGGGATGTTc                                                                                                                                                                           >  2:838282/1‑239 (MQ=255)
                                                          tctcTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGGGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACCCGGTTCCCCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCCAAAACATGCCCTCAATTTGCTT‑GGGGTAGTTCAAATC‑g                                                                                                                                                                    >  2:656359/1‑239 (MQ=255)
                                                             ctTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATa                                                                                                                                                                 <  2:201099/239‑1 (MQ=255)
                                                                      cGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATTAACTGTAGGCAGCTCCCACTCCGCTGGAAAATCGTTCGCCCTCCCCATAGGCGGACGGTAACGCTGGTCGCCAATGTTGGCAATCAAACGGC‑GGCTT‑GGACCCCACGACAAAACGCCACAAATTTACCT‑GTTGTTGAACCAAtatgaaaagactttt                                                                                                                                                         >  1:488762/1‑225 (MQ=255)
                                                                      cGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑TAACCCAGCAACAACATTCCATCAATTTTCTT‑GGTGGTGATCAAATC‑GATAAAGGgtttt                                                                                                                                                        >  2:588660/1‑239 (MQ=255)
                                                                        gCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACAAGCCCTCAATTTTCTT‑GGGGGTGATCAAATC‑GATAAAAGtttttt                                                                                                                                                       >  2:495550/1‑238 (MQ=255)
                                                                        gCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGGGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACCCTGTTCCTCAATGCTGGC‑ATCAAACCGCCGCCTT‑GAAACCAGCAACCACAATCCATCAAATTTCTT‑TGTGGTGATCAAAACGGAAAAAAGtttttt                                                                                                                                                       >  1:604648/1‑239 (MQ=255)
                                                                        gCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGAAATGAACAGTAAGCAACTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGGATACCCTGTTCCTCAATGCTGGC‑CTCAAACGGCAACCTTGGAACCCACCACAAAAATCCCAT‑AATTTGCTT‑GGggatggaaaaaaagaaaaaaggttttt                                                                                                                                                        >  1:315095/1‑216 (MQ=255)
                                                                                            aTGCCCTTGCTCATATAAAGAATTTACTGCATCAAATGCGGCGGGCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTTCGGTGCAAATTCGTTCGCCATCCCCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑TATAAAGGTTTTTTCCTGCTGATTTTGATGCGcac                                                                                                                                  <  2:568175/239‑1 (MQ=255)
                                                                                             tGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGGGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGACTAAACC‑GTTAAAGGGTTTTTCCTGCGGATTTT‑ATgcggcaca                                                                                                                                >  2:313930/1‑234 (MQ=255)
                                                                                                 ctttctatttttaattatTTCCCGCAACAAGTGCGGTGGTCAGTTTGTCGATATGAACTTTAGGCAGCTCCAGCTGCTGTGCAAATTCGGTCGCCGTCGCCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCg                                                                                                                             <  1:733256/227‑1 (MQ=255)
                                                                                                   tGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCCGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCTTCAATTCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGGGATGATTCAAAC‑GATAAAGGGTTTTTCCTTCTGATTTTTGTTAGCACAGTCGcc                                                                                                                           >  2:460601/1‑239 (MQ=255)
                                                                                                   tGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGGTCCCCCATGCTTGC‑ATCAAACCGCACCCCT‑TGACCCCCCACACACATTCCACCAAATTTCTT‑GGGGTGGATCAAAAC‑CAAAAAGGATTTTTCCTTCTGATTTTTATGGGcacaagccca                                                                                                                           >  1:732073/1‑233 (MQ=255)
                                                                                                       cATATAAATAATTTCCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAAt                                                                                                                        <  2:796991/238‑1 (MQ=255)
                                                                                                       cATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATc                                                                                                                       <  2:419915/239‑1 (MQ=255)
                                                                                                       cATATAAATAATGTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAAt                                                                                                                        <  2:286225/238‑1 (MQ=255)
                                                                                                         tataAATAATTTACTGCATCAAATGCGGCGGTTAGATTGTCGATCTGAACTGTAGGCAACTTCAGCTCCGGTGCAAAATCGTTCGCCATCACCATCGGCGGCAGAGTACGCTGTTACTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGGGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAg                                                                                                                     <  2:867919/239‑1 (MQ=255)
                                                                                                          ataaataaTTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTTCGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAAACCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGGTTTTTCCTGCTGATTTTGATGGGGCCCGTCGCCAATCAGc                                                                                                                    >  2:104299/1‑239 (MQ=255)
                                                                                                          ataaataaTTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCCAATC‑GATAAAGGGTTTTTCCTGCCTATTTTTATTCGCACAGTCGCCAATCAg                                                                                                                     >  2:132582/1‑238 (MQ=255)
                                                                                                              ataatTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCCTCAATTTTCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATTCGCCAAGTCGCCAAACAGCAccc                                                                                                                >  2:319664/1‑238 (MQ=255)
                                                                                                                aatTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GTTAAAGGTTTTTTCCTGCTGATTTTTTTGTGCACAACTGCCAAACAccaccaaa                                                                                                              >  2:36620/1‑237 (MQ=255)
                                                                                                                  tttACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑TAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGGTTTTTCCTGCTGATTTTGATTCGCACAGTCGCCAATCAGCACCaaaaa                                                                                                            >  2:673109/1‑235 (MQ=255)
                                                                                                                  tttACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCCCAAGCGCCAATCAGCACCA‑Cat                                                                                                            >  2:313064/1‑238 (MQ=255)
                                                                                                                  tttACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGGGATGATCAAATC‑GATAAAGGGTTTTTCCTGCTGATTTTTATGCGCACAGTCGCCAATCAGCAcccggta                                                                                                            >  2:480559/1‑234 (MQ=255)
                                                                                                                  tttACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAAATGTAGGCAGCTCCCACTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTAAGCCGTTCCCCCATTCTGGC‑ACTCAACGGCAGCCTT‑GGAACCAGCAACAACATACCCTCAATTTTCTTGGGGGATG‑TCAAAAC‑GTGAAAGGGTTTTTCCGGCGGTTTTTGGTTCGCACAAGCCCCAAAAAGAACCA‑Gaca                                                                                                           >  1:776937/1‑237 (MQ=255)
                                                                                                                        gCATCAAATGCGGCGGTCAGGTTGGCGATATGAAATGTAGGCAGCTCCAGCGCCGGTGCAAATTCGTTCGCCCTCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAAACCAGCCACAACATGCCATCAATTTTCTT‑GGTGAGTATCAAATA‑GATAAAGGTTTTTTCCTGCTGATTTTGTTGCGCACAGCCGCCAATCAGCACCC‑Gaattc                                                                                                         >  2:579948/1‑231 (MQ=255)
                                                                                                                         cATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTTCTT‑GGTGATGATCCAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGTTGCGCCCAATCGCCAATAAGCACCA‑GAAATCCCTGa                                                                                                    >  2:36752/1‑239 (MQ=255)
                                                                                                                         cATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGGTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GAAATCCGTGa                                                                                                    >  2:816046/1‑239 (MQ=255)
                                                                                                                         cATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCCGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCCGCAACAACATGCCATCAATTTGCTT‑GGGGATGATCAAATC‑GATAAAGGGTTTTTCCTGCTGATTTTTGTGCGCACAAGCGCCAAACAACACCA‑GAAAACCCTGa                                                                                                    >  1:437116/1‑239 (MQ=255)
                                                                                                                            cAAATGCGGCGGTCAGATTGTCGATATGAAATGTAAGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATACCCATCGCCGGCAGAATACCCTGTTCCTCCATGCTGGC‑CTCCAACGGCAGCCTT‑TAACCCCACAACAACATGCCACCAATTTTCTT‑GGGGCTCGTCAAATC‑GATAAACGGTTTTTCCTCCTGATTTTGTTGCGCACAATAGCCAAACaaaacccaaaaacccggattt                                                                                                  >  1:857750/1‑219 (MQ=255)
                                                                                                                                   ggcggTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCt                                                                                                                                                                                                                                                  <  1:64030/92‑1 (MQ=255)
                                                                                                                                   ggcggTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCt                                                                                                                                                                                                                                                  >  2:64030/1‑92 (MQ=255)
                                                                                                                                    gcggTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTAAGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GGTAAAGGGTTTTTCCTTCTTATTTTTATGCGCACAGTCGCCAATAAGCCCCAGGATA‑CCGTGATTTTCCGCCGt                                                                                         >  2:71902/1‑239 (MQ=255)
                                                                                                                                     cggTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGt                                                                                         <  2:732073/238‑1 (MQ=255)
                                                                                                                                              tggcGATATGAACTGTAGGCAGCCCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATa                                                                               <  2:315095/236‑1 (MQ=255)
                                                                                                                                                   ataaGAACTGGAGGCCACTCCAGCTCCGGGTCAAATTCGTTCCCCATCACCATCGGCGGCAGATTCCGCTGTTCCTAAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAAGC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATAcc                                                                             <  2:91487/236‑1 (MQ=255)
                                                                                                                                                    ttggaattttgggCCGCCCCAGCCCCGGTGAAAATTCGTTCGCCATCACCATCGGCGGAAGATTACGCTTTTCCTCAATGGCGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGa                                                                         <  1:816046/228‑1 (MQ=255)
                                                                                                                                                      tGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCAATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGAACAGTTGCCAATCAACACCA‑GATATACGTGATTTGCAGCCGTAACTTCGATACCGCGaaa                                                                       >  2:351885/1‑238 (MQ=255)
                                                                                                                                                          cTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCCCAGTCGCCAATCAGCACCA‑GATATCCGTGATTTTCCCCCGTAACTTCGATACCGCGAAAAAtt                                                                   >  1:840644/1‑239 (MQ=255)
                                                                                                                                                            gTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCCCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCCTAACTTCGATACCGCGAATAAttt                                                                  >  2:742537/1‑238 (MQ=255)
                                                                                                                                                             tAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATCAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAAACAGCCCCA‑GATATCCGGGATTTGCCCCCGGAACTTCGAACCCGCGCATAAtt                                                                   >  1:91487/1‑236 (MQ=255)
                                                                                                                                                             gaGGCAGCCCCAGCTCCGGTACAAATTCCTTCCCCAACACCATCGGCGGTAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGGGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCg                                                                <  2:809763/238‑1 (MQ=255)
                                                                                                                                                               gcccgctccgtccccGGTGCAAATTTGTTCGCCATCCCCATCGGCGTCAGATTGAGCTGTTCCACAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCg                                                                <  1:167366/224‑1 (MQ=255)
                                                                                                                                                                  agctccagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGACTAATTTCGCTaaa                                                           >  2:792834/1‑239 (MQ=255)
                                                                                                                                                                    ctccagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAaaga                                                         >  2:180708/1‑239 (MQ=255)
                                                                                                                                                                    ctccagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGATTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGGGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAaaga                                                         >  2:469142/1‑239 (MQ=255)
                                                                                                                                                                     tccagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCCAAGTCGCCAATCAGCACCC‑GATATCCCTGATTTGCCGCCGTAACTTCGATACCCCGAAAAAATTCCCTaaaaaa                                                        >  1:267524/1‑239 (MQ=255)
                                                                                                                                                                       cagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAATAAgg                                                      >  2:27148/1‑239 (MQ=255)
                                                                                                                                                                       cagctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGAACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAgg                                                      >  2:663704/1‑239 (MQ=255)
                                                                                                                                                                          ctccGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTTCCGCCGTAACTTCGATACCCCGAAAAATTTCGCCAAAGAAGGGaa                                                   >  1:239557/1‑238 (MQ=255)
                                                                                                                                                                          ctccGGTGCAAATCCGTTCGCCATCACCACCGCGGGCAGATTGCTCTGTTCCGCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTG‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGAAGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGAt                                                   <  2:126516/239‑1 (MQ=255)
                                                                                                                                                                          ccccGGTGAAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGAt                                                   <  2:242513/237‑1 (MQ=255)
                                                                                                                                                                           tccGGTGCAAATTCGTTCGCCACCCCCATCGGCGGAAGATTACGCTTTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGAt                                                   <  2:367325/238‑1 (MQ=255)
                                                                                                                                                                           cccGGTTAAAATTCGTTCCCCATCCCCAAGGGTGGAAGGTTACCCTGTTCCCAAATGCTGGA‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCAACAATTTGTTT‑GGTGATGACCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGAt                                                   <  2:776937/237‑1 (MQ=255)
                                                                                                                                                                              ggTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACCTGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGACGCCGAAACTTCGCAACCGCGACTAATTTCGCTAAAGAAGGGATCGCa                                               >  2:598537/1‑239 (MQ=255)
                                                                                                                                                                              ggTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGCTTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCa                                               >  2:588318/1‑239 (MQ=255)
                                                                                                                                                                               gtgcatatttttctcctaTCCCCCTCGGCGGCAGTTTTCCCTTTTCCACAATGTTGGC‑AGAAAGCGGCAGCCTT‑GACCCCAGCAACAACATGCCAACAATTTGCTT‑GGTGATGGTCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAg                                              <  1:313930/222‑1 (MQ=255)
                                                                                                                                                                               gTGCAAATTCGTTCGCCATCACCATCGGCGGCAGCTTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGGGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATAGCAg                                              >  2:613571/1‑239 (MQ=255)
                                                                                                                                                                               gTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAgggctttc                                                >  2:167366/1‑232 (MQ=255)
                                                                                                                                                                                taaaaactcttccccccaccccACCGGCGGCAAATTCCGCTGTCCCTAAAGTCGGGC‑ATCAACCGGCACCCTT‑GAACCCAGCAACAACAGGCCATCAATTTGTTG‑GGTGATTATCAAACC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAg                                              <  1:132582/218‑1 (MQ=255)
                                                                                                                                                                                          ttcgCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATTAAATC‑GATAAAGGTTTTTTCCTGCTGCTTTTGATGGGCACAGTCGCCAACCAGCACCA‑GATATCCGTGATTTTCCGCCTGAAATACGATACCGCGAATAAATCTGCTAAAGAAGGGCCGGCAGAAATCCGGGAc                                   >  1:764157/1‑239 (MQ=255)
                                                                                                                                                                                           tcgCCATCCCCATCGGGGGCAGATTCCGTGGTTCCCCAATGCTGGC‑ATGAACCGGCAGCCTT‑GAACCCAGCAAAAGCATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACa                                  <  2:893594/239‑1 (MQ=255)
                                                                                                                                                                                           ccccactccccccggggggcAGATTCCCCGTTTCCCCAATGCTGGG‑AACAAACGGCAGCCTT‑GAACCAACCAACAACACGCCAACAATTTGCTT‑GGGGATAATAAAATA‑GATAAAGGTTTTTTCCGGCTGATTTTGATGCGCACAGTCCCCAATCAGCACCA‑TATACCCGTGATTGGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACa                                  <  1:838282/222‑1 (MQ=255)
                                                                                                                                                                                             cccATCCCCATCGGGGGAAGATTCCGGTGTTCCTCAATGCTGGC‑AACAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTGGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAAt                                <  2:485839/238‑1 (MQ=255)
                                                                                                                                                                                                     cATCGGCGGCAGATTACGGTGTTCCAAAATGCTGGC‑ATCAAACGGAAGCCTT‑GAACCCAGCAACAACAGGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGaa                        <  2:857750/239‑1 (MQ=255)
                                                                                                                                                                                                         ggcggcAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACACCATGCCATCAATTTGCTT‑GGTGTTGATCAAATC‑GATAAAGGCTTTTTCCTGCTGTTTTTGATGCGACCAGTCGCCACTCAGAGCCC‑AATATCCGTGATTTTCCGCCGTAACTTCCATAACGCGAATAATTTTCCTAAAGAAGGGATGGCAGATATCCGGGACAATACCCAAAATggg                    >  1:444662/1‑238 (MQ=255)
                                                                                                                                                                                                           cggcaatttccctttttCCTCAATGCTGGC‑ATAAAACGGCAGCATA‑GAACCCAGCAAACACATGCCATCAATTTGCTT‑GGTGATGATAAAATC‑GATAAAGTTTTTTTCCTGCTGAGTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGc                  <  1:852744/224‑1 (MQ=255)
                                                                                                                                                                                                           cggcAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAGTTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATa                                                                               >  1:325079/1‑178 (MQ=255)
                                                                                                                                                                                                           cggcAGATTACGCTGTTCCCAAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATa                                                                               <  2:325079/178‑1 (MQ=255)
                                                                                                                                                                                                           cggcAAATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGAAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTg                   <  2:658815/238‑1 (MQ=255)
                                                                                                                                                                                                               aGATTACGCTTTTCCCCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCggg               <  2:241709/238‑1 (MQ=255)
                                                                                                                                                                                                                    accgtGTTCCACAATGCTGGC‑ATCCAACGGCAGCTTT‑GAACCCAGCAAACACATGCCATCAATTTGCTT‑GGTGATTATCAAATC‑GATAAAGGTTTTTTCCGTCGGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATCTCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCAt         <  2:488762/235‑1 (MQ=255)
                                                                                                                                                                                                                    aCGCTTTTCCCCAATGGTGGC‑ATCAACCGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCAt         <  2:437116/239‑1 (MQ=255)
                                                                                                                                                                                                                           tCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCtt  >  1:859602/1‑239 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
ATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGC‑ATCAAACGGCAGCCTT‑GAACCCAGCAACAACATGCCATCAATTTGCTT‑GGTGATGATCAAATC‑GATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCA‑GATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTT  >  NC_000913/4123812‑4124269

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