Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,832 T→C 4.4% H312R (CAT→CGT)  nanM ← N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,8320TC 4.4% 322.7 / 2.7 91H312R (CAT→CGT) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base T (48/39);  new base C (2/2);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTT  >  NC_000913/4537550‑4538114
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                          
cAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaaa                                                                                                                                                                                                                                                                                      >  2:319633/1‑289 (MQ=255)
 accGGCAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaa                                                                                                                                                                                                                                                                                       <  2:804349/285‑1 (MQ=255)
 aaCGGTAAGAACATAGTAATAATCTGGCGATATTATTGCGTTCATATAAATTTCGCGTACTTGAAATGAGAAATTGATAGCAGGGCCGCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGGGACTTTATTATCCTTAACAGTGATCAAAACTGAATCCGTCACCGCGTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAg                                                                                                                                                                                                                                                                                     <  1:637073/289‑1 (MQ=255)
 aaCGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaa                                                                                                                                                                                                                                                                                       <  1:393556/287‑1 (MQ=255)
 aaCGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaa                                                                                                                                                                                                                                                                                       <  1:416771/287‑1 (MQ=255)
    ggTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATg                                                                                                                                                                                                                                                                                  <  1:574669/289‑1 (MQ=255)
          aaCATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCAATTTCCCGTTATGCCAAAGATGAattcc                                                                                                                                                                                                                                                                            >  1:222073/1‑286 (MQ=255)
           aCATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAatat                                                                                                                                                                                                                                                                             <  1:319633/287‑1 (MQ=255)
              tAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTACGCCAAAGATGAATATCAGt                                                                                                                                                                                                                                                                         >  2:498402/1‑288 (MQ=255)
                       tCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGAt                                                                                                                                                                                                                                                               >  2:563592/1‑289 (MQ=255)
                         gggCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTTCGATAATGCACCCGATTTATCCCATTTCCCGTTATGCCAAAGCTGCACATCAGTGCTCTATGAttt                                                                                                                                                                                                                                                             >  1:662638/1‑289 (MQ=255)
                           gCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTACATGAttttt                                                                                                                                                                                                                                                           >  2:17564/1‑289 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACACATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCCTTAGTTTTGTACTGTGGCTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGACGTTTCTCCGCCAATAATCAATCGACTATTATTCCAGGGCAATGCTACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCGGTTATGCCAAAGATGAATATCAGTGCTATATGTTTTTTCCAGGCCTTCATgc                                                                                                                                                                                                                                             >  2:157512/1‑289 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACTGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATgc                                                                                                                                                                                                                                             >  2:820674/1‑289 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAg                                                                                                                                                                                                                                                                                     >  1:321535/1‑249 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAg                                                                                                                                                                                                                                                                                     <  2:321535/249‑1 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATgc                                                                                                                                                                                                                                             >  2:5479/1‑289 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATgc                                                                                                                                                                                                                                             >  2:234599/1‑289 (MQ=255)
                                         ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATg                                                                                                                                                                                                                                              >  1:424719/1‑288 (MQ=255)
                                           cATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATgcg                                                                                                                                                                                                                                            >  2:388806/1‑288 (MQ=255)
                                            atatAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTACGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCa                                                                                                                                                                                                                                          <  1:498402/289‑1 (MQ=255)
                                                  aTTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGt                                                                                                                                                                                                                                      >  1:89157/1‑287 (MQ=255)
                                                     tCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTa                                                                                                                                                                                                                                 <  1:563592/289‑1 (MQ=255)
                                                      cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTa                                                                                                                                                                                                                                 >  1:135757/1‑288 (MQ=255)
                                                      cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACACGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATTCGCATAGTTCttc                                                                                                                                                                                                                                 >  1:509640/1‑287 (MQ=255)
                                                           aaTTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTc                                                                                                                                                                                                                           <  2:141105/289‑1 (MQ=255)
                                                                        aaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTtctc                                                                                                                                                                                                              <  1:708501/289‑1 (MQ=255)
                                                                           tGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTg                                                                                                                                                                                                           <  2:837402/289‑1 (MQ=255)
                                                                               aGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGaa                                                                                                                                                                                                         >  2:727132/1‑287 (MQ=255)
                                                                                                ggggCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCt                                                                                                                                                                                       >  1:576917/1‑288 (MQ=255)
                                                                                                   gCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGCTTTTTTCCGGCCTTCCTGCGCATAGTTCTTACCGTTCTGGTATTTTTCTCGGTAATCTTTTAATCCGTCCCCCCCg                                                                                                                                                                                    >  1:385486/1‑288 (MQ=255)
                                                                                                    cccTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  <  1:727132/289‑1 (MQ=255)
                                                                                                     ccTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACTGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  <  1:820674/288‑1 (MQ=255)
                                                                                                              ttAGTGTTGTACTGTGAGTTTATCATCCTTCACCGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCATTTTCACCGCAAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCTGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAag                                                                                                                                                                         <  1:62876/288‑1 (MQ=255)
                                                                                                                aGTTTTGTACTGTGACTTTTTTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGACTGCAGTTTCATCGCCAATAATCAACAGACTACTAATCCCGGGCAATGATACCCCGTAGGCCCGACTTTGGTATAATTTACCCGACTTATCCGATTTCCCGTTATGCTAATGTTGAGTATCAGTGCTATATGTTTTTTTCAGTCCTTCCTGCGCAAAGTTCTGCCAGTACTGGTATTTTTCTCGTGACCCTTTGAATCCGGCTCCACCGGTAAACTTACgaga                                                                                                                                                                      >  1:396962/1‑289 (MQ=255)
                                                                                                                    ttGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCa                                                                                                                                                                  <  1:834132/289‑1 (MQ=255)
                                                                                                                           gtgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  >  2:4804/1‑266 (MQ=255)
                                                                                                                           gtgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  <  1:4804/266‑1 (MQ=255)
                                                                                                                            tgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAATTATAAGAGAATCATTGCTTa                                                                                                                                                           >  1:138054/1‑288 (MQ=255)
                                                                                                                                      atCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCATTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCCCaaa                                                                                                                                                 >  1:616743/1‑288 (MQ=255)
                                                                                                                                        ccTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTa                                                                                                                                                                                                                                 >  1:728362/1‑206 (MQ=255)
                                                                                                                                        ccTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTa                                                                                                                                                                                                                                 <  2:728362/206‑1 (MQ=255)
                                                                                                                                          ttCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAACCCg                                                                                                                                             >  1:657662/1‑288 (MQ=255)
                                                                                                                                          ttCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCg                                                                                                                                             >  2:654068/1‑288 (MQ=255)
                                                                                                                                            cacaGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCg                                                                                                                                                                                             >  1:669306/1‑238 (MQ=255)
                                                                                                                                            cacaGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCg                                                                                                                                                                                             <  2:669306/238‑1 (MQ=255)
                                                                                                                                                   aTCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACg                                                                                                                                    >  2:297353/1‑288 (MQ=255)
                                                                                                                                                      aaaaCTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTTCTATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCAt                                                                                                                                <  2:662638/289‑1 (MQ=255)
                                                                                                                                                         aCTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTg                                                                                                                             <  2:705298/289‑1 (MQ=255)
                                                                                                                                                             aaTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCa                                                                                                                                                                                                                                                                                        <  1:434336/130‑1 (MQ=255)
                                                                                                                                                             aaTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCa                                                                                                                                                                                                                                                                                        >  2:434336/1‑130 (MQ=255)
                                                                                                                                                             aaTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGtga                                                                                                                         >  1:829915/1‑289 (MQ=255)
                                                                                                                                                              aTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGtgat                                                                                                                        <  2:612972/289‑1 (MQ=255)
                                                                                                                                                                  gTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCa                                                                                                                                                                  <  2:597839/243‑1 (MQ=255)
                                                                                                                                                                  gTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCa                                                                                                                                                                  >  1:597839/1‑243 (MQ=255)
                                                                                                                                                                    cACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGAc                                                                                                                   <  2:657662/288‑1 (MQ=255)
                                                                                                                                                                                  cTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTAtt                                                                                                    <  2:135757/289‑1 (MQ=255)
                                                                                                                                                                                         ttCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCAttt                                                                                              <  2:273182/288‑1 (MQ=255)
                                                                                                                                                                                                  aataatCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTTATGAGACGGGAGCAAGCTTATTCCATTTTAAattatt                                                                                     >  1:595016/1‑288 (MQ=255)
                                                                                                                                                                                                   ataatCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTAc                                                                                   <  2:576917/289‑1 (MQ=255)
                                                                                                                                                                                                              aCTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCa                                                                        >  2:48439/1‑289 (MQ=255)
                                                                                                                                                                                                               cTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATc                                                                         <  2:385486/287‑1 (MQ=255)
                                                                                                                                                                                                                tattatTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAg                                                                      >  2:672003/1‑289 (MQ=255)
                                                                                                                                                                                                                tattatTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAg                                                                      >  2:24087/1‑289 (MQ=255)
                                                                                                                                                                                                                tattatTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAg                                                                      >  1:443987/1‑289 (MQ=255)
                                                                                                                                                                                                                   tattCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  <  2:448488/178‑1 (MQ=255)
                                                                                                                                                                                                                   tattCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGc                                                                                                                                                                                  >  1:448488/1‑178 (MQ=255)
                                                                                                                                                                                                                        cAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCa                                                                                                                                                                                                                                                                                        <  1:26779/71‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCa                                                                                                                                                                                                                                                                                        >  2:26779/1‑71 (MQ=255)
                                                                                                                                                                                                                        cAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAt                                                               >  2:730235/1‑288 (MQ=255)
                                                                                                                                                                                                                          gggCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTACGCCa                                                                                                                                                                                                                                                                                        <  1:368305/69‑1 (MQ=255)
                                                                                                                                                                                                                          gggCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTACGCCa                                                                                                                                                                                                                                                                                        >  2:368305/1‑69 (MQ=255)
                                                                                                                                                                                                                               aTGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAt                                                        >  2:689641/1‑288 (MQ=255)
                                                                                                                                                                                                                                     cTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTc                                                  >  2:477085/1‑288 (MQ=255)
                                                                                                                                                                                                                                     cTCCGTAGGCCCGACCTTGCCATTATTCACCCGCTTTATCCCAGTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGTCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCg                                                 <  1:672003/289‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCaa                                              >  2:536842/1‑289 (MQ=255)
                                                                                                                                                                                                                                          tAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAt                                             >  1:342606/1‑288 (MQ=255)
                                                                                                                                                                                                                                          tAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAt                                             >  2:282593/1‑288 (MQ=255)
                                                                                                                                                                                                                                           aGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAt                                             >  2:381473/1‑287 (MQ=255)
                                                                                                                                                                                                                                                  aCCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGc                                                                                                             >  2:659821/1‑216 (MQ=255)
                                                                                                                                                                                                                                                  aCCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGc                                                                                                             <  1:659821/216‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATc                                            <  2:856302/281‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATc                                            >  1:856302/1‑281 (MQ=255)
                                                                                                                                                                                                                                                           tAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGcc                            <  2:89157/288‑1 (MQ=255)
                                                                                                                                                                                                                                                             aTTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCAtt                         >  2:637697/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  cccGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAg                                                                                                                                                                                                                                                                          <  1:210532/43‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cccGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAg                                                                                                                                                                                                                                                                          >  2:210532/1‑43 (MQ=255)
                                                                                                                                                                                                                                                                              cATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTAt        <  1:730235/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   cccGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt   <  2:424719/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCtt  <  1:536842/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCtt  <  1:654068/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                          
CAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTGCTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTT  >  NC_000913/4537550‑4538114

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