Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,032 C→T P173P (CCG→CCA)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,0320CT98.9% 297.8 / ‑4.5 89P173P (CCG→CCA) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (54/34);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

CAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCG  >  NC_000913/65794‑66278
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
cAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGGTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGaa                                                                                                                                                                                                                                            <  1:570850/250‑1 (MQ=255)
 aaaaaaCCAGGCTTGATTATAGCCTGGTTTCGTTTGTTTGGCTGAGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGAAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGACTATCCCCCTATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCTCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTgagactg                                                                                                                                                                                                                                          >  2:362852/1‑245 (MQ=255)
    aaaCCAGGCTTGATTATATCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAGTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGAc                                                                                                                                                                                                                                       <  2:379925/251‑1 (MQ=255)
         aGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAgg                                                                                                                                                                                                                                   >  1:37232/1‑249 (MQ=255)
          ggCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGATTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAggac                                                                                                                                                                                                                                 >  2:307966/1‑248 (MQ=255)
              tgattataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCg                                                                                                                                                                                                                             >  2:575058/5‑251 (MQ=255)
                attataGCCTGGTTTCTTTTGATTGGCTGTGGTTTTATACAGTCATTACTGTCCGTAATAGGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCGGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCggg                                                                                                                                                                                                                           <  2:196205/249‑1 (MQ=255)
                 ttataGCCTGGTTGCGTTGGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGGAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATGTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCATGACGCCgtg                                                                                                                                                                                                                           <  2:183663/249‑3 (MQ=255)
                       ccTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgc                                                                                                                                                                                                                    >  1:60906/1‑251 (MQ=255)
                       ccTGGTGTCGTTTGAGTGGCTGTGGTTGTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTGGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgc                                                                                                                                                                                                                    <  2:293315/251‑1 (MQ=255)
                         tGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgcgc                                                                                                                                                                                                                  >  1:81994/1‑251 (MQ=255)
                          ggTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCt                                                                                                                                                                                                                 >  1:432448/1‑251 (MQ=255)
                               cGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGcggc                                                                                                                                                                                                            >  2:562051/1‑248 (MQ=255)
                               cGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATc                                                                                                                                                                                                            <  1:211556/251‑1 (MQ=255)
                                gTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCg                                                                                                                                                                                                           <  1:223510/251‑1 (MQ=255)
                                    gATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGCCGCAGGGCATTTGCGCTGCATCGAtcc                                                                                                                                                                                                       >  2:176328/1‑249 (MQ=255)
                                    gATTGGCTGTGGTTTTATACAGTCATTACTACCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTTGTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGTCCTCTTCCAGCACGATGGCGTTATTCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATa                                                                                                                                                                                                        <  1:700651/250‑1 (MQ=255)
                                               gTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTcg                                                                                                                                                                                            >  1:490986/1‑250 (MQ=255)
                                                ttttATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTcg                                                                                                                                                                                            >  2:313620/1‑249 (MQ=255)
                                                      aCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGttt                                                                                                                                                                                     >  1:666789/1‑251 (MQ=255)
                                                      aCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGttt                                                                                                                                                                                     >  1:569142/1‑251 (MQ=255)
                                                           cagtACTGCCCGTAATATGCCTTCGTGCCATGCTTACGCAGATAGGGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAAATGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTAc                                                                                                                                                                                 <  2:376790/247‑1 (MQ=255)
                                                           cATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACg                                                                                                                                                                                <  1:198701/251‑1 (MQ=255)
                                                                  gCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGCCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑AGGCAAATTGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACg                                                                                                                                                                          <  1:668198/250‑1 (MQ=255)
                                                                       tAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCg                                                                                                                                                                    >  1:589576/1‑251 (MQ=255)
                                                                            tGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                               >  1:483868/1‑251 (MQ=255)
                                                                                       cATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGAAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCGTCGCTACCCTGTTTTTCGAAGGTTTCTACGATTACGTTACCGGTTTCCCACTCATAtt                                                                                                                                                    >  2:580880/1‑251 (MQ=255)
                                                                                             tACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGt                                                                                                                                              >  1:465789/1‑251 (MQ=255)
                                                                                                gCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTgt                                                                                                                                           >  2:574158/1‑250 (MQ=255)
                                                                                                   gataggttttataccgcggcGTGTGCTGAATATTCGGTAATCGCGGCGCT‑AACTGACGGCAGAATATCCCCATTTAAGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGTCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCTGGCATTTGCGCGGCATCGATACCCTGTCTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACGCATATTCACCGTTGAttt                                                                                                                                        <  2:705456/233‑1 (MQ=255)
                                                                                                     tAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                                      >  1:189197/1‑251 (MQ=255)
                                                                                                     tAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                                      >  1:403651/1‑251 (MQ=255)
                                                                                                         gtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGt                                                                                                                                  >  1:14487/1‑251 (MQ=255)
                                                                                                          tttATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTc                                                                                                                                 >  1:72530/1‑251 (MQ=255)
                                                                                                           ttcttcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGAAGCCGAGCATTTGCGCTGCTTCGATACCATGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCTTTGATTTCTGCGTCg                                                                                                                                >  2:600551/6‑251 (MQ=255)
                                                                                                                   cagcGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCTGCATTTTTGCCA‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACAGGTTTCCCACTCATATTAACCGTTGATTTCTGCGTCGGGCAtttt                                                                                                                        >  2:121359/1‑251 (MQ=255)
                                                                                                                      cGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCg                                                                                                                     >  1:169451/1‑251 (MQ=255)
                                                                                                                       gTTTGCTGCATATCAGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCACCCTCCTCTTCCCCGTTGATTTCTGCGTCGGTCATTTTTCgg                                                                                                                    >  2:565022/1‑251 (MQ=255)
                                                                                                                        tttGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCggg                                                                                                                   >  2:692410/1‑251 (MQ=255)
                                                                                                                           gctgcatagctggcaccTGCGGCGCTGAAC‑GACGGAAGAATATCCCCATATAAGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGAAAATGGGCCGTGAGAATGGATCAGGACGCCGGGTATTGGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTGTCCCACTCATATTCACCGTTGATTGCTGCGTAGGTCATTTTTCGGGTGc                                                                                                                <  2:589576/235‑1 (MQ=255)
                                                                                                                           gctgcATAACCGGTAACTGGGGCGCT‑AACTGACGGCAGAATACCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                                                                                                 <  2:214703/251‑1 (MQ=255)
                                                                                                                                  ttCCGGTACCTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTGCGGCATTTCTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                          <  1:370001/249‑1 (MQ=255)
                                                                                                                                  atCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                          <  1:71802/250‑1 (MQ=255)
                                                                                                                                  atCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAt                                                                                                         <  1:147062/251‑1 (MQ=255)
                                                                                                                                       gTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGc                                                                                                    >  1:225871/1‑251 (MQ=255)
                                                                                                                                       gTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGc                                                                                                    >  1:461571/1‑251 (MQ=255)
                                                                                                                                         aaCTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGCCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTAGGCATTTTTGCCC‑ATGCAAATGGGACGTGAGACTGGACCAGGCCGCCGGGCATTTGCTCTGCATCGATACCCAGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTAATCTTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCg                                                                                                  >  2:610535/1‑251 (MQ=255)
                                                                                                                                         aaCTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCg                                                                                                  >  1:18934/1‑251 (MQ=255)
                                                                                                                                                gcgcT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCTGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAAt                                                                                           <  1:363485/251‑1 (MQ=255)
                                                                                                                                                         tGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCAACGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGACTGGACCAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAACTAGTCGGCg                                                                                   >  1:460931/1‑251 (MQ=255)
                                                                                                                                                          gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGc                                                                                    <  1:313620/249‑1 (MQ=255)
                                                                                                                                                           aCGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtg                                                                                 <  1:233684/251‑1 (MQ=255)
                                                                                                                                                           aCGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGAAGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGACGTTACCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATATTCGGCgtg                                                                                 >  1:234427/1‑251 (MQ=255)
                                                                                                                                                              gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggg                                                                              >  1:703321/1‑250 (MQ=255)
                                                                                                                                                              gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATATTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTATGCGTCGGTCATTTTGCGGGTGCAGGGAACGGTGCCGTAGAAATAGTCGGCgtgggt                                                                              >  2:471535/1‑251 (MQ=255)
                                                                                                                                                               cAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGGCGGCgtggggg                                                                             >  2:219508/1‑249 (MQ=255)
                                                                                                                                                               cAGAATATCCCCATATAAGCGACCTCTACCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCAGATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtg                                                                             <  2:473885/251‑1 (MQ=255)
                                                                                                                                                               cAGAATATCCCCACATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCC‑TGCAAATGGGCCGTGAGAATGGCCCAGGACGCCGCGCATTTGCGCTGCCTCGATACCATGTTTTTCGACGGTTTCTACGCTGGCGTTACCGGTTTCCCAATCATATTCCCCGTTGATTTCTGCTTCGGTCCTTTTGCGGGTGCAGGGCTTGGTGCCGTAGAAATAGTCGGCgtggg                                                                               >  2:251328/1‑249 (MQ=255)
                                                                                                                                                                aGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtgg                                                                            >  1:190757/1‑251 (MQ=255)
                                                                                                                                                                 gAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtggg                                                                           >  1:353425/1‑250 (MQ=255)
                                                                                                                                                                  aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGgtggtg                                                                          <  2:190757/251‑1 (MQ=255)
                                                                                                                                                                    tatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCATGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTTCGGGTGCAGGGAATGGTGCCGAAGAATTAGTCGGCGTGGGTGGTGcc                                                                        >  2:225507/1‑251 (MQ=255)
                                                                                                                                                                     atCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCgggggggggtccg                                                                       >  2:140388/1‑243 (MQ=255)
                                                                                                                                                                       ccccACATAAGCGACCGCTTCCAGCTCGATGGCGTTATGCACCGCATCCTCGGCATTTTTGCCC‑ATTCAAATGGGCCGTTAGAACGTACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGGTACCGGTGTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGGAGAAATAGTCGTCGTGGGTGGTGCCgg                                                                      <  2:709255/250‑1 (MQ=255)
                                                                                                                                                                         ccATATAAGAGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGATGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGCCGGTCATTTTGCGGGTGCGGGGGATGGTGCCGAGGAAATAGGCGGCGTGGGTGGTGCCGGTTg                                                                   >  2:177392/1‑251 (MQ=255)
                                                                                                                                                                            tataAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGAACAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGGTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTg                                                                >  1:51207/1‑251 (MQ=255)
                                                                                                                                                                               aaGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGa                                                              >  1:391007/1‑250 (MQ=255)
                                                                                                                                                                                aGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAAt                                                            <  1:196256/251‑1 (MQ=255)
                                                                                                                                                                                     ccTCTTCCAGCACGATGTCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACTCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTTCCGGTTTCCCACTCATATTCTCCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                       <  2:707532/251‑1 (MQ=255)
                                                                                                                                                                                     ccTCGTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                       <  2:265468/251‑1 (MQ=255)
                                                                                                                                                                                         ttCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTggcc                                                   >  1:606641/1‑248 (MQ=255)
                                                                                                                                                                                            cAGCACGATGGCGTTATGCACCGCATCTTCGACCTTTTTGCCC‑ATGCAAATGGGCCGTGAGCATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                                >  2:296800/1‑251 (MQ=255)
                                                                                                                                                                                                 cGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                               <  1:620503/135‑1 (MQ=255)
                                                                                                                                                                                                 cGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                               >  2:620503/1‑135 (MQ=255)
                                                                                                                                                                                                   aTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGcc                                         >  1:119474/1‑251 (MQ=255)
                                                                                                                                                                                                         ttATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCa                                                                                                                                                            >  2:50392/1‑130 (MQ=255)
                                                                                                                                                                                                         ttATGCACCGCATCTTCGGCATGTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCa                                                                                                                                                            <  1:50392/130‑1 (MQ=255)
                                                                                                                                                                                                           aTGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGAtagg                                 >  2:343054/1‑248 (MQ=255)
                                                                                                                                                                                                           aTGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGAtag                                  <  1:653594/250‑3 (MQ=255)
                                                                                                                                                                                                             gCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCCGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccccggtcgtgg                               >  2:671725/1‑242 (MQ=255)
                                                                                                                                                                                                                 cgcaactgcGGCATTTTTGCCC‑AGTCAAATGGGCCGTGAGACTGGACCAGGCCGCCTGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTGGATTTCTGCGTCGGGCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                           <  2:379234/243‑1 (MQ=255)
                                                                                                                                                                                                                   cATCTTCGGCATTTTTGCCC‑ATGCAAATTGGCCGCTAGAATGGACCAGGACGCCGGGCGTTTGCGCTTCTTCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATTTTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGATTGATGCCGTAGAAATAGTCGGCGTGGGTGGTTCGGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAgtggcgtggg                         <  2:501655/251‑2 (MQ=255)
                                                                                                                                                                                                                   cATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggc                         <  1:540536/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tgcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgtgtgc           <  1:307966/248‑8 (MQ=255)
                                                                                                                                                                                                                                      ccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgtgtgccca        >  1:637661/1‑240 (MQ=255)
                                                                                                                                                                                                                                       cATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtgg                          <  1:334162/231‑1 (MQ=255)
                                                                                                                                                                                                                                       cATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtgg                          >  2:334162/1‑231 (MQ=255)
                                                                                                                                                                                                                                           cAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGTTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccg  <  1:488008/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGACGTTACCGGTTTCCCACTCATATTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGCATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGCGGCGCGGGTgggggggcccaatgccc  >  2:32353/1‑235 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCG  >  NC_000913/65794‑66278

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