Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,509 A→G A391A (GCA→GCG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,5090AG97.3% 224.0 / ‑3.3 73A391A (GCA→GCG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/36);  total (36/37)
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.52e-01

TTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCACTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGA  >  NC_000913/1263‑1747
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ttCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTgg                                                                                                                                                                                                                                                     >  1:552754/1‑250 (MQ=255)
 tCCAATCTGAATAATATGGCAAGGTTCATCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTTCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGa                                                                                                                                                                                                                                                     <  1:394236/250‑1 (MQ=255)
  ccAATCTTAATAATATGTCAATGTGCAGCGTTTCCGGCCCGGGGATTAAAGGGATGGTTTGCATGGGGGCGCGCGTGATTTCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCACTCATCATCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTTGCGGTGAc                                                                                                                                                                                                                                                    <  2:594500/250‑1 (MQ=255)
    aaTCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAC‑ggg                                                                                                                                                                                                                                                >  1:340465/1‑250 (MQ=255)
           ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTg                                                                                                                                                                                                                                         <  1:35666/251‑1 (MQ=255)
           ataataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTg                                                                                                                                                                                                                                         <  2:266950/251‑1 (MQ=255)
              ataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGGATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGcc                                                                                                                                                                                                                                      >  1:251315/3‑251 (MQ=255)
               taTGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGCACGGGCTGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCCCGCGCCCGTATTTCCGTGGTGCTGGTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGGACTTAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACTGGAACgcttgccc                                                                                                                                                                                                                                      >  1:254445/2‑244 (MQ=255)
                  ggCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCAtt                                                                                                                                                                                                                                   <  1:604702/250‑1 (MQ=255)
                    cAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGC‑GGAAAGGcggggcactttt                                                                                                                                                                                                                                >  2:237356/1‑240 (MQ=255)
                       tGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGTTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGAAGGAAGAGTTCTACCAGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGGC‑GGGACGGCTgggcacttatt                                                                                                                                                                                                                               >  2:570280/1‑240 (MQ=255)
                           cAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAGATGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGGACGGCAGGCCATTATCTCgaggg                                                                                                                                                                                                                         >  1:24906/1‑247 (MQ=255)
                            agcgtgtccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCggtggt                                                                                                                                                                                                                        <  2:302104/239‑1 (MQ=255)
                                   ccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGAt                                                                                                                                                                                                               <  1:123992/246‑1 (MQ=255)
                                     ggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGGCATTATTTCGGTGGG‑AGG‑GGATg‑g                                                                                                                                                                                                            >  1:374563/4‑251 (MQ=255)
                                       cccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGAC‑GGAACGGCTGGCCATTCACTCGGTGGT‑AGG‑GGATg‑gg                                                                                                                                                                                                           >  1:198460/2‑249 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTAt                                                                                                                                                                                                         <  1:3999/251‑1 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTAt                                                                                                                                                                                                         <  1:323048/251‑1 (MQ=255)
                                        ccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAC‑TGGACGGCTGGCCCTTATCTCGGTGGT‑AGG‑GGATG‑GTAt                                                                                                                                                                                                         >  1:319608/1‑251 (MQ=255)
                                                gAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTTTCTCGGGGGT‑AGG‑TGATG‑GTATGCGCACCt                                                                                                                                                                                                 >  1:606164/1‑251 (MQ=255)
                                                         gggtggCAGGGCGGCGCGCGTGGTTGCTGAGATATGACGCGTCCGTATTGCCGTGGTGCTGATTGCGCAATCATCTTCCGAATACAGTATCAGGTTCCGCGTTCCGCAAAGCGACTGGGTGCTAGCTGAACGGGCCATGCAGGAGGAGTTCTACTGGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGa                                                                                                                                                                                        <  2:162740/247‑1 (MQ=255)
                                                             ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGAtct                                                                                                                                                                                     <  1:570280/250‑1 (MQ=255)
                                                             ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCACGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGAtctc                                                                                                                                                                                    <  1:343346/251‑1 (MQ=255)
                                                                    cggcgCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑aaa                                                                                                                                                                            <  1:196652/251‑1 (MQ=255)
                                                                          gcgTGTTTGCAGCGATGTCACGCGCTCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACCGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCttt                                                                                                                                                                      <  1:643932/251‑1 (MQ=255)
                                                                              gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGGACGGATGGCCATTAACACGG‑GGTGAGGGTG‑TT‑GTATGCGCACCATGGGTGGGATTTCGGCGCAAATT‑TTTGccc                                                                                                                                                                  >  2:268561/2‑250 (MQ=255)
                                                                               tttGCAGCGATGTCACGCGCCCGTATGTCCTTGGTGGTGATTACGCAAGCGTCTACCGAATGCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCGGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCg                                                                                                                                                                  <  2:164807/250‑1 (MQ=255)
                                                                                             aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GTATGCGGAACTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGccc                                                                                                                                                   >  1:10537/1‑250 (MQ=255)
                                                                                                     gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAAc                                                                                                                                          <  1:150850/251‑1 (MQ=255)
                                                                                                         ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GGATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTg                                                                                                                                      >  1:93488/1‑251 (MQ=255)
                                                                                                           ccGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATGTGT‑TGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCACCATTGTc                                                                                                                                    >  2:110947/1‑251 (MQ=255)
                                                                                                             gtggtgCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCAACTTGCGTGGGGTCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTTTC‑Gc                                                                                                                                 >  1:655808/1‑251 (MQ=255)
                                                                                                                         aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAgg‑g                                                                                                                    <  1:8296/251‑1 (MQ=255)
                                                                                                                          cGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAgg‑g                                                                                                                    >  1:288402/1‑250 (MQ=255)
                                                                                                                          cGCAATCATCTGCCGAATACGGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑Ga                                                                                                                   <  2:340465/251‑1 (MQ=255)
                                                                                                                            cAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑GGATG‑GTCAGGGCACCCTTCGGGGGATCTCGGCG‑AAAATCTTTGCCGCCCTGGCCCCCCCCA‑AGAACCACATTGTCCGCC‑TTGCTCAgg‑ggt                                                                                                                  >  2:290569/1‑248 (MQ=255)
                                                                                                                                 atcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GGATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCCTTGCTCAGG‑GATCTTCTg                                                                                                            >  1:400836/1‑251 (MQ=255)
                                                                                                                                    ctCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCCACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAAc                                                                                                         >  1:193722/2‑251 (MQ=255)
                                                                                                                                     tCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGa                                                                                                           <  1:352537/248‑1 (MQ=255)
                                                                                                                                      ccGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACg                                                                                                        <  1:274611/250‑1 (MQ=255)
                                                                                                                                      cagaatgcttgatctgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGTGCTGAACGGGCACTGCAGTAAGTGTTCTTCCTGTATCTGACAGAATGCTTCCTGGTGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCGTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATACCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGc                                                                                                       <  2:514732/236‑1 (MQ=255)
                                                                                                                                         aataaagtaTCAGGTTCTGCGTTCGGCAAAGCGACGGTGTGCGAGCGGAACGGGCAATGCAGGAAGAGTTCTACCAGGAACTGAAAGAAGGCTTGCTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCa                                                                                                    <  2:251315/243‑1 (MQ=255)
                                                                                                                                         aaTACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGATCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑GGAGG‑GGAAGAGAACCCTGAGTGGGCTATCGGCG‑AAAATCTTTGCCGCGCTGGCCCGCGACA‑AAATAAACAATGTC‑GCACTTGCTCAGGAGA‑CTTCTGAAcccacc                                                                                                    >  2:439439/1‑246 (MQ=255)
                                                                                                                                            tcagtaTCAGATTCTGCGTTCCGCAAAGCGACTGTGTGAGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACAGCTGGCAATTATCTCGGTGGT‑AGG‑CGGTA‑GTAAGCGCACCTTGCGTGGGATATCGAGG‑AAATTCTTTGCAGGTCTGGCCCGC‑CAACATTTCAACCTTTTC‑GCCCTAGCACAGG‑GATCTTCTTAACGcccaca                                                                                                  >  2:510935/6‑245 (MQ=255)
                                                                                                                                            aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAAtc                                                                                                 <  1:237356/251‑1 (MQ=255)
                                                                                                                                             cagtatccgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GGTCTTCTGAACGCTCAAtct                                                                                                >  1:240711/9‑251 (MQ=255)
                                                                                                                                              agtaTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGGGGG‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCC‑ATATCCACATTGTC‑GCCATTGCTCCGG‑GATCTTCTGAACGCTCAAtctc                                                                                               >  1:375278/4‑251 (MQ=255)
                                                                                                                                               gtaTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGGACGGCTGGCCATTATCTCGGTGGT‑AGG‑CGATG‑GTATGCGCACCTTTCGGGGGGTCCCCGCG‑GAAATCTTTGGCGCGCTGGCCCCCGCCA‑ATATCAACGTTTTC‑GCCCTTGCTCCGG‑GGTCTTCTggctgaccaatctct                                                                                              >  2:202015/3‑237 (MQ=255)
                                                                                                                                                      tttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTgg                                                                                       >  1:420413/1‑251 (MQ=255)
                                                                                                                                                       ttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑CGATG‑GGATGCGCACACTGCGTAGGATCTCCGCG‑CAATTTTTTGCCCCGCAGGCCCGCGCCC‑ATATCAACATTTCC‑GCCCTTGCCCAGG‑GACCTTCTGAACGCACACTCTCCGTCGTggg                                                                                      >  2:253589/1‑250 (MQ=255)
                                                                                                                                                           gCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑Taaa                                                                                  <  1:202015/250‑1 (MQ=255)
                                                                                                                                                                 cgcaAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTACCTGGAACTGAAAGAAGGCTTACGGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACgat                                                                           <  1:233340/249‑1 (MQ=255)
                                                                                                                                                                        cGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGC‑GGAAAGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GGATGCGGACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTTTC‑GCCATTTCTCCGG‑GATCTTCTGAACGCTCAATATCTGTCGTGG‑TAAAAAACGATGATGCg                                                                     >  1:356202/1‑250 (MQ=255)
                                                                                                                                                                         gACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCATG‑GATCTTCTGAACGCTCAATCTCTGTCGTGGGTAAA‑AACGATGATGCGAc                                                                   >  1:109989/1‑251 (MQ=255)
                                                                                                                                                                         gACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGAc                                                                   >  1:348264/1‑251 (MQ=255)
                                                                                                                                                                          aCTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGAcc                                                                  <  1:357840/251‑1 (MQ=255)
                                                                                                                                                                                ggCGAGCTGAACGGGCAATGCAGGAAGATTTCTACCTGGAACGGAAAGAAGTCTGACTGGAGCCGCGGGCAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCGCGGCG‑AAATTCTTTGCCGCACTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGc                                                            <  2:422321/250‑1 (MQ=255)
                                                                                                                                                                                 gCGAGCTGAACGGGCAATGCAGGAAGAGTGCTACCTGGAACTGAAAGAAGGCTTACGGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGAGGATGCGACCACTGGc                                                            <  2:288402/250‑1 (MQ=255)
                                                                                                                                                                                   gAGCTGAACTGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCTATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTg                                                         <  1:96304/251‑1 (MQ=255)
                                                                                                                                                                                             ggCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGGTG‑GTATGCGCACCGTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCCTTGCTCAGG‑GATCTTCTGAACGCCCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTCGCGTGCGCGTTACtt                                               >  1:616517/1‑250 (MQ=255)
                                                                                                                                                                                                  ggCAGGAAGAGTTCTACCTTGAACGGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCACTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCTCCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                                           <  2:10537/249‑1 (MQ=255)
                                                                                                                                                                                                           agTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑GGATG‑GGATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTggtc                                 >  1:510813/1‑248 (MQ=255)
                                                                                                                                                                                                               atcccTGGAACTGACAGAAGGCTTATTGGAGGCGGTGGTGGGGAC‑GGAACGGGTGGGCTTTATCTCGGTGGT‑AGG‑GGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAAGTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACAGTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTTATCAGATGTTGTTCAAt                              <  2:606164/247‑1 (MQ=255)
                                                                                                                                                                                                                tACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGGGGT‑AGG‑GGGTG‑GTATGCGCAACTTGCGGGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GGTCTTTTGAACCCACCATCCCTGTCGGGG‑TAAATAACGCTAATGCGACCACTGGCGTGCGCGGTTCtcatccaggtcttgtcaatac                            >  2:156721/1‑235 (MQ=255)
                                                                                                                                                                                                                 aCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCGCGGTGGT‑AGG‑TGATG‑GTATTCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATAc                            <  1:428889/250‑1 (MQ=255)
                                                                                                                                                                                                                      gAACTGAAAGAAGGATTACTGGAGCCGCTGACAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGAAG‑GTATGCGCACCCTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCACTGGCCCGCTCCA‑ATATCAACATTGTA‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTAATCATCAGATGCTGTTCAATACCGATc                       >  1:345450/1‑250 (MQ=255)
                                                                                                                                                                                                                              agaagGCTTCTTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCg              <  2:93488/251‑1 (MQ=255)
                                                                                                                                                                                                                              agaagGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCg              <  1:13668/251‑1 (MQ=255)
                                                                                                                                                                                                                                 agGCTTACTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGG‑AGG‑GGATG‑GCATGCGCACCTTGGGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GGTCTTCTGGACGCCCAATCTCTGTCGTGG‑GAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGATGGTCAATACCGATAGGGTTATCGa             >  1:289873/1‑249 (MQ=255)
                                                                                                                                                                                                                                       aCTGGAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGttt      <  1:439439/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aCTGGAGCCGCTGGCGGGGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑CGG‑TGCTG‑GGATGCGCCCCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCC‑ATATCAACATTTTC‑GCCCTTGCTCAGG‑GTTCTTCTGAACGCCCACTCTCTGTCGGTG‑TAAATAACGCTGATGCGACCAATGGCGGGCGCGTTACTCATCAGGTGCTGTTCAATTCCGGTCCGGTTATCGGAGTttttt     >  2:343543/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gggtgCCGTTGGCGGTGAC‑GGAACGGCGGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTTGGATCTCGTCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGAGCAGGTTATCGAAGTGTTtgtg   <  2:109989/247‑1 (MQ=255)
                                                                                                                                                                                                                                          ggAGCCGCTGGCGGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCGCTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGa  <  1:87618/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGAC‑GGAACGGCTGGCCATTATCTCGGTGGT‑AGG‑TGATG‑GTATGCGCACCTTGCGTGGGATCTCGGCG‑AAATTCTTTGCCGCACTGGCCCGCGCCA‑ATATCAACATTGTC‑GCCATTGCTCAGG‑GATCTTCTGAACGCTCAATCTCTGTCGTGG‑TAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGA  >  NC_000913/1263‑1747

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