Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1710AG57.1% 36.9 / 98.0 77E41E (GAA→GAG) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (19/14);  new base G (27/17);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAG  >  NC_000913/1428927‑1429414
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCg                                                                                                                                                                                                                                               >  1:317236/1‑251 (MQ=255)
 ggCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCgg                                                                                                                                                                                                                                              >  1:346426/1‑251 (MQ=255)
    agttgCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAGAAATTCTGATTAGCCAGGTAACCCAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATAGGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATg                                                                                                                                                                                                                                           <  2:185104/247‑1 (MQ=255)
     gTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCAACACGGTG‑GATGG‑GAACACGGTGGGCTCAGAGAATCCGGATgc                                                                                                                                                                                                                                          >  1:495757/1‑250 (MQ=255)
                cggaCTGACGGATGCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAATCCGGGCGtt                                                                                                                                                                                                                               <  2:460696/248‑1 (MQ=255)
                 gtaCTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGtt                                                                                                                                                                                                                               >  1:175382/3‑250 (MQ=255)
                      gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTGCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGCACCACTGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATTAATCCGGTCGTTACAGc                                                                                                                                                                                                                          <  2:378134/250‑1 (MQ=255)
                         ggattCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATgg                                                                                                                                                                                                                      >  1:64105/4‑251 (MQ=255)
                                 cGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGAGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGa                                                                                                                                                                                                               <  1:139373/250‑1 (MQ=255)
                                          cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAg                                                                                                                                                                                                     <  2:171667/251‑1 (MQ=255)
                                              tataAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAc                                                                                                                                                                                                  <  1:110908/250‑1 (MQ=255)
                                                    ttCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGCCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGCACGGCCGGTACAgtgtcc                                                                                                                                                                                           >  2:315031/1‑246 (MQ=255)
                                                        gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTACCAGCACCACGGAG‑G‑TGGTGAACACGGTGGGCTAAGAGAATCCGGATGAAGACGGGCGTTACAGCATGGATGGGGAGGACGGTCAGTACAgtggcaccct                                                                                                                                                                                       >  2:349461/1‑242 (MQ=255)
                                                          ttAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGt                                                                                                                                                                                     >  2:125431/1‑251 (MQ=255)
                                                              ccgggTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCGTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAGCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTTAGTACGGTCAGTACAGCGCTATTCTGTtggt                                                                                                                                                                                 <  2:481989/248‑1 (MQ=255)
                                                               cAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggtt                                                                                                                                                                                >  2:163313/1‑235 (MQ=255)
                                                                aGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttg                                                                                                                                                                               <  1:163313/251‑18 (MQ=255)
                                                                    aaCACAGTGTGATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAgtgtcatcctgcaggttgacg                                                                                                                                                                            <  2:62126/250‑21 (MQ=255)
                                                                       acaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggttt                                                                                                                                                                        >  1:448740/1‑227 (MQ=255)
                                                                                         gccgGTTCAGGCGGGTTTTTTTGGGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAACCGGGACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGGGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGATTACGGTCATTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCgg                                                                                                                                                      <  2:447809/251‑1 (MQ=255)
                                                                                             gtccAGGCGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCGGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc                                                                                                                                                  <  1:537475/248‑1 (MQ=255)
                                                                                               gcgggtgggCTTTTGTGTGGGGGGAATAGGTCAGGAAAGATTTCAGTTGTACTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCGGCTGAAAGAAAAACGAAACAGCCCCA‑GGTGTG‑GGGTGAACGCGCTGGCTTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAGGGTCATTTCAGCGGGATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt                                                                                                                                                <  1:34094/245‑1 (MQ=255)
                                                                                                  ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                                                                             >  1:320049/1‑251 (MQ=255)
                                                                                                     gggCTTTTTTGTGGGGTGAAGATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGAAAAACGTAACAACACCACGGTG‑G‑TGGTTAACACGCTGGCCTCGGAAAATCCGGATGAACCCGGGCGTGACAGCATGGACGTTGAGTAGGGTAAGTACAGCGTTATTCTGTTGGTGGAAGTATTCCCGCCGTCACATGCCGGGACCATTACCgt                                                                                                                                          <  1:369497/251‑1 (MQ=255)
                                                                                                      ggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCgtg                                                                                                                                         >  2:744/1‑251 (MQ=255)
                                                                                                       gCTTTTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCCCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑AGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTCCCgtgt                                                                                                                                        >  2:284960/1‑251 (MQ=255)
                                                                                                            tttGTGGGGTGAATATGGCAGTAAATATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCTCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGTGCTCATAGAATCCGGATGAAGCCGGGCGGTACAGCATGGATGTGTTGTACGGTCAGTACAGTGTCATCCTGCGGGGTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                                                                                                                    <  1:43405/250‑1 (MQ=255)
                                                                                                             ttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGaa                                                                                                                                   <  1:366256/250‑1 (MQ=255)
                                                                                                                tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTCCTGAAAGCCGGAACAGGAACACAGGCACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCAAGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACGGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGAtt                                                                                                                               >  2:561195/1‑251 (MQ=255)
                                                                                                                gggggTGAATAGGGCAGTAAAGATTTCAGGTGGACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGATGAAAGCAAAACGTAACAGCACCGCGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAAACCGGTTGAAGCCGGGCGTTACAGCATGGAGGTGGAGTACGTTCAGTACAGCGTTATTTTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGAtg                                                                                                                               <  2:182601/250‑2 (MQ=255)
                                                                                                                  gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACCGGAAAACCGGTAAAGAACTGCACCATTCAGCTGCAAGCCAGACGTAACAGCCCCACGGTG‑A‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGCGTACGGCAAGTACAGTGTCATCCTGCCGGCTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTtc                                                                                                                              >  2:493375/1‑250 (MQ=255)
                                                                                                                  gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAACATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGCTGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATtcc                                                                                                                             >  2:502073/1‑250 (MQ=255)
                                                                                                                    gTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAACACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGCAGATtctc                                                                                                                            >  1:26086/1‑250 (MQ=255)
                                                                                                                     tGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGACAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                                                                                                                           >  1:502116/1‑247 (MQ=255)
                                                                                                                         aatGGCTGTAAAGATTGCAGGTTAACTGAAAGCCGGCACAGGAAAACCGGTACAGAACTGCAATACCCAGCGGAAAGCAAAACTGAAAAGAACCACGGTG‑G‑GGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGCGTTACAGCATTGACGTGGAGTCGGGTCAGTCCAGCGTTATTCTTTTGGTGAAATTATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATGCTCAACCCg                                                                                                                      <  2:510030/250‑1 (MQ=255)
                                                                                                                              cAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAGGAAATCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTG‑G‑TGGTGACCCCGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                                                                 <  1:137197/251‑1 (MQ=255)
                                                                                                                              cAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggacg                                                                                                                 <  2:198451/251‑8 (MQ=255)
                                                                                                                                 tAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTg                                                                                                              >  2:156236/1‑251 (MQ=255)
                                                                                                                                  aaaGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCAGCGCACGCCGGGACCATCACCGGGTATGAAGATTCACAAccgtggacgctga                                                                                                             <  2:88738/251‑12 (MQ=255)
                                                                                                                                     gATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATg                                                                                                          >  1:180217/1‑251 (MQ=255)
                                                                                                                                         tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                                                                                      >  1:517741/1‑251 (MQ=255)
                                                                                                                                          caggagtccTGAAAGACGGCACAGACAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCCGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAtttt                                                                                                     >  2:258665/9‑251 (MQ=255)
                                                                                                                                                 ccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCtgtg                                                                                              >  1:138193/2‑251 (MQ=255)
                                                                                                                                                    gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                                                                                             <  1:175182/249‑1 (MQ=255)
                                                                                                                                                    gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa                                                                                           <  2:446098/251‑1 (MQ=255)
                                                                                                                                                    gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGTTGAACACGCTGGCCTCATAAAATCCGGAGGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGTTCAGTACGGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc                                                                                            <  2:341157/250‑1 (MQ=255)
                                                                                                                                                       aGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATg                                                                                         >  2:161128/1‑250 (MQ=255)
                                                                                                                                                          cGGCACAGGCAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGCTTTTCTCGGTGCCATGACgg                                                                                     >  1:454630/1‑251 (MQ=255)
                                                                                                                                                          cGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGCCgg                                                                                     >  1:220903/1‑251 (MQ=255)
                                                                                                                                                           ggCAAAGGAAAACCGGTACAGAACTGCACAATCCAGATGAAAGAAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAAAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACAGGCAGGGACCATTAACGTGTATGAAGATTCTCAACCCGGTAAGCTGAATGATTTTAACGGTGCCATGACgg                                                                                     >  1:279228/1‑250 (MQ=255)
                                                                                                                                                             cacaGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggagg                                                                                  >  1:19581/1‑251 (MQ=255)
                                                                                                                                                                        cGGTAAAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTG‑G‑TGGTGAACACTTGGGCCTCAGAAAAGCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCGGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGt                                                                       <  2:320049/251‑1 (MQ=255)
                                                                                                                                                                         ggTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggc                                                                      >  2:166983/1‑249 (MQ=255)
                                                                                                                                                                                  cggccaaaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGTGCTTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCAc                                                             <  1:200182/245‑1 (MQ=255)
                                                                                                                                                                                  cTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgc                                                             >  1:56678/1‑248 (MQ=255)
                                                                                                                                                                                     cacaATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGCTACGCTGAATGATTTTATCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgc                                                          >  2:4505/1‑251 (MQ=255)
                                                                                                                                                                                          tCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGt                                                     >  1:191609/1‑251 (MQ=255)
                                                                                                                                                                                           ccAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtt                                                    >  1:188442/1‑251 (MQ=255)
                                                                                                                                                                                                   aagccagaCGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTg                                            <  1:177645/244‑1 (MQ=255)
                                                                                                                                                                                                   aaGCAAAACGTAACAGCACCACGGTG‑G‑GGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCTCTGCGCCGTTTTGCACTg                                            >  2:456390/1‑251 (MQ=255)
                                                                                                                                                                                                      ccagaCGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGGACTGAtg                                         >  2:364392/5‑251 (MQ=255)
                                                                                                                                                                                                       aaaaCGTAACAGCACCACGGTG‑G‑TGGTGCACACGCTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGAtgg                                        >  2:470524/1‑251 (MQ=255)
                                                                                                                                                                                                           cGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg                                                                   <  2:392506/220‑1 (MQ=255)
                                                                                                                                                                                                           cGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg                                                                   >  1:392506/1‑220 (MQ=255)
                                                                                                                                                                                                            gTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGaa                                   <  1:156236/251‑1 (MQ=255)
                                                                                                                                                                                                            gTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                    >  1:281099/1‑250 (MQ=255)
                                                                                                                                                                                                                    accacGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGACGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGAACGCTGCATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGCGGCACTGCGCCGTTTTGAACTGATGGTGGAAGGGGTGgc                           >  2:183974/1‑251 (MQ=255)
                                                                                                                                                                                                                    accaaGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAGTGCTTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGCGGCCCTTCTCCGTTTTGAACTGATGGTGGAAGAGGTGgc                           >  2:320214/1‑251 (MQ=255)
                                                                                                                                                                                                                     ccacGGGG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGgcg                          <  2:223634/251‑1 (MQ=255)
                                                                                                                                                                                                                      cacGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgcgc                         >  1:28404/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACGGATGGTGGAATAGGTGGCGCGt                       <  2:56678/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggtg‑g‑tggtGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGt                       >  2:509459/1‑250 (MQ=255)
                                                                                                                                                                                                                          gtg‑g‑tggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACACCCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGCGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTa                      >  1:511099/1‑250 (MQ=255)
                                                                                                                                                                                                                          gtg‑g‑tggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTa                      >  1:46331/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cacGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgtgg          <  2:448740/251‑4 (MQ=255)
                                                                                                                                                                                                                                       cacGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgtgg          >  2:562228/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  cAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAg  >  1:139319/1‑248 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAG  >  NC_000913/1428927‑1429414

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