Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,134 G→A G139G (GGC→GGT)  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,1340GA96.7% 271.4 / NA 95G139G (GGC→GGT) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (0/2);  major base A (51/37);  minor base C (3/0);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCAC  >  NC_000913/65887‑66379
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
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:461571/1‑251 (MQ=255)
                                          gTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGc                                                                                                                                                                                                         >  1:225871/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)
                                                                                                      aTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGcc                                                                                                                                              >  1:119474/1‑251 (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                                                                                                                               >  2:334162/1‑231 (MQ=255)
                                                                                                                                          cATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtgg                                                                                                                               <  1:334162/231‑1 (MQ=255)
                                                                                                                                              cAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGTTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccg                                                                                                       <  1:488008/251‑1 (MQ=255)
                                                                                                                                              cAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGACGTTACCGGTTTCCCACTCATATTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGCATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGCGGCGCGGGTgggggggcccaatgccc                                                                                                       >  2:32353/1‑235 (MQ=255)
                                                                                                                                               aaatgggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGTGTGCACAATgccgc                                                                                                      >  2:653512/5‑251 (MQ=255)
                                                                                                                                                 atgggCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCa                                                                                                    <  1:343054/249‑1 (MQ=255)
                                                                                                                                                                ggACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACTATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGc                                                                                      <  1:368701/250‑1 (MQ=255)
                                                                                                                                                                 cagacggtgtgctcttccgatctcTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCt                                                                                    <  1:535673/228‑1 (MQ=255)
                                                                                                                                                                  aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGTCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAGTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                                                                                                              <  2:84239/160‑1 (MQ=255)
                                                                                                                                                                  aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                                                                                                              >  1:84239/1‑160 (MQ=255)
                                                                                                                                                                       gaCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCGCCAATGGACGGGGATGCCTGAt                                                                                 >  1:115399/2‑248 (MQ=255)
                                                                                                                                                                             gggCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCGGCGCCCAGATAGGGGCGTGGCGCGAGGGGGTGTGCCCAATGCCGCCAATGGCCGGGAATGCCttgtagaggaggg                                                                        >  2:11960/1‑239 (MQ=255)
                                                                                                                                                                              ggactttgctctgcatggtaccccTGTGTTTCGAAGGTTTCTACGAGAACGTTACCGGTTTCCCCCTCATATTCGCCGTTTATTTCTGGGTCGGTCATTTTGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGGCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGagcagtcg                                                                       <  2:569142/230‑4 (MQ=255)
                                                                                                                                                                              ggCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGagcagtcg                                                                       <  1:140388/251‑4 (MQ=255)
                                                                                                                                                                                       gcTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCCGTCGGTGAGTTgg                                                              >  1:235005/1‑251 (MQ=255)
                                                                                                                                                                                       gcTGCATCGATACCCTGTTTTTCGAAGGTTGCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTgg                                                              <  1:143621/251‑1 (MQ=255)
                                                                                                                                                                                        cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGTGTGTGTGCACAATGCCGCCAATGGACGGGAAGGCCtcgatcggaagagctgcgggtagg                                                             >  2:535673/1‑228 (MQ=255)
                                                                                                                                                                                        cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATTGAGCAGTCGGGGGGTTggg                                                             >  1:290831/1‑250 (MQ=255)
                                                                                                                                                                                            aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat                                                         <  1:575058/251‑3 (MQ=255)
                                                                                                                                                                                              cGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGTCCGCCTGCGCCCAGAtag                                                                                                                                       >  1:396011/1‑169 (MQ=255)
                                                                                                                                                                                              cGATACCCTGTTTTTCGAAGGTTTCTACGAGGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAGTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGTCCGCCTGCGCCCAGAtag                                                                                                                                       <  2:396011/171‑3 (MQ=255)
                                                                                                                                                                                                 tACCCTGTTGTTCGAAGGTTTCTACGATGACGTTGCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGAGGCGTGGCGCGGGCGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGTTAGAGCAGTCGGGGAGTTGGCggatcgggg                                                    >  1:204996/1‑243 (MQ=255)
                                                                                                                                                                                                  aCCCTGTTTTCCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggag                                                    <  1:397596/250‑1 (MQ=255)
                                                                                                                                                                                                           tcagagggTTTCTACGATGTCTTTTCCGTTTTCCCACACATATTCACCTTTGATTTCTGCGTCGGTCTTTTTGCGGGTGCGGGGAATGGTGCCGTAGAAATAGTCGGCGTAGGTGGTGCCGGTTGCTGGAATCGACTGCCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGAGAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCttt                                           <  2:666789/244‑1 (MQ=255)
                                                                                                                                                                                                             cgaaGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTc                                        <  1:121359/249‑1 (MQ=255)
                                                                                                                                                                                                             cgaaGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTc                                        <  1:219508/249‑1 (MQ=255)
                                                                                                                                                                                                               aaGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATTGAGCAGTCGGTGGGGTGGCGGATCGGGGGAGGGCTTTTTCgc                                      >  1:275749/1‑250 (MQ=255)
                                                                                                                                                                                                                  gTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCAcc                                   <  1:562051/251‑1 (MQ=255)
                                                                                                                                                                                                                     tCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCtt                                 <  1:574158/250‑1 (MQ=255)
                                                                                                                                                                                                                       tACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATc                                                                                                                                                                >  2:646685/1‑121 (MQ=255)
                                                                                                                                                                                                                       tACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGGCTGGCCCGCCTGCGCCCAGATCGTGGCGTGGCGCGGGTGTGTGTGCACCATGCCGCCAATGGACGGGAATGCCTGATCGAGAAGTCGGTGGGTTGGCGTATCGGAGGAGGGATTTTTCGCCCCTTCa                               >  2:224452/1‑250 (MQ=255)
                                                                                                                                                                                                                       tACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAGTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATc                                                                                                                                                                <  1:646685/121‑1 (MQ=255)
                                                                                                                                                                                                                        aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCa                               <  1:251328/249‑1 (MQ=255)
                                                                                                                                                                                                                          gATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTgg                                                              >  1:57077/1‑216 (MQ=255)
                                                                                                                                                                                                                          gATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTGGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGTTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTgg                                                              <  2:57077/216‑1 (MQ=255)
                                                                                                                                                                                                                           aTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGGCGGTGGGGTGGCGTATCGGAGGGGGGCTTTTTCTCACCTTCAACCAc                          >  2:253672/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCgg                  >  1:173880/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGttt               >  1:514228/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGAAGTCGGGGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCg             >  1:408807/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 cATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGCGGTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTTCACAGGTTTCGATGCTAAcc    >  1:651535/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    atTCTCCGTTGATTTCGGCGGCGGTCAGTTGGCGGGTGCAGGGAATGGGGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCAc  <  2:461571/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
TACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCT‑AACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCAC  >  NC_000913/65887‑66379

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