Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,890,626 G→A L425F (CTC→TTC)  cysJ ← sulfite reductase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6260GA96.8% 276.8 / ‑4.9 93L425F (CTC→TTC) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (51/39);  total (54/39)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGC  >  NC_000913/2890380‑2890870
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCgaagaa                                                                                                                                                                                                                                                   <  1:214545/250‑1 (MQ=255)
tGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGATATCGTTATGTTCGATACATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCGGCGAgaa                                                                                                                                                                                                                                                    >  2:1001/1‑247 (MQ=255)
    gTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCACTACCGGTGCCTGGGCCAATCATACTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTgg                                                                                                                                                                                                                                              >  1:38088/1‑251 (MQ=255)
           tGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGCCTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGGGGcacc                                                                                                                                                                                                                                       >  2:22036/1‑251 (MQ=255)
           tGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGcacc                                                                                                                                                                                                                                       >  1:507994/1‑251 (MQ=255)
            gCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGcac                                                                                                                                                                                                                                        <  1:365553/249‑1 (MQ=255)
                cTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGCCGCAATACCGGCGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTCTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCgg                                                                                                                                                                                                                                  >  1:391694/1‑251 (MQ=255)
                  ccACCTGGTACAGGAAGTCTTCCGTAAAGGGCGGATTCCCAAAGAACCGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCa                                                                                                                                                                                                                                <  2:238071/251‑1 (MQ=255)
                   cacatGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTCTGCCTGGAGCTTCGTCGGCGGAGCGTTGCTGCATAGGGGCGCGGCACGGCGCAATGCAGGTGACTGGGCCAATCATAATCACCCGGGTTACTGGATTGTCTGGCGGGAGGGAAATATCGTTATGTTTGCTAACATCACGGGCTTAGCCCTCTTCTTCCACGAGGCTAGCGAGAGAGCTCGGGGCACGACCGGCAc                                                                                                                                                                                                                               >  2:48369/5‑251 (MQ=255)
                        ggaaCAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGCCGGAAGTTATCGTTATGTTCGATAAATACGCGGCATTCGCCCTCTTCTTCCACGCGGTCAGCGCAGAAGCTGGAGGCACCACCCGCACGGgcc                                                                                                                                                                                                                          >  2:195391/4‑250 (MQ=255)
                        ggTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAATACCAGCCAGTTTGTCCCGGGCGCTTCGGCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTGCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGgcg                                                                                                                                                                                                                          <  2:476863/251‑1 (MQ=255)
                        ggTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACggg                                                                                                                                                                                                                            <  1:190782/249‑1 (MQ=255)
                          tACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGgcacg                                                                                                                                                                                                                        >  2:120235/1‑248 (MQ=255)
                          tACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGGGGCACCCCCGGCCCGGGcgcg                                                                                                                                                                                                                        >  2:243429/1‑251 (MQ=255)
                           aCAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCgg                                                                                                                                                                                                                       <  1:302280/251‑1 (MQ=255)
                               gAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCtt                                                                                                                                                                                                                   <  2:432933/251‑1 (MQ=255)
                                 cccctacacgacgcccttccgacctaaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTcc                                                                                                                                                                                                                 <  2:476421/226‑1 (MQ=255)
                                  gTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTcc                                                                                                                                                                                                                 >  1:451096/1‑250 (MQ=255)
                                   tCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAAGAAGCTGGAGGCACCCCCGGCACGGGCGCGGCCTTccgc                                                                                                                                                                                                               >  1:657442/1‑249 (MQ=255)
                                      tCCGTAAATTGCGGTGTACACAAGAACAGCCATTTTTTCCCTGGCGCTGCGTCGGTGGCTCGTTGCTGCATAAAGGCGCGGAACGGCGCAAGACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTGCTTCCACGCGGGCAGCGAAGAAGCTGGAGGCAACACCGGCACGGGCGCGGCCTTTCAcgt                                                                                                                                                                                                             <  2:115540/250‑1 (MQ=255)
                                      cccggaaaGTGCGAATTACCAAATAACAGCCAGTTTGTGCCCGGCGCTTCGGCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTGCCCCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTTCAcgtc                                                                                                                                                                                                            <  2:596520/246‑1 (MQ=255)
                                        cGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                                                                                                                                          >  1:397513/1‑251 (MQ=255)
                                         gTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                                                                                                                                          >  1:660136/1‑250 (MQ=255)
                                         gTAAAGGGGGGATAACCAAAGAAAAGCCAGTTTTTGCCTGGCGCGTCGTTGGCTGCGCGGTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCTTAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTGCCACGCGGGCAGCTAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTa                                                                                                                                                                                                         <  2:680062/251‑1 (MQ=255)
                                           aaaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                                                                                                                                       >  1:150232/1‑251 (MQ=255)
                                           aaaGTGCAGATTACCAAAAGGCAGCCAGTTTTTACCTGGCGCTAAGTCGGCGGCGCGTTGCTGAATAAAGGCGCGGAAAGGCGCAGGACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTCGCAGGAGGAAGTTATCGTTATGTTCGATAAATACGCGGACATCGACCTCCTCTTCCCAGCGGTCAGCGAAGCAGCTGGACGCGCCACCGGCACGGGCGCGGCCATCCAccgcgaac                                                                                                                                                                                                        >  2:95239/1‑243 (MQ=255)
                                                gCGGATTACCAAAGAACAGCCAGGTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGGTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGTACGGGCGCGGCCTTCCACGTCGTAACGcacc                                                                                                                                                                                                  <  2:663365/251‑1 (MQ=255)
                                                 cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGCTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACAGGACCGGGCGCGGCCTTCCACGTCGTAACGcccc                                                                                                                                                                                                  >  2:33588/1‑247 (MQ=255)
                                                      tACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacc                                                                                                                                                                                             <  1:1001/250‑1 (MQ=255)
                                                          aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTcc                                                                                                                                                                                                                 >  1:476421/1‑226 (MQ=255)
                                                            agcccAGCCAGTTTTTACCTGTGTCGTCGTCGGCTGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCTGAAGTTATCGTTATGTTCGATAAATACGCTGACTGCGCCCTCTTCGTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACTCACCACACCAACGGTg                                                                                                                                                                                      <  2:16049/247‑1 (MQ=255)
                                                              aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                                                                                                                                    >  1:457382/1‑251 (MQ=255)
                                                               aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACg                                                                                                                                                                                   >  1:394773/1‑251 (MQ=255)
                                                                   ccAGGTTTTACTTGGTGTTTCGGCGGCGGGGCGTTGCTGCATAAAGGCGCGGAACGGTGCAATACCGGTGCCTGGGCCAATCATAATCGCCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAGCGTTATGTGCGACAAATACGCGGACTGCGCCCTCTTCTTCCCCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGTGCGCGGCCTTCCACGTCGTAACGCACCACACCAACgtggatgtgta                                                                                                                                                                               <  2:150232/251‑11 (MQ=255)
                                                                           tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTCTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                                                                                                                                       >  1:405106/1‑251 (MQ=255)
                                                                               tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGGGTACTTCGTTCTCg                                                                                                                                                                    >  1:50745/1‑250 (MQ=255)
                                                                                   gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGc                                                                                                                                                                                                                                                 <  2:414553/169‑1 (MQ=255)
                                                                                   gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGc                                                                                                                                                                                                                                                 >  1:414553/1‑169 (MQ=255)
                                                                                   gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTc                                                                                                                                                               <  1:290668/251‑1 (MQ=255)
                                                                                                ggcgtggttgcATAAAGGCGCGGAACGGCGCAAATCCTGTGCCTTGGCTAATCATAAGCACCGGGGTTGCTGGATTTGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGAGGACTTCGCCCTCTTCTTCCACGCGGGCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCTGCCTTCCACGGCGTAACGCACCACACCAACGTTGACTTGTACTTCGTTCTCGACTTCCGCATGCgaggag                                                                                                                                                  <  2:569596/243‑1 (MQ=255)
                                                                                                   gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCg                                                                                                                                               >  1:344568/1‑251 (MQ=255)
                                                                                                   gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCg                                                                                                                                               >  1:133498/1‑251 (MQ=255)
                                                                                                           aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAAt                                                                                                                                        >  1:415174/1‑250 (MQ=255)
                                                                                                           aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATa                                                                                                                                       <  1:22036/251‑1 (MQ=255)
                                                                                                             aaaGGCTCGGAACGGCGCAATACCGGTGCCTGGGCTAATCATAATCACCGGTGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                                                                                                                     <  1:32921/251‑1 (MQ=255)
                                                                                                             aaaGGCTCGGAACGGCGAAATACCGTTGCCGGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAGCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGAACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                                                                                                                                      <  2:527205/250‑1 (MQ=255)
                                                                                                                 gcgcGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                                                                                                                                 <  1:381209/251‑1 (MQ=255)
                                                                                                                  cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGACTACAGAcgc                                                                                                                                >  1:387585/1‑251 (MQ=255)
                                                                                                                     ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                                                                                                                             >  1:286115/1‑251 (MQ=255)
                                                                                                                     ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTAGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACGTCCGCCTGCGCGGAGGCGATGGAACACAgcccaggc                                                                                                                             >  2:181532/1‑244 (MQ=255)
                                                                                                                          ggCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTc                                                                                                                          >  1:461408/1‑249 (MQ=255)
                                                                                                                           gcgcAATACCGGTGCCTGGGCCAATCAAAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGc                                                                                                                       <  2:451096/251‑1 (MQ=255)
                                                                                                                               aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgc                                                                                                                   <  1:87692/251‑1 (MQ=255)
                                                                                                                                   ccGGTGCCTGGGCCAATCATAATCTCCGGGGTTTCTGGATTGGCGGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCTCCGGCACGGGCGCGGTCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgca                                                                                                                <  1:608714/250‑1 (MQ=255)
                                                                                                                                        gCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGCGGCGATGGAATACAGCCGCGGCGTCAGCGGGCgcagc                                                                                                              >  1:654578/1‑247 (MQ=255)
                                                                                                                                         ccTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGatt                                                                                                         >  1:545474/1‑251 (MQ=255)
                                                                                                                                                cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGATATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGCTGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAgc                                                                                                  >  2:15428/1‑251 (MQ=255)
                                                                                                                                                 aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCCGCAGATTAATTAgcg                                                                                                 >  1:487892/1‑251 (MQ=255)
                                                                                                                                                  aTCAGAATCACCGGGGTTTCTGGATTGGCTGGCAGGAGGAAGTTACCGTTATGTTCGATAAATCCGCGGACTTCGCCCTCTTCATCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAgcgc                                                                                                <  2:185403/251‑1 (MQ=255)
                                                                                                                                                        aTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGc                                                                                          <  1:212474/251‑1 (MQ=255)
                                                                                                                                                               gggTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACAACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGc                                                                                   >  1:432318/1‑251 (MQ=255)
                                                                                                                                                                 gTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCGTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCtg                                                                                 <  1:274329/251‑1 (MQ=255)
                                                                                                                                                                                   aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCAc                                                               >  1:643979/1‑251 (MQ=255)
                                                                                                                                                                                    ggCGGGAGTTAACTTTATTTGCGATAAATACGCGGCCCTCGCCCCCTTCTGCCACGCGGTCAGCGAAGAAGCTGGAGGCACCTCCGGCACGGGCGCGGCCTGCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGTGAAACGCa                                                                <  2:138191/249‑1 (MQ=255)
                                                                                                                                                                                       gggagtgaTCGTTAAGATCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAATCTGGAGTCCCCACCGGCACGGGCGCGGCCTTCCTCGGCTTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCGTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATg                                                           <  2:491974/244‑1 (MQ=255)
                                                                                                                                                                                       ggAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATg                                                           >  1:365235/1‑251 (MQ=255)
                                                                                                                                                                                       ggAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTGCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATg                                                           <  1:195391/251‑1 (MQ=255)
                                                                                                                                                                                              aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCTGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTcaaca                                                     >  1:273960/1‑250 (MQ=255)
                                                                                                                                                                                              aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGCGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGACACGAACCATGTcaaca                                                     >  1:677267/1‑250 (MQ=255)
                                                                                                                                                                                              aTCGTTATGTTCGATAAATACGCGGACGTCGCCCGCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTcaacaa                                                    <  1:97270/251‑1 (MQ=255)
                                                                                                                                                                                                cGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATc                                                  >  1:94456/1‑251 (MQ=255)
                                                                                                                                                                                                 gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCg                                                 >  1:124239/1‑251 (MQ=255)
                                                                                                                                                                                                    aTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcg                                              >  1:100015/1‑251 (MQ=255)
                                                                                                                                                                                                       ttCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAAGGAGTACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCCTGTAAACACTCGGcgtcg                                           >  2:632210/1‑251 (MQ=255)
                                                                                                                                                                                                       gttcgataacgaggcgtccTTCGCCCCCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACGTCGTTCTCGACTGCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACATTCGGcgtc                                            <  2:665849/233‑1 (MQ=255)
                                                                                                                                                                                                         cGATAAATTCGCGGATGTCGCCCTCCTCTGCCCCGCGGTCAGCGAAGAAGCGGGAGGCACCACCGGAGCGGGCGCGGCCTTCCACGGCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTGCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGTCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCgtcgtc                                         <  2:87113/251‑1 (MQ=255)
                                                                                                                                                                                                             aaaTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCgcgg                                     >  1:150621/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAAt                               <  1:243429/251‑1 (MQ=255)
                                                                                                                                                                                                                     cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACAACACCAACGGTGACGTGTACTTAGTTCTCGGCTTCCGCCTACGCGGAGGCGATGGAATACCAACGCGGCGTCAACGGGCGCAGCAGATTAATTAGCGCTTCTGCCTCAAGCTGTGACGGGGAGAAACGCCCCATGTCAACAATCGGCGTCTGCgcggcggaagg                              >  2:516352/1‑246 (MQ=255)
                                                                                                                                                                                                                       gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCt                            >  1:382931/1‑251 (MQ=255)
                                                                                                                                                                                                                       gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCt                            <  1:33588/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCCCCGGCACGTGCGCGGCCTTCCACGACGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATACTGTAAc                       >  2:48207/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGCCCCCTGCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTaa                        <  1:351511/250‑1 (MQ=255)
                                                                                                                                                                                                                             gCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACt                      >  1:232356/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccccTCTTCTTCCAGGCGGTCTGCTAAGAGGCTGGAGGCACCTCCGGCCCGGGCGCGGCCCTCCACGTCGTACCGCACCACACCAACGGTGCCGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACt                      <  2:550965/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tcttctTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGc                  >  1:271969/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tcttctTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGc                  >  1:679180/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tcttctTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCAAGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGc                  >  1:560350/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cttcCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCCGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAAAGCTGTAACTTCGAtttt              >  2:459625/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ccACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCg           >  1:255041/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cgcgGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCCCCACACCACCGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCATCGGCATCAAGCTGTGCCGGGGCGAAACGCACCATGTCAACAATCGGCGACGTCGAGGCGTAATGATGTCACTTCACTTTCTCGcc         >  2:515552/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cgcgGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGccc        >  1:40008/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cAGCGAAGAAGCTGGATGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCCCCACACCAACGGTGACGTGTACTTCGTTCGCGACTTCCGCCTGCGAGGAGGCGATGGAATTCAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGc  <  1:403264/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAGCGAAGAAGCTGGAGGCACCACCGGGCCGGGCGCGGCCTCCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTCCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGc  <  2:461408/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACAAGc  >  2:463110/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGC  >  NC_000913/2890380‑2890870

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