Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,053 A→G V166V (GTT→GTC 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,0530AG98.9% 287.2 / ‑4.2 89V166V (GTT→GTCaraDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (53/35);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

TGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGC  >  NC_000913/65808‑66295
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
tgattataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCg                                                                                                                                                                                                                                              >  2:575058/5‑251 (MQ=255)
  attataGCCTGGTTTCTTTTGATTGGCTGTGGTTTTATACAGTCATTACTGTCCGTAATAGGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCGGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCggg                                                                                                                                                                                                                                            <  2:196205/249‑1 (MQ=255)
   ttataGCCTGGTTGCGTTGGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGGAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATGTTTGCCCATGCAAATGGGCCGTGAGAATGGACCATGACGCCgtg                                                                                                                                                                                                                                            <  2:183663/249‑3 (MQ=255)
         ccTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgc                                                                                                                                                                                                                                     >  1:60906/1‑251 (MQ=255)
         ccTGGTGTCGTTTGAGTGGCTGTGGTTGTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTGGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgc                                                                                                                                                                                                                                     <  2:293315/251‑1 (MQ=255)
           tGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgcgc                                                                                                                                                                                                                                   >  1:81994/1‑251 (MQ=255)
            ggTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCt                                                                                                                                                                                                                                  >  1:432448/1‑251 (MQ=255)
                 cGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGcggc                                                                                                                                                                                                                             >  2:562051/1‑248 (MQ=255)
                 cGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATc                                                                                                                                                                                                                             <  1:211556/251‑1 (MQ=255)
                  gTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCg                                                                                                                                                                                                                            <  1:223510/251‑1 (MQ=255)
                      gATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGCCGCAGGGCATTTGCGCTGCATCGAtcc                                                                                                                                                                                                                        >  2:176328/1‑249 (MQ=255)
                      gATTGGCTGTGGTTTTATACAGTCATTACTACCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTTGTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGTCCTCTTCCAGCACGATGGCGTTATTCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATa                                                                                                                                                                                                                         <  1:700651/250‑1 (MQ=255)
                                 gTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTcg                                                                                                                                                                                                             >  1:490986/1‑250 (MQ=255)
                                  ttttATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTcg                                                                                                                                                                                                             >  2:313620/1‑249 (MQ=255)
                                        aCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGttt                                                                                                                                                                                                      >  1:666789/1‑251 (MQ=255)
                                        aCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGttt                                                                                                                                                                                                      >  1:569142/1‑251 (MQ=255)
                                             cagtACTGCCCGTAATATGCCTTCGTGCCATGCTTACGCAGATAGGGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAAATGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTAc                                                                                                                                                                                                  <  2:376790/247‑1 (MQ=255)
                                             cATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACg                                                                                                                                                                                                 <  1:198701/251‑1 (MQ=255)
                                                    gCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGCCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCAGGCAAATTGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACg                                                                                                                                                                                           <  1:668198/250‑1 (MQ=255)
                                                         tAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCg                                                                                                                                                                                     >  1:589576/1‑251 (MQ=255)
                                                              tGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                                                >  1:483868/1‑251 (MQ=255)
                                                                         cATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGAAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCGTCGCTACCCTGTTTTTCGAAGGTTTCTACGATTACGTTACCGGTTTCCCACTCATAtt                                                                                                                                                                     >  2:580880/1‑251 (MQ=255)
                                                                               tACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGt                                                                                                                                                               >  1:465789/1‑251 (MQ=255)
                                                                                  gCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTgt                                                                                                                                                            >  2:574158/1‑250 (MQ=255)
                                                                                     gataggttttataccgcggcGTGTGCTGAATATTCGGTAATCGCGGCGCTAACTGACGGCAGAATATCCCCATTTAAGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGTCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCTGGCATTTGCGCGGCATCGATACCCTGTCTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACGCATATTCACCGTTGAttt                                                                                                                                                         <  2:705456/233‑1 (MQ=255)
                                                                                       tAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                                                       >  1:189197/1‑251 (MQ=255)
                                                                                       tAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                                                       >  1:403651/1‑251 (MQ=255)
                                                                                           gtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGt                                                                                                                                                   >  1:14487/1‑251 (MQ=255)
                                                                                            tttATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTc                                                                                                                                                  >  1:72530/1‑251 (MQ=255)
                                                                                             ttcttcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGAAGCCGAGCATTTGCGCTGCTTCGATACCATGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCTTTGATTTCTGCGTCg                                                                                                                                                 >  2:600551/6‑251 (MQ=255)
                                                                                                     cagcGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCTGCATTTTTGCCAATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACAGGTTTCCCACTCATATTAACCGTTGATTTCTGCGTCGGGCAtttt                                                                                                                                         >  2:121359/1‑251 (MQ=255)
                                                                                                        cGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCg                                                                                                                                      >  1:169451/1‑251 (MQ=255)
                                                                                                         gTTTGCTGCATATCAGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCACCCTCCTCTTCCCCGTTGATTTCTGCGTCGGTCATTTTTCgg                                                                                                                                     >  2:565022/1‑251 (MQ=255)
                                                                                                          tttGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCggg                                                                                                                                    >  2:692410/1‑251 (MQ=255)
                                                                                                             gctgcatagctggcaccTGCGGCGCTGAACGACGGAAGAATATCCCCATATAAGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGAAAATGGGCCGTGAGAATGGATCAGGACGCCGGGTATTGGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTGTCCCACTCATATTCACCGTTGATTGCTGCGTAGGTCATTTTTCGGGTGc                                                                                                                                 <  2:589576/235‑1 (MQ=255)
                                                                                                             gctgcATAACCGGTAACTGGGGCGCTAACTGACGGCAGAATACCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                                                                                                                  <  2:214703/251‑1 (MQ=255)
                                                                                                                    ttCCGGTACCTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTGCGGCATTTCTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                                           <  1:370001/249‑1 (MQ=255)
                                                                                                                    atCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                                           <  1:71802/250‑1 (MQ=255)
                                                                                                                    atCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAt                                                                                                                          <  1:147062/251‑1 (MQ=255)
                                                                                                                         gTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGc                                                                                                                     >  1:225871/1‑251 (MQ=255)
                                                                                                                         gTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGc                                                                                                                     >  1:461571/1‑251 (MQ=255)
                                                                                                                           aaCTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGCCCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTAGGCATTTTTGCCCATGCAAATGGGACGTGAGACTGGACCAGGCCGCCGGGCATTTGCTCTGCATCGATACCCAGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTAATCTTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCg                                                                                                                   >  2:610535/1‑251 (MQ=255)
                                                                                                                           aaCTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCg                                                                                                                   >  1:18934/1‑251 (MQ=255)
                                                                                                                                  gcgcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCTGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAAt                                                                                                            <  1:363485/251‑1 (MQ=255)
                                                                                                                                           tGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCAACGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGACTGGACCAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAACTAGTCGGCg                                                                                                    >  1:460931/1‑251 (MQ=255)
                                                                                                                                            gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGc                                                                                                     <  1:313620/249‑1 (MQ=255)
                                                                                                                                             aCGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtg                                                                                                  <  1:233684/251‑1 (MQ=255)
                                                                                                                                             aCGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGAAGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGACGTTACCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATATTCGGCgtg                                                                                                  >  1:234427/1‑251 (MQ=255)
                                                                                                                                                gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggg                                                                                               >  1:703321/1‑250 (MQ=255)
                                                                                                                                                gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATATTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTATGCGTCGGTCATTTTGCGGGTGCAGGGAACGGTGCCGTAGAAATAGTCGGCgtgggt                                                                                               >  2:471535/1‑251 (MQ=255)
                                                                                                                                                 cAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGGCGGCgtggggg                                                                                              >  2:219508/1‑249 (MQ=255)
                                                                                                                                                 cAGAATATCCCCATATAAGCGACCTCTACCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCAGATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtg                                                                                              <  2:473885/251‑1 (MQ=255)
                                                                                                                                                 cAGAATATCCCCACATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCTGCAAATGGGCCGTGAGAATGGCCCAGGACGCCGCGCATTTGCGCTGCCTCGATACCATGTTTTTCGACGGTTTCTACGCTGGCGTTACCGGTTTCCCAATCATATTCCCCGTTGATTTCTGCTTCGGTCCTTTTGCGGGTGCAGGGCTTGGTGCCGTAGAAATAGTCGGCgtggg                                                                                                >  2:251328/1‑249 (MQ=255)
                                                                                                                                                  aGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtgg                                                                                             >  1:190757/1‑251 (MQ=255)
                                                                                                                                                   gAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtggg                                                                                            >  1:353425/1‑250 (MQ=255)
                                                                                                                                                    aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGgtggtg                                                                                           <  2:190757/251‑1 (MQ=255)
                                                                                                                                                      tatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCATGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTTCGGGTGCAGGGAATGGTGCCGAAGAATTAGTCGGCGTGGGTGGTGcc                                                                                         >  2:225507/1‑251 (MQ=255)
                                                                                                                                                       atCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCgggggggggtccg                                                                                        >  2:140388/1‑243 (MQ=255)
                                                                                                                                                         ccccACATAAGCGACCGCTTCCAGCTCGATGGCGTTATGCACCGCATCCTCGGCATTTTTGCCCATTCAAATGGGCCGTTAGAACGTACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGGTACCGGTGTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGGAGAAATAGTCGTCGTGGGTGGTGCCgg                                                                                       <  2:709255/250‑1 (MQ=255)
                                                                                                                                                           ccATATAAGAGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGCTGATGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGCCGGTCATTTTGCGGGTGCGGGGGATGGTGCCGAGGAAATAGGCGGCGTGGGTGGTGCCGGTTg                                                                                    >  2:177392/1‑251 (MQ=255)
                                                                                                                                                              tataAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGAACAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGGTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTg                                                                                 >  1:51207/1‑251 (MQ=255)
                                                                                                                                                                 aaGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGa                                                                               >  1:391007/1‑250 (MQ=255)
                                                                                                                                                                  aGCGACCTCTGCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAAt                                                                             <  1:196256/251‑1 (MQ=255)
                                                                                                                                                                       ccTCTTCCAGCACGATGTCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACTCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTTCCGGTTTCCCACTCATATTCTCCGTTGATTTCTGCGTCGGTCATGTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                                        <  2:707532/251‑1 (MQ=255)
                                                                                                                                                                       ccTCGTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                                        <  2:265468/251‑1 (MQ=255)
                                                                                                                                                                           ttCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTggcc                                                                    >  1:606641/1‑248 (MQ=255)
                                                                                                                                                                              cAGCACGATGGCGTTATGCACCGCATCTTCGACCTTTTTGCCCATGCAAATGGGCCGTGAGCATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                                                 >  2:296800/1‑251 (MQ=255)
                                                                                                                                                                                   cGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                                                <  1:620503/135‑1 (MQ=255)
                                                                                                                                                                                   cGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTc                                                                                                                                                                                >  2:620503/1‑135 (MQ=255)
                                                                                                                                                                                     aTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGcc                                                          >  1:119474/1‑251 (MQ=255)
                                                                                                                                                                                           ttATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCa                                                                                                                                                                             >  2:50392/1‑130 (MQ=255)
                                                                                                                                                                                           ttATGCACCGCATCTTCGGCATGTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCa                                                                                                                                                                             <  1:50392/130‑1 (MQ=255)
                                                                                                                                                                                             aTGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGAtagg                                                  >  2:343054/1‑248 (MQ=255)
                                                                                                                                                                                             aTGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGAtag                                                   <  1:653594/250‑3 (MQ=255)
                                                                                                                                                                                               gCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCCGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccccggtcgtgg                                                >  2:671725/1‑242 (MQ=255)
                                                                                                                                                                                                   cgcaactgcGGCATTTTTGCCCAGTCAAATGGGCCGTGAGACTGGACCAGGCCGCCTGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTGGATTTCTGCGTCGGGCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                                            <  2:379234/243‑1 (MQ=255)
                                                                                                                                                                                                     cATCTTCGGCATTTTTGCCCATGCAAATTGGCCGCTAGAATGGACCAGGACGCCGGGCGTTTGCGCTTCTTCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATTTTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGATTGATGCCGTAGAAATAGTCGGCGTGGGTGGTTCGGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAgtggcgtggg                                          <  2:501655/251‑2 (MQ=255)
                                                                                                                                                                                                     cATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggc                                          <  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)
                                                                                                                                                                                                                              aaatgggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGTGTGCACAATgccgc                  >  2:653512/5‑251 (MQ=255)
                                                                                                                                                                                                                                atgggCCGTGAGAATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCa                <  1:343054/249‑1 (MQ=255)
                                                                                                                                                                                                                                               ggACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACTATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGc  <  1:368701/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGTCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCAGTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                          <  2:84239/160‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                          >  1:84239/1‑160 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
TGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGC  >  NC_000913/65808‑66295

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