Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,797 G→A Q487Q (CAG→CAA)  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,7970GA98.8% 252.0 / NA 83Q487Q (CAG→CAA) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/49);  total (34/49)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  NC_000913/1558‑2040
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ccctttcgggggATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCGCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAAGACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTTAAGCGTCAACAAAgctggct                                                                                                                                                                                                                                           <  2:319608/242‑1 (MQ=255)
aCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGACCAGGTTCTCGAAGGGTTTGTGATTGGCGTCGGGGGCGGTGGCGGGGCGCTGCGGGAGCCACTGACGCGTAAACAAAgctggct                                                                                                                                                                                                                                           >  2:17787/1‑251 (MQ=255)
  cTTGCGTGGGATCTCGGCGAAATTCTTTGCAGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGAGTTACTCATCAGATGCTGTTCAATACCGATCAGGATATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAAATGAAGAGTCAACAAAgctggctg                                                                                                                                                                                                                                          >  1:504413/1‑250 (MQ=255)
     gCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaag                                                                                                                                                                                                                                       <  2:504413/250‑1 (MQ=255)
      cgggtgATCTCGGCGAAATTCTTTGCCGCGCGGGCCCGCGCCAATATCAACATCGTCGCCATTGCTCAGGGAGCTGCGGGACGCTCAATCTCTGTCGTGGTAAATAACGAAGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCTATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaagaa                                                                                                                                                                                                                                     <  2:198460/246‑1 (MQ=255)
       gTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAAt                                                                                                                                                                                                                                    >  1:471246/1‑251 (MQ=255)
          ggATCTCGGCTAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCGCAATCTCTGGCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATaaa                                                                                                                                                                                                                                 <  2:510813/251‑1 (MQ=255)
                cGGCGAAATGCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATGGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACata                                                                                                                                                                                                                             <  1:278449/249‑1 (MQ=255)
                  gCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGa                                                                                                                                                                                                                         <  1:268561/251‑1 (MQ=255)
                          cTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACgtg                                                                                                                                                                                                                 >  1:648993/1‑251 (MQ=255)
                             tGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCt                                                                                                                                                                                                              >  1:163702/1‑251 (MQ=255)
                             ggCCGCGCTGGCCCGCGACCATTTCAACATTGTCGCCATGGCTCAGGGATCTTCTTAACGCTCTATCGCTGGCGTGGTAAATAACGAGGATTCGACCACTGGCGTGCGCGTTACTCATCAGAGGCTGTTCAATGCCGATCAGGTTAGCGAAGTGTTTGTGAGTGTCGCCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAATCTTCAACAAAGCTGGCTGAAGAAGAAACATATCGACTTACgtggca                                                                                                                                                                                                              <  2:24906/250‑4 (MQ=255)
                              gCCGCGCTGGCCTGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCGCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCt                                                                                                                                                                                                              <  2:471246/250‑1 (MQ=255)
                                      ggCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAGCGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGgtgt                                                                                                                                                                                                       <  1:260233/249‑1 (MQ=255)
                                         ccGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGGCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCa                                                                                                                                                                                                  <  1:137679/251‑1 (MQ=255)
                                            cgcCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACt                                                                                                                                                                                               <  1:314583/251‑1 (MQ=255)
                                             gcCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTc                                                                                                                                                                                              <  1:354511/251‑1 (MQ=255)
                                              ccAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCg                                                                                                                                                                                             <  1:407077/251‑1 (MQ=255)
                                                aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGa                                                                                                                                                                                            >  1:125033/1‑250 (MQ=255)
                                                         aTTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCt                                                                                                                                                                                  <  2:648993/251‑1 (MQ=255)
                                                          gtggcgacatggcgtatggatctggtGAACGCTCAAACGCTGTCGTGGTAAATAACGATGATGCGACCACTGGTGTGCGCGTTTCTCATCAGATGCTGTTCAATCCCGATCAGGTTATCGATGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGCAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTc                                                                                                                                                                                 <  2:324984/226‑1 (MQ=255)
                                                              cGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGCTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCa                                                                                                                                                                             <  2:301528/251‑1 (MQ=255)
                                                                ccATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCAATGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTTGAGCAACTGAAGCGTCCACAAATCTGGCTGAAGACACAACATATCGACTTCAGTGTCAGCGGTGTTGcacaatcgtagacacggctaccaaat                                                                                                                                                                           >  2:278493/1‑226 (MQ=255)
                                                                 cATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑Tg                                                                                                                                                                         >  1:58181/1‑251 (MQ=255)
                                                                      cataGGGATCTTCGGAAAGCTCAATCTCTGTCGGGGTAAATAACGAAGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCAt                                                                                                                                                                    <  2:125033/248‑1 (MQ=255)
                                                                         aGGGATCTTCTGCACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCCGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGGGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgct                                                                                                                                                                                                                                               >  2:240974/1‑174 (MQ=255)
                                                                         aGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgct                                                                                                                                                                                                                                               <  1:240974/174‑1 (MQ=255)
                                                                          gggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGGAGCGTCCACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACCCGAAGGCAATACCCACCAA‑TGTGCATGGcc                                                                                                                                                                >  2:463057/1‑251 (MQ=255)
                                                                           ggATCTTCTGAAAGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGATATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCt                                                                                                                                                               >  1:112668/1‑251 (MQ=255)
                                                                                 tctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATc                                                                                                                                                          <  1:253589/250‑1 (MQ=255)
                                                                                    gAACGCTCAATCGCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGTCTAAATCTg                                                                                                                                                        <  1:107304/249‑1 (MQ=255)
                                                                                     aaCGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGa                                                                                                                                                      >  1:666637/1‑250 (MQ=255)
                                                                                      aCGCTCAATCTATGTCGTGGGAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATACTGGTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGAGACGGTGGCGTTGACGGTACGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGATGAAGAATAAACATATCGCCTTAAGTGTCTGCGGTGTTGACAACTCGACGGACCTGCTCACCAA‑Tgggcagggcctacaagtggta                                                                                                                                                     >  2:420680/1‑230 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                 <  1:227325/249‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTgg                                                                                                                                               <  1:256849/251‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTg                                                                                                                                                <  1:51793/250‑1 (MQ=255)
                                                                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCGGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTTTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                 <  2:374563/249‑1 (MQ=255)
                                                                                                   gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAAAATAAACATATCGACTTACGTGTCTGCGGGGTTTACAACTCGAAGGCACTGCTCACACA‑TGTGCATGGCCTAAATCTGGAAAACTGCCAGga                                                                                                                                          >  2:576554/1‑248 (MQ=255)
                                                                                                     cGTGGTAAATAACGATGATGCGACCACTGGCTGGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCGGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaag                                                                                                                                                                                                                                       <  2:429854/154‑1 (MQ=255)
                                                                                                     cGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaag                                                                                                                                                                                                                                       >  1:429854/1‑154 (MQ=255)
                                                                                                        ggTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTgg                                                                                                                                  <  1:626377/251‑1 (MQ=255)
                                                                                                             aTAACGATGAAGCGACCACTGGCATGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGGGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTATCTGCGGTGTTGCCAACCCGAAGTACCTGCTCACCAAATGGAC‑TGGACCAAATCTGGAAAACTGGCAGgaaaaaacggcccc                                                                                                                              >  2:444113/1‑239 (MQ=255)
                                                                                                               aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGGACTGGCGCAAGc                                                                                                                           <  1:501096/251‑1 (MQ=255)
                                                                                                                 cGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCa                                                                                                                         <  1:290569/251‑1 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                           >  2:3764/1‑212 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                           <  1:3764/212‑1 (MQ=255)
                                                                                                                      atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                                                                    <  1:268555/251‑1 (MQ=255)
                                                                                                                           aCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATTGCCTAAATCTGGAAAACTGGCCGGAAGAACTGGCGCAAGCCAAAGAGCCGtt                                                                                                               >  2:20377/1‑251 (MQ=255)
                                                                                                                              aCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTaa                                                                                                            <  1:20377/251‑1 (MQ=255)
                                                                                                                               cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACaaa                                                                                                                                                                                                                                                  >  2:600620/1‑117 (MQ=255)
                                                                                                                               cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACaaa                                                                                                                                                                                                                                                  <  1:600620/117‑1 (MQ=255)
                                                                                                                                 gtcgttctcgTTACTCTGCAGATGCTGTTCAAGTCCGATCAGGTTATCGACGTGTTTGTGATTGGCGTCGGTGTCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACTTGTCTTCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCGAAATCTTGAAAACTGGCAGGAAGAACTTGCGTAAGCCAAATAGCCGTTTAAt                                                                                                           <  2:34004/241‑1 (MQ=255)
                                                                                                                                    gTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                      <  1:249673/251‑1 (MQ=255)
                                                                                                                                          gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTAGCGTCGGTGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACAAATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGACCAAGCCAAAGAACCGTTTAATATCGGGCGCtt                                                                                                >  2:477947/1‑251 (MQ=255)
                                                                                                                                              cTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATCTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTAGGGCGCTTAAtt                                                                                            >  2:163745/1‑251 (MQ=255)
                                                                                                                                                catctGATGCTGGTCAATTCCGATCAGGTTATCGAAGTGTTTGGGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAATAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGTACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCTATAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                           <  2:655808/250‑1 (MQ=255)
                                                                                                                                                catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGTTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                           <  1:388932/250‑1 (MQ=255)
                                                                                                                                                catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                          <  1:247855/251‑1 (MQ=255)
                                                                                                                                                catcaGAGGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                          <  1:617089/251‑1 (MQ=255)
                                                                                                                                                    aGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTc                                                                                      <  1:609608/251‑1 (MQ=255)
                                                                                                                                                     ggtgtggTTCAATACCGAGCAGGTTTTCGAAGTGGTTGTGATTGGCTGCGGTGGCGTGGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCCAAAGAGCTGTTTAATCTCTGGCGCTTAATTCGCCTCg                                                                                     <  2:400836/245‑1 (MQ=255)
                                                                                                                                                                  aCCGATCAGGTTCTCGTAGTGTTTGTGATTCACGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAGAGCTGGCTGGGGAATAACCATATCGGCTTACGTGTCTGCGGTGTTGCCAACTCGCAGGCACTGCTCCCCAC‑TGAGCATGGCCTAAATCTGGaa                                                                                                                                                     >  2:73224/1‑174 (MQ=255)
                                                                                                                                                                  aCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGaa                                                                                                                                                     <  1:73224/174‑1 (MQ=255)
                                                                                                                                                                         aGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgctg                                                                 >  1:41477/1‑251 (MQ=255)
                                                                                                                                                                            ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGCATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACAAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGAGCAAGCAAAAGCGCAGTTTTATATCGGGCGCTTTATTCGCCTCGTGAAAGAATATCATCTGATGAAc                                                              >  2:29525/1‑251 (MQ=255)
                                                                                                                                                                            ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                                                              >  1:439176/1‑251 (MQ=255)
                                                                                                                                                                             tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                           <  1:125305/205‑1 (MQ=255)
                                                                                                                                                                             tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                           >  2:125305/1‑205 (MQ=255)
                                                                                                                                                                              aTCGAAGTGTTTGGGATTGGCGGCGGGGGAGTAGGCGGTGCGCTTCTGGTGCAACTGAATCGTCAACAAAGCGTGCTGATGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCATTGCTCACCAA‑TGTGCATGGCCTAAATTTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                             <  2:375278/250‑1 (MQ=255)
                                                                                                                                                                                cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                           <  1:156721/250‑1 (MQ=255)
                                                                                                                                                                                     tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCCAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGCAAGAATATCATCTGCTGAACCCGGTCAtt                                                     >  2:563405/1‑251 (MQ=255)
                                                                                                                                                                                      gtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttg                                                    <  1:545447/251‑1 (MQ=255)
                                                                                                                                                                                            gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt                                              >  2:176529/1‑251 (MQ=255)
                                                                                                                                                                                            gATTGGAGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCCCCCCGTGAAAAAATATCATCTTCTGACCCCGGTCATTGTTGACt                                              >  2:280962/1‑251 (MQ=255)
                                                                                                                                                                                             aTTGGTGTCTGTGGCGTTGGGGGTTCGTTGATGGAGCAACTGAGGCGTCAATAACGCTGGCTGAAGAAGAAACATATCGTCTGACGTGTCTGCTGTGTTGCCAACTCGAAGGCACAGCACACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAAATAGCCGTTTAATCACGGGCGCGTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGGCATTGTTGACt                                              <  2:543366/250‑1 (MQ=255)
                                                                                                                                                                                                              tGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGACTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGCACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGccg                            >  2:261741/1‑249 (MQ=255)
                                                                                                                                                                                                                            tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGCAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGt                                                                                                                                                                        >  2:352448/1‑97 (MQ=255)
                                                                                                                                                                                                                            tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGt                                                                                                                                                                        <  1:352448/97‑1 (MQ=255)
                                                                                                                                                                                                                               aGCAACTGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc           <  1:611862/251‑1 (MQ=255)
                                                                                                                                                                                                                               aGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc           <  1:344392/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATgacgg         >  1:89258/1‑247 (MQ=255)
                                                                                                                                                                                                                                     tGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGACCAATATGCCgccctc     >  2:673986/1‑246 (MQ=255)
                                                                                                                                                                                                                                        aGCGTCAACAAATATGGCTGAAGAATAATCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAACGAGACGTTTAATCTCGGGCGCTTAATTAGCCTCGTGAACGAATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGGTCAATATGCCGACtttctt  >  2:378470/1‑247 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  NC_000913/1558‑2040

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