Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,4680TG35.1% 121.1 / 14.3 80Y211S (TAC→TCC) cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (10/40);  new base G (27/0);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 7.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTC  >  NC_000913/2889219‑2889716
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
cccTGAGTGGACACGACGGCTTCGGCCCCCGCCCGCGTATTCTCCAGCGGCAGATAGCCAAACCCACTCGCGGTGCGGGCGTAGGTTTTCTTGGTGCCGTGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCCCCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCaggta                                                                                                                                                                                                                                                             <  2:535157/249‑1 (MQ=255)
     aGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAACCTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑TCCGTGG‑TTTT‑GAATTTTCGCGGCAGGGAGGTCt                                                                                                                                                                                                                                                        >  1:576253/1‑251 (MQ=255)
      gTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAAATTACGCGGCAAGGAGGTCTg                                                                                                                                                                                                                                                       >  1:476343/1‑251 (MQ=255)
         gtcccgcggctttCGGCCTCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGTTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGTCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGGGG‑TTTT‑GAAGTTACGCGGCAGGTAGGTCTGGcc                                                                                                                                                                                                                                                    <  2:379024/240‑1 (MQ=255)
         gtcccgacgtcctCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACGCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGGTCTGTGGCAGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTTCGCGGCAGGTAGGTCTGGcc                                                                                                                                                                                                                                                    <  2:379060/239‑1 (MQ=255)
               acgGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTT‑ACT‑ACCGTGG‑GTTT‑GAATTTTCGCGGCAGGGAGGGCTgggcgggggg                                                                                                                                                                                                                                              >  1:56540/1‑243 (MQ=255)
                       ggCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCTCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCTGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑ttctt                                                                                                                                                                                                                                     <  2:53418/251‑1 (MQ=255)
                          cACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑CAT‑ACCGTGG‑TTTT‑GAATTTTAGCGGCAGGGAGGGCTGGCCGAGGAGCGG‑gtctttat                                                                                                                                                                                                                                  >  1:1125/1‑248 (MQ=255)
                             cGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTTCGCGGCCGGTAGGTCTgggcgcggggcgggtctttattcg                                                                                                                                                                                                                                >  2:59952/1‑229 (MQ=255)
                              gCCAGCGTATGCTCCAGCGGCAGATAGCCAACCTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCagt                                                                                                                                                                                                                              <  2:21250/251‑1 (MQ=255)
                              gCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCagt                                                                                                                                                                                                                              <  1:563815/251‑1 (MQ=255)
                                      aTGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGAGCGATATCGTTCTGTGTCGGGATC‑ACT‑TCCGTGG‑TTTT‑GAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGa                                                                                                                                                                                                                      <  1:540166/251‑1 (MQ=255)
                                       tGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑TAATTTACGCGGCAGGGATGGGTGGCCGAGGGTCGG‑ttttttttcagttgtgggggc                                                                                                                                                                                                                     >  1:315583/1‑232 (MQ=255)
                                                   aGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGGGG‑GTTT‑TAATTTTCGCGGCCGGGAGGGCTTGGCGGGGGTCGG‑TTTTTTCTCAGGAGTGGCGG‑CT‑TTTTTCTGGTc                                                                                                                                                                                                       >  1:334288/1‑251 (MQ=255)
                                                           aaaCTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTC‑CCT‑ACCATGG‑TGTT‑GTATTTAACCGGCAGGGAGGGGTGGGCGGGGGTAGG‑TTTTTCTTCAATAGTGGCGAGCT‑TTTTC‑TTGTCGAGCCac                                                                                                                                                                                                >  2:411472/1‑249 (MQ=255)
                                                           aaaCTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGa                                                                                                                                                                                               <  1:428723/251‑1 (MQ=255)
                                                                       ggTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATT‑ACT‑ACCGTGG‑TTTT‑TGATTTACGCGGCAGGGAGGGATGGGCGGGGATAGG‑TTCTTTATTAGTAGTGGGGA‑CT‑TTTTTCTTGTCGGGCCAGATCTCCGCCTAcc                                                                                                                                                                                   >  1:441532/1‑250 (MQ=255)
                                                                          gCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATAcgcg                                                                                                                                                                                 <  2:491299/250‑1 (MQ=255)
                                                                           cGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGAGCGCCAGGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑AGCGTGG‑TGTT‑GAAATTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TGCTTCATCAGTAGTGGCGT‑CT‑TTTTCCTGGTCGAACCAGATCTCCGCATGcgcgcg                                                                                                                                                                               <  2:376264/251‑1 (MQ=255)
                                                                            gggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATT‑ACT‑ACCGTGG‑GTTT‑TTATTTTCGCGGCAGGGGGGTCTGGCCGAGGATCCG‑GTCTTCTTCAGTAGTGGCGC‑CC‑TTTTTCTCGTCGAGCCAGAGCTCCGCATACGCGCgg                                                                                                                                                                              >  2:272977/1‑251 (MQ=255)
                                                                            gggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGGAGGTCTGGGCGAGGATCGG‑TTCTTCATCAGTAGGGGCGA‑CT‑TTTTTCTGGTCGAGCCAGATCTCCGCATACGCGCgg                                                                                                                                                                              >  1:15901/1‑251 (MQ=255)
                                                                            gggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCCTGG‑TTTTTGATTTTA‑GCCGCAGGGAGGGCTGGCCGAGGGTCGG‑GTCTTCATTAGTTGGGGCGG‑CT‑TTTTTCTGGTCGAGCCAGATCTCCGCATACCCGCgg                                                                                                                                                                              >  2:549275/1‑251 (MQ=255)
                                                                             ggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGGGG‑TTTT‑GAATTTACGCGGCAGGGAGGTCTGGCCGAGGGTCGG‑TTTTTTATCAGTCGTGGCGG‑CT‑TTTTCCTGGTCGAGCCAGATCCCCGcctccgccccgg                                                                                                                                                                             >  2:287502/1‑240 (MQ=255)
                                                                              gCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCCCCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCTTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTGCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGt                                                                                                                                                                            <  2:56972/251‑1 (MQ=255)
                                                                                gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACC‑ACCGTGG‑GTTT‑GGATTTACGCGGCAGGGAGGGCTGGGCGAGGGTCGG‑TTCTTCCTCAGGAGTGGCGG‑CT‑TTTTTCTGGGCGGGCCAGATCTCCCCCTACGCGCggggac                                                                                                                                                                          >  1:504673/1‑248 (MQ=255)
                                                                                gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTCAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTTCGCGGGGGGGAGGGCTGGGCGGGGGTCGGGTTTTTC‑TCAGTAGGGGGGA‑CT‑TTTTCCTGGTGGAGCCAGATCTCCCCCTATGCGCgggga                                                                                                                                                                           >  1:540795/1‑248 (MQ=255)
                                                                                  aGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCAAGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑CCC‑ATCGTGG‑GTTT‑GAATTTACCCGGAAGGCAAGTCTGGCCGAGGGCCGG‑GTTTTCATTAGTAGTGGCGA‑CTCTTTTCCTG‑TCGAGCCAGATCTCCGCATaatcccgggcacgc                                                                                                                                                                        >  2:239144/1‑238 (MQ=255)
                                                                                     ttttCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTATGCGGCAGGGAGGTCTGGCCGAGGGTCGG‑GTCTTTATCAGTAGGGGGGA‑CT‑TTTTCCTGGTCGGGCCAGATCTTCGCATACGCGCGGGGACGAggg                                                                                                                                                                     >  1:547490/1‑250 (MQ=255)
                                                                                            gttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCCCGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGGGTTTT‑GA‑TTTTCGCGGCAGGGAGGGCTGGCCGGGGGGCGG‑GTTTTCTTCAGGAGTGGGGC‑CT‑TTTTCCTTGTCGAGCCCGATCTCCCCATACGCGCGGGTGGGCGGCAACAGa                                                                                                                                                               >  2:279426/1‑250 (MQ=255)
                                                                                              tgCCGGGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCGGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATg                                                                                                                                                             <  2:576253/250‑1 (MQ=255)
                                                                                               gCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCt                                                                                                                                                           <  1:21892/251‑1 (MQ=255)
                                                                                                 cGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCa                                                                                                                                                         <  1:341717/251‑1 (MQ=255)
                                                                                                  gtgtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAg                                                                                                                                                        <  2:376014/251‑1 (MQ=255)
                                                                                                             aaaGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCCAC‑‑ACCGTGG‑TTTT‑GAATTTACGCGGCCGGTGGGGCTGGCCGGGGGGCGG‑TTCTTCATCAGTAGTGGGGA‑CT‑TTTTTCTGGTCGAGCCAGATCTCCGCACACACGCGGGTACGGGGCAAcgactgtcaagaaaacttcctc                                                                                                                                             >  2:251436/1‑230 (MQ=255)
                                                                                                             aaaGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCttcttc                                                                                                                                             <  2:365063/251‑1 (MQ=255)
                                                                                                              aaGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCg                                                                                                                                            <  1:18878/251‑1 (MQ=255)
                                                                                                                    accgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTACGCGGCAGGGAGGTCTTGCCGAGGATCGG‑TTCTTTATCAGTAGTGGCGA‑CT‑TTTTCCTTGTCGAGCCAGATCTCCCCATACGCGCGGGTACGAGGCAACAAATGCTCCGAAATCTTCTTCCCCCAcc                                                                                                                                      >  1:156964/1‑250 (MQ=255)
                                                                                                                    accgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTACGCGGCAGGGAGGTCTGGCCGAGGGTCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTTCTGGTCGAGGCAGATCTCCGCATACGCGCGGGTACGAGGCACAAGAAGGCCCGAAATCTTCTTTGCCCAcc                                                                                                                                      >  1:478450/1‑250 (MQ=255)
                                                                                                                        cccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGACGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTC‑ACT‑ACCGGGG‑TTTT‑GAATTTTCGCCGCAGGGAGGGCTCGCCGCGGCTCGG‑TTCCTCATTCGTTGGGGGGA‑CT‑TTTTTCTGGTCGAGCCAGATCTCCGCCTACGCGCGGGGACGAGGCAACAGATGCCCAGCAATCTTCTTCGCCCACTCgg                                                                                                                                   >  1:297829/1‑249 (MQ=255)
                                                                                                                               aCAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGGGG‑GTTT‑GAATTTTCGAGGCAGGGAGGTCTGGCCGGGGATCGG‑GTCTTCATCCTTAGTGGCGA‑CT‑TTTTTCTGGTCTAGCCAGATCTCCTCATACGCGCGGGTAAGAGGCAACAGATTCTCCGCAAACTTCTTCTCCCACTCGTACCCTTcc                                                                                                                           >  1:173339/1‑251 (MQ=255)
                                                                                                                                cAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGTCCGAGGATCGG‑TTCGTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCcg                                                                                                                          <  1:528027/251‑1 (MQ=255)
                                                                                                                                cAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCcg                                                                                                                          <  1:163360/251‑1 (MQ=255)
                                                                                                                                cAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCcg                                                                                                                          <  2:489053/251‑1 (MQ=255)
                                                                                                                                     tAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTg                                                                                                                      <  2:147448/250‑1 (MQ=255)
                                                                                                                                                  cTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATA‑ACT‑ACCGTGG‑TTTT‑GAATTTACTCGGAAGGAAGGTCTGGCCGAGGCTCGG‑TTCTTTATCAGTAGAGGCGA‑CT‑TTTTCCTGGTCGGGCCAGATCTTCGCATACCCGCGGGGACGAGGGAAAAGATGCCTAGAAATCTTCTTCCCCCACTCGTACGCTTgcgccgtgagctggcatccgt                                                                                                        >  1:571742/1‑234 (MQ=255)
                                                                                                                                                    ggCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTaa                                                                                                      <  2:56540/250‑1 (MQ=255)
                                                                                                                                                           ttCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTc                                                                                                <  2:312427/250‑1 (MQ=255)
                                                                                                                                                                cGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATT‑ACC‑ACCGTGG‑GTTT‑GAATTTACACGGCCGGTCGGCGTGGCCCCGGATCGG‑TTCTTCATCAGCCATGGCGC‑CT‑TTTTTCTGGTCCAGCCAGCTCTCCGCATACGCGCGGGGACGCGGGGACAGATGCTCAGAAATCTTCTTCGCCCACTCCTACGCTTCCGCGCGCAGCTGCGACTCGTCAGGGGTCgggggg                                                                                          >  1:472219/1‑246 (MQ=255)
                                                                                                                                                                  aTCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACA‑ACCGTGG‑GTTT‑TAATTTACGCGGCAGGGAGGTCTGGGCGAGGATCGG‑GTCTTCCTCAGTAGGGGGGG‑CT‑TTTTCCTGGTCGAGCCAGATTTCCGCATACGCGCGGGGACGGGGCAAAAGATGGCCAAAACTCTTCTTTGCCCCCTCGGACGCTTCTGCGTGCAGCTGCGGCTccttaggggtcggggggg                                                                                         >  1:409010/1‑233 (MQ=255)
                                                                                                                                                                       cACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑CAT‑ACCGTGG‑GTTT‑GAATTTTCGCGGCAGGGAGGGCCGGCCGCGGTGCGG‑TTTTTTTTCCGTAAGGGGGG‑CT‑TTTTTCTGGGCGAACCAGATCTCCGCCTACGGGCGGGGACGAGGCAAAAGATGCTCAGAAAAATTCTTCGGCCCCTCGGACGCTTCCCCGCGCCGCCGCGGCTCGTTAgggggcgagggggcgagg                                                                                    >  2:311376/1‑235 (MQ=255)
                                                                                                                                                                       cACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑TGATTTACGCGGCAGGGAGGGCTGGCCCAGGATCGG‑GTCTTCATCAGTAGTGGCGG‑CT‑TTTTTCTTGTCGAGCCAGATCTCCGCCTACGCGCGGGGGCGGGGCAACAGATGCTCAGAAATCTTTTTCGCCCACTCTTACGCTTCCCCCTGCCGCTCCGGCACGGAAGGGTTCGCGGTGCagcggg                                                                                   >  2:126399/1‑247 (MQ=255)
                                                                                                                                                                             gTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑TGATTTTCGCGGCCGGGAGGTCTGGGCGAGGGTCGG‑TTCTTCATCAGTAGTGGGGG‑CT‑TTTTTCTGGTCGGGCCAGATCTCCCCATACGCGCGGGTAAGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCCTCCGCGTGCCGCTGCGGCTTGGAAGGGGTCGAGGGGGAGGGTACGTT‑Ac                                                                            >  1:75978/1‑251 (MQ=255)
                                                                                                                                                                                      gTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑CCTCACCG‑GG‑ATTT‑TAATTTTCGCGGCAGGGAGGGCTGGCCGCGGAACGG‑GTCTTCATCAGTAGTGGCGA‑CC‑TTTTTCTGGGCGAGCCAGATCTCCCCAAACGCGCCGGGACGAGGCAACCGATGCTCAGAAATCTTCTTCGCCCACACGCACGCCTCCGCGTGCAGCTGCGACTCGTAAGGGGTTGAGGTGCAGAGGACGTTTACG‑TTCATGt                                                                   >  1:577582/1‑251 (MQ=255)
                                                                                                                                                                                         ggCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGt                                                                <  1:558114/251‑1 (MQ=255)
                                                                                                                                                                                         ggCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAAGTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCg                                                                 <  2:472502/250‑1 (MQ=255)
                                                                                                                                                                                             tGCAGTTCGTTATCGTTCTGTGGCGGGATT‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CG‑TTTTCCTGGTCGAGCCAGAGCTCCGCATACGCGCGGGTACGAGTCAACAGATGCTCAGAAATCTTCTTCTCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAATGGTTCGAGTTGCAGAGGACGTT‑CCCGTTCATGTCGTTa                                                              <  2:199753/249‑1 (MQ=255)
                                                                                                                                                                                                aGATCGATATCTTTCTGTGGCGTGATC‑TCT‑ACCGTTG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTGGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGt                                                         <  1:65762/251‑1 (MQ=255)
                                                                                                                                                                                                 gATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑GAATTTACGCGGCCGGGAGGGCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTTCTGGTCGCGCCAGATCTCCGCATACGCGCGGGGACGAGGCAAAAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTTGAGGTGCAGAGGACGTT‑ACGGGTCATGTCGCTAGCTGTc                                                        >  1:453320/1‑251 (MQ=255)
                                                                                                                                                                                                     gATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCa                                                    <  1:142728/251‑1 (MQ=255)
                                                                                                                                                                                                      atatCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑GTTT‑TAATTTACGCGGCAGGGAGGTCTGGCCGAGGATCGG‑TTTTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCCGATCTCCGCATACGCCCGGGTACGAGGCAACAAATGCTCCGAAATCTTCTTCCCCCACTCCTACGCTTCCGCGGGCAGCTGCGACCCGTAAGGGGTCGTGGGGCAGAGTAAGTT‑ACGGTTTATGTCGTTAGGTGTCGcccg                                                   >  1:116078/1‑249 (MQ=255)
                                                                                                                                                                                                            ttCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGt                                                                                <  1:393788/217‑1 (MQ=255)
                                                                                                                                                                                                                 tgGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgac                                        >  2:438534/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgacc                                       <  1:556470/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgacc                                       <  1:445205/251‑1 (MQ=255)
                                                                                                                                                                                                                     gggATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGTTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCGTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCGGTCGCCAGCGCATCAAgaccgac                                    <  2:391221/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgaccgacc                                   >  1:239063/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAgacctcccga                                 <  2:159723/251‑7 (MQ=255)
                                                                                                                                                                                                                          c‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAg                                <  1:125384/250‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgca                            <  1:21363/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cGTGG‑TTTT‑GACTTTACGCGGCAGGTCGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcagc                          >  2:86690/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ttt‑GAATTTACGCGGCAGGCAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTCGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGATGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTTATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGCCCGAGTGCAGCATCTggg                   >  2:62963/1‑250 (MQ=255)
                                                                                                                                                                                                                                           t‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGACAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                 >  2:174592/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaTTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGATGCCACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGGGCCGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGCCCTAGTGCAGCACCTGGTGCAc               >  2:489998/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tACGCGGCAGGTAGGTCTGGCCGAGGCTCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGt           >  2:265874/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGt           >  1:555927/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tACGCGGCAGGTAGGTCGGGCCGAGTGTCGG‑TTCTTCAGCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGGGCAGCATCTGGTGCACCgg            <  2:339734/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt          <  2:360806/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgcgGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCAGCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGttt         <  2:533916/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     gcgGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTg        <  1:390942/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     gcgGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTg        <  2:105305/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc      <  2:156964/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           ggtaggtCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTc  <  1:438281/251‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
CGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATC‑ACT‑ACCGTGG‑TTTT‑GAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGG‑TTCTTCATCAGTAGTGGCGA‑CT‑TTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTT‑ACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTC  >  NC_000913/2889219‑2889716

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