Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913107,5530T.35.4% 211.7 / 111.4 82intergenic (+79/‑152)lpxC/secMUDP‑3‑O‑acyl‑N‑acetylglucosamine deacetylase/SecA translation regulator
Reads supporting (aligned to +/- strand):  ref base T (22/31);  new base . (17/12);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCACAATATTATTGGTGCATTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTGCGGCGAGTTTA  >  NC_000913/107309‑107791
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
tCACAATATTATTGGTGCATTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTACAGGTTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAAttt                                                                                                                                                                                                                                          <  2:262869/250‑1 (MQ=255)
  acaATATTATTGGTGCATTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAAtttt                                                                                                                                                                                                                                         >  1:708709/1‑250 (MQ=255)
     atatTATTGGTGCATTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAag                                                                                                                                                                                                                                    <  2:279074/251‑1 (MQ=255)
          attGTTGCATTTACCGCTTATAAATCCGGTCATGCAGTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACGGTCCTGGACTACTGAGAACTGCCGGTGGCGTTCAAAGCGCCTTCAGCTGTACGGGCGAAACTACATTGAGACTGTCGTATAAAAGTCGTCTGGCAAATCTGGCACTCCCACGTTCGACGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCt                                                                                                                                                                                                                                <  1:839578/251‑1 (MQ=255)
            tGGTGCATTTCCCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCta                                                                                                                                                                                                                               <  2:307552/250‑1 (MQ=255)
                  aTTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATaaaaaa                                                                                                                                                                                                                        <  2:688849/251‑1 (MQ=255)
                       ccGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGTCCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACGGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGc                                                                                                                                                                                                                  <  1:633484/251‑1 (MQ=255)
                           ttATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                                                                                                                                                                              >  2:155168/1‑251 (MQ=255)
                            tataAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGACGCATGGGAATATGTGACCTTCAAGGACGACGCAGCACTGCCGTTGGCCTTCAAAGCGACTTACGCTGTACTGGCATAACGACATTGATACTGTCGTATAAAATTCGACTGGCAAATCTGGCCCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GATTTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                                                                                                                                                                             >  2:379076/1‑251 (MQ=255)
                            tataAATCCGGTCATGCACTGAATAACAAAATGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCGCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                                                                                                                                                                              <  2:826199/250‑1 (MQ=255)
                                 aTCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                                                                                                                                                                         >  2:495706/1‑251 (MQ=255)
                                      gTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                                                                                                                    >  2:454736/1‑251 (MQ=255)
                                              cTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATaa                                                                                                                                                                                                                                     <  2:697299/210‑1 (MQ=255)
                                              cTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATaa                                                                                                                                                                                                                                     >  1:697299/1‑210 (MQ=255)
                                                gAATAACAAACTGCTGCACCCTGTCCTGGCGAAACAGGACGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGACGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTACGAACGCTGTTTTCTTAAGCACttt                                                                                                                                                                                          >  2:859034/1‑251 (MQ=255)
                                                    aacaacCTGCTGCAGGCTGTCCTGGCGACACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATAATGTCGTATAAAATTCGACTGGCCAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAACCGGTGCGAAc                                                                                                                                                                                                              >  2:77083/1‑227 (MQ=255)
                                                    aaCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAGGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCGTGCACTCTCTCCTGCCAGGTTAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                                                                                                                                                                              <  1:77083/227‑1 (MQ=255)
                                                     aCAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                                                                                                                                                                                     >  1:189466/1‑251 (MQ=255)
                                                        aaCTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGCATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTCAGCACTTTTCCGcac                                                                                                                                                                                   >  2:426856/1‑250 (MQ=255)
                                                        aaCTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                                                                                                                                                                                  >  2:283981/1‑219 (MQ=255)
                                                        aaCTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                                                                                                                                                                                  <  1:283981/219‑1 (MQ=255)
                                                            gctgcAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                                                                                                                                                                              <  1:495706/251‑1 (MQ=255)
                                                              ggcAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                                                                                                                                                                            <  2:708709/250‑1 (MQ=255)
                                                               gcAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                                                                                                                                                                            <  1:260592/250‑1 (MQ=255)
                                                                 aGGCTGTCCTGGCGCACCAGGCAGCCTGGGCATATGTGACCTTCCAGGACGACGCAGACCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                                                                                                                                                                         >  2:884049/1‑251 (MQ=255)
                                                                      gTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt                                                                                                                                                                      <  1:426252/249‑1 (MQ=255)
                                                                      gTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt                                                                                                                                                                     <  1:12905/249‑1 (MQ=255)
                                                                      gTCCTGGCGAAACAGGAAGCCGGGGAATATGTGACCTTCCAGGACGACGCATAACTGCCTTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc                                                                                                                                                                    <  2:32912/251‑1 (MQ=255)
                                                                           ggCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCTAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCt                                                                                                                                                               <  1:859034/251‑1 (MQ=255)
                                                                            ccgatcgaGGAAGCCGGGGCAGATGGGACCTGCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAGCTTCATTCGTGCtg                                                                                                                                                              <  2:309621/244‑1 (MQ=255)
                                                                                   aGGAAGCCTGGGCATATGTGACCTTCCAGGACGACGCCGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCACATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGCACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg                                                                                                                                                              >  1:309621/1‑244 (MQ=255)
                                                                                     gAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                                                                                                                                                                       >  1:882871/1‑200 (MQ=255)
                                                                                     gAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGGTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                                                                                                                                                                       <  2:882871/200‑1 (MQ=255)
                                                                                       aGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg                                                                                                                                                                 >  1:206059/1‑236 (MQ=255)
                                                                                       aGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg                                                                                                                                                                 <  2:206059/236‑1 (MQ=255)
                                                                                            gggAATATGTGACCTTCCAGGACGACGCCGAACTGCCGTTGGACTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCttt                                                                                                                                               >  1:801556/1‑250 (MQ=255)
                                                                                                atatGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACt                                                                                                                                                                               >  1:78005/1‑213 (MQ=255)
                                                                                                atatGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTGTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACt                                                                                                                                                                               <  2:78005/213‑1 (MQ=255)
                                                                                                agatGTGACCTTCAAGGACGACGCACAACGGCCGTTGGCCTTCAAAGCGACTTCAGCGGTACTGGCATAACGACATTTATACTGTCGGATAAAATTCGACTGGCAAATCTGGCACTCTTTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCATGCTTTAATGAt                                                                                                                                         <  2:909549/249‑1 (MQ=255)
                                                                                                   tgtgACCTTACAGGACGACGCCGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCCGGCTTTACTGATa                                                                                                                                        >  2:580419/1‑250 (MQ=255)
                                                                                                      gACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGAACAACTTATTTTCATTCGTGATGGGGACTGCAGGCTTTAATGATACGAt                                                                                                                                    >  2:944910/1‑251 (MQ=255)
                                                                                                       aCCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGAtt                                                                                                                                   <  1:208720/251‑1 (MQ=255)
                                                                                                          ttCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATTAGAGCTATTAAAAACCGGGCGAAAGCTGTTTTCTTTAGCCCTTTTCCGCACCACCTATCTTTCTTCCTGGTGTGGACTGCCGGCCTTAATGAttagcattgtg                                                                                                                               >  1:54919/1‑241 (MQ=255)
                                                                                                            ccAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgc                                                                                                                              <  1:454736/251‑1 (MQ=255)
                                                                                                            ccAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATACGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgc                                                                                                                              >  1:643072/1‑251 (MQ=255)
                                                                                                              aGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgcgc                                                                                                                            <  1:524754/251‑1 (MQ=255)
                                                                                                                gacgacGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTa                                                                                                                          >  1:355354/1‑250 (MQ=255)
                                                                                                                     cgCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATGCGACTGGCAAATCTGGCGCTCTCTCCGGCCAGGTGAACCATTCGGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACGTTTCCGCACAACTTAGCTTCATTCGTGCTGTGTACTGCAGGCCTTAATGATAAGATGTGTGCGCTAAATAc                                                                                                                     <  2:89122/251‑1 (MQ=255)
                                                                                                                          aaCTGCCGTTGGCCGTCAAAGCGCCTTTAGCTGTACTGAACTAACGCCATTTATACTGTCGTATAAAATTCGACAGGCACATCTGGCAGTCTCTCCGGCCAGGAGAACCAGTCGTTTTTTTT‑GAATTTTATACGAGCTATAAACAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCCCAACTTATCTTCATTCGTGCTGTGGTCTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTg                                                                                                                >  1:592048/1‑250 (MQ=255)
                                                                                                                           acggCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGCCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTTTTTTCTTAAGCACGTTTCCGCACAACTTATCGTCATTCGTGCTGTTGACTGc                                                                                                                                                      <  2:107529/209‑1 (MQ=255)
                                                                                                                           aCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGCATCAACTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGCGCTATAAAAAACGGTGCGAACGCTGTTTTCTTCAGCACTTTTCCGCACAACTTATCTTCCTTCGTGCTGTGGACTGc                                                                                                                                                      >  1:107529/1‑212 (MQ=255)
                                                                                                                              gCCGTTGGCCTTCAAAGCGACTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAat                                                                                                             >  1:372012/1‑250 (MQ=255)
                                                                                                                                cGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATg                                                                                                          >  2:202668/1‑250 (MQ=255)
                                                                                                                                           aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCa                                                                                              <  1:202668/251‑1 (MQ=255)
                                                                                                                                            aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGCCTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGCATATGTTCGTGATGGCa                                                                                              >  2:36550/1‑250 (MQ=255)
                                                                                                                                                    tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTTGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACg                                                                                        <  1:340465/249‑1 (MQ=255)
                                                                                                                                                     cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATACGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGa                                                                                     <  2:643072/251‑1 (MQ=255)
                                                                                                                                                     cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTg                                                                                      >  2:935590/1‑249 (MQ=255)
                                                                                                                                                         tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTgg                                                                                 >  1:534640/1‑251 (MQ=255)
                                                                                                                                                          gtccTGGCCTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGTTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATGCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGGAATAACGTGAGTgg                                                                                 <  2:355354/246‑1 (MQ=255)
                                                                                                                                                            aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTGTTCTTAAGCACTTTTCCGCACAACTTATCGTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCTCTAAATGCGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGaa                                                                               <  1:637566/250‑1 (MQ=255)
                                                                                                                                                                        aCATTTATACTGTCTTATAAACTTCGACTTGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTGTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATGCGTGCTGTGGACTGCAGGCTTTAATGAGAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGTAATACTGACTCGCTg                                                                  <  2:798021/251‑1 (MQ=255)
                                                                                                                                                                                     cGTATAAAATTCGACTTGCAATTCTTGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACGTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAATATTTGTGCGCTAAATACGTTTGAATATGAACGGGATGGCAATAACGTGAGTGGAATTCTGACGCGCTGGCGACAGTTTGGt                                                     <  2:753548/251‑1 (MQ=255)
                                                                                                                                                                                      gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGGTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGAGGCGCTGGCGACAGTTTGGTa                                                    <  1:371667/251‑1 (MQ=255)
                                                                                                                                                                                          aaaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCGTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAAc                                                 <  1:580419/250‑1 (MQ=255)
                                                                                                                                                                                              ttCGACTGGCAAATCTGGCCCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACt                                           >  2:930429/1‑251 (MQ=255)
                                                                                                                                                                                                 gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCGACTTc                                         <  1:297765/251‑1 (MQ=255)
                                                                                                                                                                                                   cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTg                                       >  1:176890/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCAAATCTGGGACTCTCTCCGGCCAGGTGCACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGcc                                    >  2:410083/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCAAATCTGGCCCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCGATAAAAAACGGTTCGAACGCTGTGTTCTGAAGCACTTTTCCGCTGAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTGTAATTATAAGAGTTGTGCGCTAAAATCGATTGAATATTATCGGGATGGCAATAACCTGAGTGGGATACTGACGCGCGGGCGACAGTTGGGTGAACCCTACTCCTGGcc                                    <  2:592048/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATACAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGc                                     >  2:627357/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGc                                  >  2:506640/1‑250 (MQ=255)
                                                                                                                                                                                                                 aCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTGTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCttat                         <  1:884049/251‑1 (MQ=255)
                                                                                                                                                                                                                       tcCGGCCAGGTGAACCAGTAGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAggg                    >  2:941378/1‑250 (MQ=255)
                                                                                                                                                                                                                        ccGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAggg                    <  2:534640/249‑1 (MQ=255)
                                                                                                                                                                                                                        ccGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGGTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAggg                    <  1:434838/249‑1 (MQ=255)
                                                                                                                                                                                                                          ggCCAGGAGAACCAGTCGTTTATTTT‑GAAATTTATAAGAGCGATGAAAAACGGTGCGAAAGCTGTTTTCTTAAGGACTTTTCCGCACCACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAGGATAAGAATTGTGCGCTAATTACGTTTGAATATTGTCGGGTCGGAACTACCGTGAGTGGAGTACTGGCGGGCTGGCGAGAGTTTGGTAAACGCTACTAGTGCCCGCATCTCTTATTAGGGATgc                >  2:340360/1‑249 (MQ=255)
                                                                                                                                                                                                                          gcccAGGTGAACCAGCCGTTTTTTTT‑GGATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGt               >  2:340333/3‑251 (MQ=255)
                                                                                                                                                                                                                               ggTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTgcg           <  1:627357/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAAGCTGTTTTCTTAAGCACTTTTCCGCACAAATTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAAAGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTgcg           >  1:655642/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTGAACCAGTCGTTTTTTTT‑GAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTgcggc         >  1:99378/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTGCGGCGAGTTTa  <  1:941378/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
TCACAATATTATTGGTGCATTTACCGCTTATAAATCCGGTCATGCACTGAATAACAAACTGCTGCAGGCTGTCCTGGCGAAACAGGAAGCCTGGGAATATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCTGGCGACAGTTTGGTAAACGCTACTTCTGGCCGCATCTCTTATTAGGGATGGTTGCGGCGAGTTTA  >  NC_000913/107309‑107791

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